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Abstract 



In this thesis we first apply the 1+3 co variant description of general relativity to analyze n-fluid Friedmann- 
Lemaitre (FL) cosmologies; that is, homogeneous and isotropic cosmologies whose matter-energy content 
consists of n non-interacting fluids. We are motivated to study FL models of this type as observations sug- 
gest the physical universe is closely described by a FL model with a matter content consisting of radiation, 
dust and a cosmological constant. Secondly, we use the 1 + 3 covariant description to analyse scalar, vec- 
tor and tensor perturbations of FL cosmologies containing a perfect fluid and a cosmological constant. In 
particular, we provide a thorough discussion of the behaviour of perturbations in the physically interesting 
cases of a dust or radiation background. 



iii 



Acknowledgements 

First and foremost, I would like to thank Dr. John Wainwright for his patience, guidance and support 
throughout the preparation of this thesis. I also would like to thank Dr. Achim Kempf and Dr. C. G. 
Hewitt for their time and comments. 



iv 



Contents 



1 Introduction 

1 . 1 Cosmological Models 

1.2 Geometrical Background 13 

1.2.1 Projection and Differentiation 4 

1.2.2 Kinematic Quantities 5 

1.2.3 Decomposition of the Energy-Momentum Tensor 7 

1.2.4 The Weyl Curvature Tensor 

1.2.5 Spatial Curvature 

1.2.6 Evolution and Constraint Equations [10 

1.3 Overview of the Thesis [jj 

2 The Friedmann-Lemaitre Cosmologies 13 

2.1 Properties of the RW Metric 13 

2.2 Characterizations of the FL Cosmologies 17 

2.2.1 The Kinematic Characterization 17 

2.2.2 The Weyl Characterization [18 

3 Dynamics of re-Fluid FL Cosmologies 2C 

3.1 Behaviour and Classification of n-Fluid FL Cosmologies 20 

3.1.1 Basic Evolution Equations 21 

3.1.2 n-Fluid FL Cosmologies |23 



v 



3.1.3 Hubble-Normalized Scalars 

3.1.4 A Compact State Space for n-fluid FL Cosmologies 

3.1.5 Invariant Sets and Equilibrium Points in the Compact State Space 

3. 1.6 Classification of n-Fluid FL Cosmologies 

3.1.7 Asymptotic Behaviour of Solutions 

3.2 Parameterization of FL Models 

3.2.1 The Intrinsic Parameters 

3.2.2 The Observational Parameters 

3.2.3 Conserved Quantities 

3.3 2-Fluid FL Cosmologies - Qualitative Analysis 

3.3.1 General Features 

3.3.2 FL Cosmologies with Radiation and Dust (i?Z?C-universes) . . . 

3.3.3 FL Cosmologies with Dust and Cosmological Constant 

3.4 3-Fluid FL Cosmologies - Qualitative Analysis 

3.4.1 General Features 

3.4.2 A Model of the Physical Universe 

3.5 Explicit FL Cosmologies 

3.5.1 Single Fluid FL Cosmologies 

3.5.2 FL Vacuum Cosmologies 

3.5.3 Flat FL with a Single Fluid and A 

3.5.4 FL Cosmologies with Radiation and Dust (i?DC-Universes) . . . 

3.5.5 FL Cosmologies with Radiation and A (i?CA-Universes) 

3.5.6 Explicit Solutions in Terms of Observational Parameters 

3.6 Discussion 

4 Perturbations of FL Models 

4. 1 Historical Development 

4.1.1 Metric Approach 

4.1.2 Geometrical Approach 

4.2 The Linearized Einstein Field Equations 

4.2.1 The Linearization Process 

4.2.2 The Linearized Evolution and Constraint Equations 



vi 



4.3 Harmonic Decomposition 

4.4 Scalar Perturbations 

4.4.1 The Governing DE 

4.4.2 Kinematic Quantities and Weyl Curvature 

4.5 Vector Perturbations 

4.5.1 The Governing DE 

4.5.2 Kinematic Quantities and Weyl Curvature 

4.6 Tensor Perturbations 

4.6.1 The Governing DE 

4.6.2 Kinematic Quantities and Weyl Curvature 

4.7 Dynamics of the Linear Perturbations: General Features 

4.7.1 Conformal Time and the Particle Horizon 

4.7.2 Governing DEs Using Conformal Time 

4.7.3 General Solutions for Scalar and Tensor Perturbations with A = 

4.7.4 Functional Dependence of First-Order Quantities 

4.8 Perturbations of FL Models with Pressure-Free Matter 

4.8.1 The General Solution in Integral Form for Scalar Perturbations . . 

4.8.2 Scalar Perturbations with A = 

4.8.3 Scalar Perturbations with A > 0, K = 

4.8.4 Tensor Perturbations for A = 

4.9 Perturbations of Radiation-Filled FL Models 

4.10 Perturbations in Asymptotic Epochs of FL Models 

4.10.1 The Flat FL Asymptotic Epoch 

4.10.2 The Milne Asymptotic Epoch 

4.10.3 The de Sitter Asymptotic Epoch 

4.11 Discussion 

A Evolution and Constraint Equations 

B Spatial Curvature 

B.l The Generalized Gauss Equation 

B.2 The Spatial Gradient of the 3-Ricci Scalar 



vii 



C Solutions of the Friedmann Equation 

C.l Solutions with a Stiff Fluid 

C.2 Existence-Uniqueness for the Friedmann Equation 

D Conformal Time 

E Covariant Differential Identities 

F Covariant Harmonics 

F 1 Scalar Harmonics 

F2 Vector Harmonics 

F3 Tensor Harmonics 

G The 3-Cotton-York Tensor 

G. 1 Evolution Equation for the 3-Cotton-York Tensor 
G2 Tensor Perturbations 



H Bessel and Associated Legendre Functions 

H. 1 Normal Forms of Bessel and Associated Legendre DEs 

H.2 Special Solutions of the Governing DEs 

H.3 Resonant Solutions of the Legendre DE 

H.4 Asymptotic Behaviour 

Bibliography 



viii 



List of Figures 



3. 1 The 2-fluid FL state space for i?DC-universes 

3.2 The 2-fluid FL state space for DCA-universes 

3.3 Thumbnail plots depicting behaviour of the length scale in FL cosmologies 

3.4 The 3-fluid FL state space for radiation, dust and cosmological constant 

3.5 3-fluid FL state space cross-section along the expanding flat FL manifold 

3.6 Past and future evolution of the RDCA model of the physical universe 

4. 1 Conformal time and the particle horizon 

4.2 Scalar perturbation modes for dust and spatial curvature 

4.3 Scalar perturbation modes for dust and cosmological constant 

4.4 Tensor modes of the Hubble-normalized electric Weyl tensor in open FL with dust . . . 

4.5 Tensor modes of the Hubble-normalized electric Weyl tensor in open FL with radiation 



105 



107 



109 



111 



ix 



List of Tables 



3. 1 The past and future attractors of generic n-fluid FL cosmologies 42 

3.2 Asymptotic behaviour of Hubble-normalized quantities in generic n-fluid FL cosmologies. |42 



4. 1 Functional dependence of first-order quantities 

4.2 Perturbations of flat FL with kr] <C 1 

4.3 Perturbations of flat FL, with kr] S> 1 .... 

4.4 Perturbations of Milne 

4.5 Perturbations of de Sitter (kfj Cl) 

4.6 Perturbations of de Sitter (kfj 3> 1) 



x 



Chapter 1 



Introduction 



1.1 Cosmological Models 

The goal of cosmology is to describe the large-scale structure, origins and evolution of the universe as a 
whole. Within the framework of general relativity, a cosmological model, denoted (M, g, u) consists of 
a differential manifold M , a Lorentzian metric g (which describes the geometry of the space-time) and a 
unit timelike vector field u, 

#juV = -1, (1.1) 

that defines the world-lines of a family of fundamental observers. In cosmology one represents the matter 
content as a smooth distribution (a continuum) whose particles represent galaxies in the present epoch. 
The time-like vector field u is the 4-velocity of the particles in this continuum. One member of the 
congruence of time-like curves defined by u, namely the world-line of our galaxy, plays a particularly 
important role, since it is the position from which we make observations. 

The evolution of the universe is assumed to be governed by the Einstein field equations (EFEs), 

R a b - \Rg a b = T a b, (1.2) 

where R a i is the Ricci tensor, R is the Ricci scalar (R = g ab R a b) an d T a b is the energy -momentum tensor 
of the distribution of matter and energy in the universe. Here and throughout we will use geometerized 
units so that c = 1 and 8irG = 1, where c is the speed of light in vacuum and G is the gravitational 
constant. 

There are two fundamental assumptions that are used to construct models of the universe: 
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I. The universe is isotropic about our galaxy. More precisely, at any event on the world-line of our 
galaxy, all directions orthogonal to the world-line are physically and geometrically equivalent, when 
the universe is viewed on a sufficiently large scale. 

A variety of cosmological observations {e.g. number counts of distant galaxies and peculiar velocity 
surveys) provide support for assumption I in the present epoch. The strongest support is provided 

by the temperature of the cosmic microwave background radiation (CMBR) which is isotropic to 

AT 
T 



within lfT 5 (i.e. 4£ ~ 1(T 5 ). 



II. The Copernican Principle. This principle states that our galaxy is not in a preferred position in the 
universe. Unlike assumption I, this assumption is philosophical, rather than observational. 

Assumptions I and II together imply that the universe satisfies the so-called Cosmological Principle, 
namely 

i) the universe is isotropic about all fundamental world-lines, and 

ii) the universe is spatially homogeneous, i.e. all physical and geometric quantities do not vary in 
space. 



It is a well-known consequence of assumptions I and II (we give a concise proof in Chapter |2j see 
Theorem |2.1.1| ) that the metric g a j, can be written in the form of the Robertson-Walker (RW) metric]. 



l 



ds 2 = -dt 2 + e 2 (t)g al3 (x X )dx a dx 13 , (1.3) 

where g a p is a 3-metric of constant curvature. The preferred timelike vector field u is given by 

d 

u =a? (L4) 

The Robertson- Walker metric provides the basis for the Friedmann-Lemaitre (FL) cosmological mod- 
els that play a central role in modern cosmology. A Friedmann-Lemaitre cosmolog]^ is a cosmological 
model (M., g, u) that satisfies the Einstein field equations (11.21) with a suitable energy-momentum tensor 
T a b, with metric g and fundamental congruence u having the Robertson- Walker form (11.31 ) and (11.4I ). 

In this thesis we study FL cosmologies and their perturbations within a general mathematical frame- 
work. In section [Ol we provide the necessary geometrical background for this study. 



Here we use Greek indices to denote spatial indices. They can take on the values 1, 2 or 3 
2 Tht 
(1989)) 



2 The first explicit FL cosmological models were first discovered by Friedmann in 1929 (for some historical details, see Ellis 
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1.2 Geometrical Background 



In this thesis we will make use of the so-called 1 + 3 covariant description of general relativity^ which 
has been developed for use in spacetimes in which there is a preferred timelike congruence u, as in the 
case of a cosmological model {M., g, u). The "1+3" refers to the fact that one performs a "time + space" 
decomposition relative to u by projecting tensors and tensorial equations parallel to u and orthogonal to u. 
The second aspect of the 1 + 3 description is to write a tensor as a sum of algebraically simpler parts, i.e. 
to give an algebraic decomposition. For example, an arbitrary rank 2 contravariant tensor can be written 
as the sum of a symmetric tracefree tensor, an anti-symmetric tensor and the trace multiplied by the metric 
tensor. 

The first step is to introduce the basic variables in the 1 + 3 description, which are defined by ap- 
plying projection relative to u and algebraically decomposing the covariant derivative VbU a , the energy- 
momentum tensor T a b and the Riemann-Christoffel tensor as follows: 

i) The covariant derivative of the fundamental 4-velocity, VbU a , describes the kinematics of the con- 
tinuum described by the vector field u. The so-called kinematic quantities are defined by writing 
VbU a as a sum of simpler parts and projecting relative to u. We define these quantities in section 



ii) The energy-momentum tensor T a b describes the matter and energy content of the universe and so 
gives rise to a set of quantities known as matter variables when decomposed relative to u. We 
discuss this decomposition in section 11.2.31 

iii) The Weyl curvature tensor C a b c d, which is defined as the tracefree part of the Riemann curvature 



tensor, is decomposed relative to u into two tracefree symmetric tensors, the so-calledj electric part 
E a b and magnetic part H a b. We discuss this decomposition in section [T. 2 .4 1 

The second step is to derive evolution and constraint equations for the basic variables by using the 
Ricci identities and Bianchi identities in conjunction with the Einstein field equations. We will discuss 
this process in section [T~.2.6i 



OS 



3 See, for example, Ellis (1973), Wainwright and Ellis (1997), pl-30, or Ellis and van Elst (1998), Chapter 2. 
4 There is a formal analogy with the decomposition of the electromagnetic field tensor F a t into an electric field E a and a 
magnetic field H a - 
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1.2.1 Projection and Differentiation 

Projection parallel to u is achieved by simply contracting a vector or tensor with u a . Projection into the 
3-space orthogonal to u is defined by the projection tensor 



which satisfies the identities 



Kb = 9ab + U a U b , 



h a \ h a b u b = 0, K = 3- 



(1.5) 



Corresponding to these two types of projection, two derivatives are defined. 

The overdot ' denotes differentiation along the fundamental congruence, defined for scalars as 



and for tensors as 



/ = f,aU a , 



Tab-c — U d VdT a b... c , 



(1.6) 



(1.7) 



where V a denotes covariant differentiation. 

Following Maartens (1997), we use angle brackets to denote the orthogonal projections of vectors, the 
orthogonally projected symmetric trace-free part of tensors, and their time derivatives, i.e. 



X {a) = h a b X b , 
X {a) = h a b X\ 



rj-i(ab 



hth b \ - \h ah h cd 



c"- d 3 

= \h^ a c h b) , - \h ab h cd 



'icd 



(1.8) 
(1.9) 



Further, we will use the projection symbol _L to identify orthogonal projection over all tensor indices, i.e. 



a ub ur us 



(1.10) 



The spatial covariant derivative V a is defined by projecting all free indices of the covariant derivative 
(compare Ellis (1973), Ellis and Bruni (1989)). For example, 



Vo/ = K b V b f, V a T bc de = Kh b h r c hih^V p T q \ 



sf 



(1.11) 



The covariant derivative of a scalar can then be decomposed into a derivative along the fundamental 
congruence u and a spatial covariant derivative via 



V a / = V a / - U a f. 



(1.12) 
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The totally skew pseudotensor is defined by 

^abcd = ^abcd] ^ ^0123 = (_^-l/2 9 (L13) 

where 

g = det{g ab ). (1.14) 
The projected permutation tensor is then defined by (compare Maartens et al. (1995)), 

^abc = riabcdU d ■ (1.15) 

and is used to define the spatial curl of a vector by 

curl X a = e afec V b X c . (1.16) 

and the spatial curl of a rank 2 tensor by 
1.2.2 Kinematic Quantities 



curl T ab = e cd(a V c T b) d . (1.17) 



We now focus on the decomposition of the covariant derivative of the fundamental 4-velocity. Following 
Ellis (1973, p8-l 1), we decompose the covariant derivative of a unit timelike vector field u as 

V b u a = a ab + u ab + Hh ab -ii a u b) (1-18) 

with a ab symmetric and trace-free, ui ab antisymmetric, and 

a ab u a = 0, LO ab u a = 0, u a u a = 0. (1.19) 

It follows from (fTTT8T) and (fTTT9l that 

H = lV a u a , (1.20) 

it a = u b V b u a , (1.21) 

Oab = V( a n b) - Hh ab + ii( a u b ), (1-22) 

Uab = V[6« tt ] +U[ a U b ]. (1.23) 



1.2. Geometrical Background 



6 



Here, H is the Hubble scala^, u a is the acceleration vector, a a b is the rate of shear tensor and uj a b is the 
vorticity tensor. In practice we will usually replace uj a t, by the vorticity vector, defined by 

w a = l^bcd^^ (L24) 

which satisfies u a cu a = 0. We can write the vorticity tensor in terms of the vorticity vector using the 
relation 

u ab = ^bcd^^ (1 25) 

The magnitude of the rate of shear tensor and the vorticity tensor are defined by 

2 1 ab 2 a 

We now state two useful propositions concerning an irrotational vector field u a , i.e. u a = 0. 

Proposition 1.2.1 Given a vector field u, there exist scalar fields f(x a ) and t(x a ) such tha$u a = —ft a 
if and only ifiv a = 0. 

Proof: Given u a = —ft >a , it follows by differentiating that Ui aib u c i = and that u a = 0, by d 1.24b . The 
converse is a special case of Frobenius' Theorem, which is usually stated in terms of differential forms 
(see for example Flanders (1989), see the theorem on page 94). □ 



Comment: The significance of this result is that if uj = 0, then the 3-spaces orthogonal to u a at each 
point form spacelike hypersurfaces given by t = constant, i.e. there is a one-parameter family of space- 
like hypersurfaces orthogonal to the timelike congruence defined by u. 

If we can further impose the condition of zero acceleration (u a = 0), we obtain the following result. 

Proposition 1.2.2 u> a = 0, ii a = <^=^ u a = — t iQ , for some scalar field t(x a ). 



5 The rate of expansion scalar B, defined by O = 3H, is often used in the expansion i 1.1 8b . 

6 The minus sign ensures that i > 0, i.e. t increases into the future along the fundamental congruence. 
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Proof: Given uj a = ii a = 0, equation dl.231 ) implies u\ a . b ] = U[ a ,6] = 0. It then follows from Proposition 
ll.2.1l that f \ a t h = 0. This in turn implies that / = g(t) and hence that one can redefine t to set / = 1. □ 

In this case, the scalar t that labels the hypersurfaces is clock time along the normal congruence. 

1.2.3 Decomposition of the Energy-Momentum Tensor 

The energy-momentum tensor T ab is a symmetric tensor that specifies the total matter-energy content of a 
cosmological model. We use the standard decomposition of the energy-momentum tensor with respect to 
a timelike vector field u (see Ellis (1973), p7): 

Tab = ^u a u b + q a u b + q b u a + ph ab + 7r ab , (1.26) 

where 

q a u a = 0, n ab u b = vr a a = 0, n ab = ir ba . (1.27) 

In this equation, \i is the total energy density measured by an observer moving with fluid 4- velocity u a , 
q a is the energy flux relative to u a (representing heat conduction and diffusion), p is the isotropic pressure 
and ir ab is the trace-free anisotropic pressure (due to processes such as viscosity). 
A perfect fluid with 4- velocity u is defined by 

q a = 0, TT ab = 0, (1.28) 

in which case the energy-momentum tensor (11.261 ) simplifies to 

T ab = li.u a u b + p(g ab + u a u b ). (1.29) 

1.2.4 The Weyl Curvature Tensor 

The Weyl conformal curvature tensor (or simply the Weyl tensor) is defined in terms of the Riemann tensor 
and the Ricci tensor according to 

Cabcd = Rabcd ~ \ {9acRbd + 9bdRac ~ QbcRad ~ QadRbc) + \R{9ac9bd ~ 9ad9bc)- (1.30) 

(see, for example, Ellis 1973, p7). We note that the Weyl curvature tensor has all the symmetry properties 
of the Riemann tensor and is in addition trace-free (C a b d = 0). 
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It is helpful to view dl.301 ) as a decomposition of the Riemann tensor R a b c d, with its 20 independent 
components, into a tracefree part C a b c d and the trace R a b'. 

Rabcd < ► {Cabcd, Rab}, (1-31) 

each of which have 10 independent components. The Ricci tensor is determined at each point by the 
energy-momentum tensor through the field equations (1 1 - 2b - On the other hand, as we shall see later, the 
Weyl tensor governs the propagation of changes in the gravitational field, i.e. gravitational waves. 

Within the 1 + 3 formalism the Weyl tensor can be decomposed further into an "electric part" E a b and 
a "magnetic part" H a b by projecting C a b c d and its dual orthogonal to u: 

E ac = C abcd U b U d , H ac = hab fC efcdU b U d . (1.32) 

We note that E a b and H a b are symmetric, trace-free and orthogonal to u: 

E ab u h = 0, H ab u b = 0. (1.33) 
In addition, this decomposition describes the Weyl tensor completely sincd^ 

E ab = = H ab <=> C abcd = 0. (1 .34) 

Note that through (11.301 ) E a b, H a b and R a b provide a full description of the space-time curvature as 
described by the Riemann tensor. 

1.2.5 Spatial Curvature 

If u is irrotational, we have seen that there exists a family of hypersurfaces (i.e. 3-dimensional manifolds) 
orthogonal to u. In this situation the projection tensor h a b can be viewed as the 3 -metric induced on these 
hypersurfaces and the spatial covariant derivative V a can be viewed as the covariant derivative determined 
by this 3-metric. 

Analogous to the spacetime Ricci identities, the 3-space Ricci identities are given by 

V c V d X a - V d V c X a = ^R abcd X b , (1.35) 



7 The Weyl tensor can, in fact, be written in terms of E a b and H a b (see, for example, Ellis and van Elst (1998) eq. (22)). 
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where ^'R a bcd * s tne spatial Riemann tensor and X a is any vector field in the hypersurfaces, i.e. X a u a = 
0. It can be shown by applying (11.351) to an arbitrary vector field X that is orthogonal to u, expanding the 
left hand side using d 1 - 1 lb and d 1 - 51 ) and introducing the spacetime Riemann tensor via 

V c V d X a - V d V c X a = R abcd X b , (1.36) 

that ^Rabcd is related to R abcd via 

^ Robed = (Rabcd)± ~ @ac<3>bd + @ad@bc, (1-37) 

where Q ab is the expansion tensor of the congruence u, defined by 

<5>ab = O-ab + Hh ab . (1.38) 

Equation dl.371 ) is the well-known Gauss equation (see, for example, Ellis (1973), p33-4). 
The spatial Ricci tensor and spatial Ricci scalar are then defined in the usual manner: 

®R ab = h cd ^R acbd , MR = h ab ^R ab . (1.39) 

Further, the trace-free spatial Ricci tensor is defined as 

^S a b = (3) Rab ~ I {3) Rh ab . (1-40) 

We can obtain an expression for the spatial Ricci scalar and trace-free spatial Ricci tensor in terms of the 
spacetime quantities from (1B.5I ). incorporating the Weyl tensor through (11.301 ) and (11.321 ): 

(3) i? = R + 2R bd u b u d -6H 2 + 2a 2 , (1.41) 
{3) S a b = E ab + \R {ab) -H<j ab + a c {a <j b)c . (1.42) 

Then using the field equations (11.21 ) in the form (IA.2b - (IA.5b to eliminate R and R ab , we obtain 

= -QH 2 + 2fi + 2a 2 , (1.43) 
(%o6 = E ab + \-K ab - Ha ab + a c {a a b)c . (1.44) 

We conclude this section with a brief digression on spaces of constant curvature. A differentiable 
manifold M. of dimension n with metric g is said to be a space of constant curvature if and only if 

Rabcd = C(g ac g bd - g a d9bc)- (1-45) 
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It follows by contraction that 

R = n(n-l)C, and R ab - ±Rg ab = 0, (1.46) 

i.e. the tracefree Ricci tensor S ab = R ab — \Rg ab 1S zero. In addition it follows from the twice contracted 
Bianchi identities that if n > 2 then R is constant (see proposition [L23]). If n = 3, the Riemann tensor is 
completely determined by the Ricci tensor according tqj 

Rabcd = 4/i[ rf [ b S a ] c ] + \Rg a [c9b}d- (1-47) 

It follows that ifn = 3 then (A4, g) is a space of constant curvature if and only if S ab = 0. Within the 
framework of the 1 + 3 covariant formalism this result can be stated as follows. 

Proposition 1.2.3 Let u be an irrotational timelike vector field. The hypersurfaces orthogonal to u are 
spaces of constant curvature if and only if 

{3) S ab = 0. (1.48) 

Further, if this condition holds, then 

V a (3) i? = 0. (1.49) 
1.2.6 Evolution and Constraint Equations 

In the 1 + 3 covariant formalism, the information contained in the EFEs is displayed in an indirect manner 
by using the Ricci identities applied to the fundamental 4-velocity u, 

V c V d u a - V d V c n a = R abc dU b , (1.50) 

and the Bianchi identities, 

V [a R bc]de = 0. (1.51) 

The latter can be written in divergence form by contracting once and then twice and using (11.301 ) to 
introduce the Weyl tensor§| 

V d C abc d = -V [a R c b] - ±5 c [a V b] R, (1.52) 



8 See, for example, Ellis (1973) eq. (80). 
9 Also see Wainwright and Ellis (1997), p27. 
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and 

V b {R a b - \R8*) = 0. (1.53) 

The Riemann tensor in (11.501 ) is likewise expressed in terms of the Weyl tensor and Ricci tensor using 
(11.30I ). In all these identities, the Ricci tensor is expressed in terms of the energy-momentum tensor (11.261 ) 
using the EFEs in the form 

Rab = T a b — \g a bT c c , (1-54) 

which introduces the fluid 4-velocity u into d 1.521 ) and d 1.531 ). Equation dl.181 ) is used to replace all 
occurrences of V u a by the kinematic quantities and d 1 .321 ) is used to replace C a b c d by E ab and H^. 
Finally, one projects all three identities parallel to and orthogonal to u. The result is a set of evolution 
equations giving the derivatives of H, a a t, u> a , E ao + ^ir a b, H a b, H and q a along the congruence u, and 
a set of constraint equations that contain only the spatial covariant derivative V a . We refer the reader to 
Appendix E] for the complete list of equations that arise from this process. 

1.3 Overview of the Thesis 

The main goal of this thesis is twofold. Firstly, we analyze the dynamics of so-called n-fluid FL cos- 
mologies, that is, homogeneous and isotropic cosmologies whose matter-energy content consists of n 
non-interacting fluids. Secondly, we give a unified coordinate-independent discussion of perturbations 
of FL cosmologies containing a perfect fluid and a cosmological constant. Some parts of this thesis are 
new, while other parts provide a unified presentation of known results for these cosmologies. Our treat- 
ment throughout is completely coordinate-independent since it is based on the so-called 1 + 3 covariant 
description of general relativity and in this respect it differs from much of the current literature. 

In Chapter 1 , we introduce the 1 + 3 covariant description of general relativity and in Appendix [A] 
we give the full system of evolution equations and constraint equations. In this sense the thesis is self- 
contained. However, a reader who has no prior knowledge of the 1+3 covariant description will encounter 
an initially steep learning curve. We assume the reader has a basic knowledge of general relativity and 
some familiarity with dynamical systems theory (see, for example, Wainwright and Ellis (1997), Chapter 
4). 

Chapter 2 plays a supporting role. In this chapter we describe the geometrical, i.e. coordinate indepen- 
dent properties of the Robertson- Walker metric, and give two characterizations of the Friedmann-Lemaitre 
models that highlight their very special nature. 
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In Chapter 3, we give a complete qualitative description of the dynamics of n-fluid FL cosmologies. 
A central role is played by a new formulation of the evolution equations for these models as a dynamical 
system on a compact n-dimensional state space. The most important example is the 3-fluid model in 
which the three fluids are radiation, pressure-free matter and a cosmological constant, since it provides a 
model of the physical universe. Our discussion of this model is one of the highlights of this chapter. 

In Chapter 4, we give a unified derivation of the governing equations for scalar, vector and tensor per- 
turbations of FL cosmologies containing a perfect fluid and a cosmological constant, using the so-called 
geometrical approach pioneered by Ellis et al (1989). We then analyze the behaviour of the perturbations 
in various regimes using both exact and approximate solutions. In particular we determine how the pertur- 
bations affect the so-called observational quantities that determine the extent to which a perturbed model 
deviates from an exact FL model. 

The Appendices contain a wealth of additional material relating to perturbations of FL universes, not 
all of which is needed for the purposes of this thesis. 



The Friedmann-Lemaitre 
Cosmologies 



In this chapter we first derive the geometrical (i.e. coordinate independent) properties of the Robertson- 
Walker metric within the framework of the 1 + 3 description of general relativity. We then give two 
characterizations of the Friedmann-Lemaitre cosmologies, firstly using the kinematic quantities and sec- 
ondly using the Weyl tensor. 

2.1 Properties of the RW Metric 

In this section we show that the high symmetry of the RW line-element leads to strong restrictions on the 
kinematic quantities of the fundamental congruence, on the Weyl tensor and on the Ricci tensor. We then 
derive the RW line-element assuming that the Cosmological Principle holds. 

Proposition 2.1.1 The kinematic quantities of the fundamental congruence &1.4\) associated with the 
Robertson-Walker metric M.3\) satisfy 

o- ab = 0, uj a = 0, u a = 0, y a H = 0. (2.1) 
Proof. For the RW line-element (1 1 - 3b . the co variant derivative 

V b U a = Ma,6 - T C ah U c (2.2) 
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simplifies to 

V 6 *u a = \g a b,o- (2-3) 



Furthermore, it follows from the RW metric (11.31) that 

£ 

9a/3,o = 2 j9a(3, and g a0)0 = 0. (2.4) 
Since the projection tensor h ab has components h a p = g a p, h a o = 0, it follows that 

£ 

V 6 n a = -h ab . (2.5) 
The decomposition dl.181 ) immediately implies a ab = 0, u a b = 0, u a = and 

H = ^. (2.6) 
Thus, is purely a function of t, which implies V a H = 0. □ 

We now consider the restrictions on the spatial curvature associated with the Robertson- Walker metric. 
Proposition 2.1.2 The spatial curvature of the RW metric satisfies 

(% a6 = 0, V a ( 3 )i? = 0. (2.7) 

Proof: Since the hypersurfaces orthogonal to u are spaces of constant curvature, the result follows 
immediately from proposition 11.2.31 □ 

We now show that the Weyl tensor of the Robertson- Walker metric is zero. 



Proposition 2.1.3 The Robertson-Walker metric di.il ) has zero Weyl tensor: 

Kb = 0, H ab = (2.8) 
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Proof: We make use of the contracted Ricci identities in the foi 




°(ab) = -2H(J ab + V( a U b ) + il( a U b ) - O £ a b)c ~ ^(a^b) - {Eab ~ 2 R (ab))> ( 2 - 9 ) 

= H ab - 2ii( a uj b ) - V( a cj b ) - curl(cr afe ). (2.10) 

We note that H ab = is an immediate consequence of ( 12.101 ), on applying proposition 12. l.lj If we apply 
proposition 12. l.ll to (|2T9l > and !2.1.2l to (11.421 ). we obtain 

E a b — \R(ab) = 0, E ab + \R(ab) = 0, (2.1 1) 

which implies E ab = 0. □ 

Finally, we derive the restrictions satisfied by the Ricci tensor of the Robertson-Walker metric d 1 .3b - 
Proposition 2.1.4 The Ricci tensor of the Robertson-Walker metric f |i. 3\) satisfies 

R ab u a h b c = 0, R {ab) =0. (2.12) 

Proof: We note that = follows immediately from (12. 111 ). In order to show R ab u a h b c = 0, we 

make use of the contracted Ricci identities in the forrro 

= V 6 a a6 - 2V a H - curIK) - 2e abc ii b oj c + R bc u b h c a , (2.13) 

Then on applying proposition 12. 1 . 1 1 to ( 12.131 ), we obtain R ab u a h b c = 0. □ 

We conclude this section by showing that the Cosmological Principle, as stated in section [TTT1 leads to 
the RW metric. 

Proposition 2.1.5 If the fundamental congruence of a cosmological model (M., g, u) satisfies 

o- ab = 0, io a = 0, ii a = 0, V a H = 0, (2.14) 



'These equations are an intermediate step in deriving iA.9l and iA. 13t . and follow from the Ricci identities, without use of 
the field equations. In fact, on applying the field equations dl.2t in the form dA.2b -l [A3t to l |2.9b and d2.10t . we obtain l lA.9l > and 

2 This equation is an intermediate step in deriving diA. lift , and follow from the Ricci identities, without use of the field 
equations. In fact, on applying the field equations dl.2t in the form JA.2l >-l lA"3t to d2.13t . we obtain dA.l It . 
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and the hypersurfaces orthogonal to u satisfy 

{3) S ab = 0, (2.15) 
then g and u can be written in the RW form di. JD awcf di.4D . 

Proof: By proposition ll.2.21 the fundamental 4-velocity u can be written as 

u a = -t a , (2.16) 

where t is a scalar field. We then introduce co-moving coordinates x a , i.e. coordinates x a that are constant 
along the world lines of the fundamental congruence. Then the metric and 4-velocity assume the form 

d 

ds 2 = -dt 2 + g a/3 (t, x x )dx a dx 13 , u = — . (2.17) 
Next, using the metric (12.171) . it can be quickly shown that 

VfeU a = \g a bfl- (2-18) 

This equation, in conjunction with dl.51 ), dl.221 ), (12.141 ) and ( 12.161 ), leads to 

g aP fl = 2H(t)g a p. (2.19) 
We introduce the length scale function £{t) via 

H(t) = j, (2.20) 

and thus can integrate (12.191 ) to obtain 

9af3(t,x x ) =£(t) 2 g a ^x x ). (2.21) 

Since we have that ^'S^ = 0, it follows from proposition 1 1 . 2 . 3 1 that g a p describes a 3-space of constant 
curvature. Thus the metric ( 12.171 ) assumes the Robertson- Walker form dl -3b - □ 



Theorem 2.1.1 (Derivation of the RW Metric) If a cosmological model (Ai, g, u) satisfies the Cosmo- 
logical Principle relative to u, then g and u can be written in the RW form di-il ) and di.4| ). 
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Proof: The Cosmological Principle implies that any geometrically defined vector orthogonal to u must 
be zero, since otherwise it would define a preferred direction orthogonal to u. In particular, we have 
w n = u a = V a H = 0. Likewise, any geometrically defined rank two symmetric tracefree tensor that 
is orthogonal to u (X a i,u b = 0) must be zero, since otherwise its eigenvectors would define a preferred 
direction orthogonal to u. In particular we have a ab = ^Sab = 0. The desired result thus follows from 
proposition 12.1.51 □ 

2.2 Characterizations of the FL Cosmologies 

In this section we present two characterizations of FL cosmologies in terms of the kinematic quantities 
and Weyl curvature. 

2.2.1 The Kinematic Characterization 

We now present the kinematic characterization of the Friedmann-Lemaitre models^ 

Theorem 2.2.1 A cosmological model (M,g,u) which satisfies the Einstein field equations with perfect 
fluid source is a Friedmann-Lemaitre cosmology if and only if the kinematic quantities of the fluid 4- 
velocity u satisfy 

a ab = 0, Lo a = 0, ii a = 0. (2.22) 

Proof. @ If we assume the model is a FL cosmology, the kinematic restrictions (12.221 ) follow immediately 
from proposition 12.1.11 

Conversely, assume that the restrictions (12.221 ) hold. Since the source is a perfect fluid with 4- velocity 
u, equation dl.281 ) implies -k^ = and q a = 0. Then the conditions (12.221 ) in conjunction with dA.91 ) imply 
= 0, which in turn implies ( 3 'S a b = via (1 1 -44b - The conditions (12.221 ) and q a = in conjunction 
with dA.l 11 ) imply V a H = 0. The desired result then follows from proposition 12.1.51 □ 



3 One can replace the assumption u a — in Theorem l2.2.1l bv the assumption that the perfect fluid has a linear equation of 
state p = (7 — l)fi. However, the proof requires a more detailed analysis and is beyond the scope of this thesis (see instead 
Collins and Wainwright (1983)). 

4 Krasinski (1997) gives the outline of a proof, referring to other sources for some parts (see page 11). 
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2.2.2 The Weyl Characterization 

We now present a characterization of the FL models using the Weyl tensor (see section [T.2.41 ). 

Proposition 2.2.1 If a cosmological model (M, g, u) satisfies the Einstein field equations with perfect 
fluid source having 4-velocity u and barotropic equation of state p = p(p) with fj, + p ^ 0, and has zero 
Weyl tensor, 

E ab = 0, H ab = 0, (2.23) 

then the kinematic quantities satisfy 

a ah = 0, u ab = 0, u a = 0. (2.24) 

Proof. On account of (12.231 ), equations dA.181 ), (IA.20I ) and dA.211 ) simplify to 

(p, + p)a ab = 0, V a ^ = 0, (n + p)u a = 0. (2.25) 
Since + p) 7^ 0, (12.251 ) implies cr a f, = and oj a = 0. We now make use of the chain rule to write 

VaP = ^V a /x = 0. (2.26) 

Then using the perfect fluid condition, i.e. 7r a i, = and q a = 0, it follows from (1A.15I ) and (12.261 ) that 
(fi + p)u a = 0, and hence u a = 0, which completes the proof. □ 

We can now present the main result of this section, which we shall refer to as the Weyl characterization 
of the Friedmann-Lemaitre models. 

Theorem 2.2.2 A cosmological model (Ai,g,u), which satisfies the Einstein field equations with perfect 
fluid source having 4-velocity u and barotropic equation of state p = p(fJ.) with /i+p 7^ 0, is a Friedmann- 
Lemaitre model if and only if the Weyl tensor is zero: 



E a b = 0, H a b = 0. 
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Proof. This result follows immediately from proposition 12. 1.31 and theorem \2.2. 11 in conjunction with 
proposition 12. 2. II □ 

We now present one final result which restricts the gradients of the physical and geometrical scalars 
in an FL cosmology. 

Proposition 2.2.2 A Friedmann-Lemaitre cosmology satisfies 

Va/J = V aP = V a H = V a (% = 0. (2.27) 

Proof: We note that V a H = V a ^R = follows immediately from propositions 12.1.11 and 12. 1 . 21 Fur- 
ther, it follows from (IA.4b . (1A.5I ) and proposition ^. 1 .4l that q a = and ir a b = 0. Hence, (IA. 15b and (IA.20I ) 
along with theorems 12.2. l| and |2.2.2| implies V a /i = V a p = 0. □ 



Chapter 3 



Dynamics of n-FLUiD FL 
Cosmologies 



The main goal of this chapter is a description of the dynamics of n-fiuid FL cosmologies. In section 13.11 
using the results from chapter [TJ we derive the basic evolution equations for n-fiuid FL cosmologies and 
describe their behaviour using dynamical systems techniques. In section 13.21 we describe the degrees of 
freedom of the ra-fluid FL cosmologies in two ways, firstly, using intrinsic parameters and secondly, using 
observational parameters. In sections [331 and l3~4l we apply the general analysis to the case of 2-fluid and 
3-fluid FL cosmologies. Finally, in section |331 we give a unified presentation of the known solutions of 
the EFEs for single fluid and 2-fluid FL cosmologies, and relate them to the state space analysis in section 
[33] 



3.1 Behaviour and Classification of n-Fluid FL Cosmologies 

In this section we derive the evolution equations for n-fluid FL cosmologies without introducing local 
coordinates and a line-element, by specializing the general evolution and constraint equations (appendix 
[At using the restrictions derived in chapter |2l We formulate the evolution equations as ODEs for a set 
of dimensionless variables in a compact state, which enables us to use dynamical systems methods to 
describe the evolution of n-fluid FL cosmologies qualitatively. 
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3.1.1 Basic Evolution Equations 

We know from theorems 12. 2. 11 and 12. 2. 21 and proposition 12. 1.41 that FL models have the following proper- 
ties: 

i) the kinematic quantities of the fundamental congruence u satisfy 

a ab = 0, uj a = 0, u a = 0, (3.1) 

ii) the Weyl curvature is zero: 

E ab = 0, H ab = 0, (3.2) 

and 

iii) the stress-energy tensor satisfies 

7T a fc = 0, q a = 0, (3.3) 

i.e. it has the form of a perfect fluid. 

The remaining non-zero physical and geometric variables that appear in the general evolution and 
constraint equations in Appendix |A] (1A.8I) - (1A.21I ) are the Hubble scalar H , the energy density fj, and the 
pressure p of the matter content. Propositions 12.1.11 and 12.1.41 show that these scalars have zero spatial 
gradient 

h*H b = 0, h*n >b = 0, h* Pib = 0, (3.4) 

and hence are purely functions of clock time t along the fundamental congruence. With the above re- 
strictions, the general evolution and constraint equations reduce to two evolution equations for p, and H, 
namely dA.8t and (IA. 14b . which read 

p = -3H(n+p), (3.5) 
H = -H 2 -^ + 3p), (3.6) 

where the overdot denotes differentiation with respect to t. The spatial geometry of the model is deter- 
mined by the Gauss equation in the form (11.431 ) and (11.44b . It follows that the tracefree Ricci tensor is 
zero, 

(3) Sab = 0, (3.7) 
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and that the spatial Ricci scalar is determined by H and p, according to 

(% = -6H 2 + 2fi. (3.8) 

In order to obtain the standard evolution equations, we introduce the length scale function £(t), which 
is defined by H, up to a constant multiple, according to 

£ 

H = -. (3.9) 



This definition is consistent with the length scale defined by the RW metric (11.3b . 
To proceed we need the following result: 

Proposition 3.1.1 In a Friedmann-Lemaitre model the quantity 

WR(t)£(t) 2 (3.10) 

is constant in time. 

Proof. A straightforward calculation using (ESJ-OU yields (^R£ 2 )' = 0. □ 

Since £ is only defined by ( 13.9b up to a multiplicative constant, we may choose it to satisfyQ 

i (z)R£ 2 =K = ±1, or 0, (3.11) 

where K is the spatial curvature parameter. The value for K indicates whether is positive, negative 
or zero. Models where K = —1, and +1 are then referred to, respectively, as open FL models, flat FL 
models and closed FL models. 

It is important to note that if K = ±1, then £{t) is completely determined via eq (13.1 lb . by the spatial 
geometry, represented by ^R. On the other hand, if K = 0, equation (13.11b does not restrict £{t), which 
means that if K = then £{t) can be rescaled with a multiplicative constant. 

We now present the standard form of the evolution equations for FL models (see, for example, Wain- 
wright and Ellis 1997, p52). Using ( |3.9b and ( 13.111 ), equations ( 13.5b , ( 13.6b and ( 13.8b are rewritten as the 
conservation equation, 

fi = -3%+p), (3.12) 



'The factor of g is chosen for compatibility with the RW metric. 



3.1. Behaviour and Classification of n-Fluid FL Cosmologies 



23 



the Raychaudhuri equation, 

\ = ~W + ZP), (3-13) 

and the Friedmann equation, 

3£ 2 = [if - 3K. (3.14) 

Note that (13.131 ) is a consequence of (13.121 ) and (13.141) provided that I ^ 0, as follows by differentiating 
(EH). 

3.1.2 rz-Fluid FL Cosmologies 

In this chapter, we consider FL cosmological models whose matter content consists of n fluids with the 
same four-velocity u. The total energy density and pressure are equal to a sum of the energy density and 
pressure of the individual components: 

n n 

fi = ^2fH, P = ^2pi- (3-15) 

i=l i=i 

The fluids are assumed to be non-interacting, which means that each fluid component individually satisfies 
the conservation equation (13- 1 2b . i.e. 

e, 

= -3-(/ij (3.16) 

We further assume that each fluid obeys a linear equation of state with equation of state parameters 73 
(i = 1, . . . , n), i.e. 

Pi = (7j - I) m, (3.17) 

with 

0<7i<2. (3.18) 

The most important physical values are 7^ = 1 for pressure-free matter and 7i = | for radiation. The 
value 7j = can be used to represent a cosmological constant and 7j = 2 is sometimes considered to 
correspond to a so-called stiff fluid. For definiteness we order the equation of state parameters according 
to 

71 > 72 > ••• > 7„, (3.19) 

and assume that 71 > 1, i.e. at least one matter component is matter with a non-negative pressure. The 
conservation equations ( 13.161 ) then assume the form 

Ai = -3~7i/Xj. (3.20) 
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As well, the Friedmann equation ( 13.141) and Raychaudhuri equation (13.131 ) assume the form 

n 

P = J2l^ 2 -K, (3.21) 

i=i 

and 

I n 

- = -i^(3 7! -%. (3.22) 
i=i 

In summary, we will use the term n-fluid FL cosmology to refer to an FL cosmological model whose 
matter content consists of n non-interacting fluids with the same 4- velocity u and equations of state given 
by (13.17l) - d3~T9l ). Following Jantzen and Uggla (1992), we will use abbreviations such as i?DC-universe 
and i?(7A-universe to identify n-fluid FL cosmologies, where R denotes radiation, D denotes dust, C 
denotes curvature and A denotes a cosmological constant. 

We have shown that the state of an n-fluid FL cosmology as a function of time t is determined by 
the functions (£, ... , fi n ) which satisfy the system of ODEs (13.201) and (13-2 1 b - The state of an n-fluid 
cosmology can equivalently be described by the variables (H, n\, . . . , /i n ), which satisfy the evolution 
equations 

fii = -3H"fim, (3.23) 

n 

H = -H 2 - \ ^(3 7i - 2)m. (3.24) 

2=1 

The first of these is a restatement of (13.201) using the definition ( 13.91 ) of H. The second follows from 
(13.61 ) on using (13.151 ) and (13.171) . Further, upon applying (13.9I ). the Friedmann equation (13.211) gives an 
additional constraint that must be satisfied, namely 

n K 

i=l * 

We note that this equation explicitly determines £ in terms of H and fi when K ^ 0. 

At this stage, we present four simple and well-known solutions of equations (I3.20I )- (I3.22I ) that play an 
important role when we formulate the evolution equations (|3.23b - d3T24T )for the n-fluid FL cosmologies as 
a dynamical system. 

1. Flat FL universe (K = 0, < 7 < 2). For a single fluid, the length scale, energy density and 
pressure are 

^ = 4 (£) 37 ' ^ = WP' P=(7-l>- (3-26) 
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In the case of dust (7 = 1), this solution is the Einstein-de Sitter universe (Einstein and de Sitter 
(1932)). 

2. de Sitter universe (K = 0). The de Sitter universe can be thought of as the special case of the flat 
FL universe when the single fluid is a cosmological constant {i.e. when 7 = 0). The length scale, 
energy density and pressure are 

£ = £ exp (J^tj > A* = A, p = -A. (3.27) 

3. Milne universe (K = —1). The length scale, energy density and pressure are 

£ = t, (j, = 0, p = 0. (3.28) 
The Milne universe can be thought of as the open FL universe with zero matter content. 

4. Einstein static universe (K = +1). This solution arises when we require that the length scale is 
constant. The length scale, total energy density and total pressure are 

£ = £ , v = 3£v 2 , P=-£o\ (3-29) 

where £q > is constant]^] The matter-energy content can be interpreted as two, or more generally, 
n non-interacting fluids. For example, for a 2-fluid model, the energy densities are given by % 

Ml = K 1 M2 = ^0 > ( 3 " 3 °) 

V71-72/ V71-72/ 

where < 72 < | < 71 < 2. Einstein's choice was 71 = 1 (pressure-free matter) and 72 = (a 
cosmological constant) so that fii = 2£q 2 , p\ = and \ii = £$ 2 = A, P2 = —A, where A is the 
cosmological constant. 



2 The values of /1 and p follow l |3.13t and d3.14| l. 

3 These expressions are obtained by solving the linear system 



Hi + M2 = 3£ 2 , (3 7 i - 2) M i + (372 - 2)^2 = 0, 



which follows from d3.2U and d3.22t with I — Iq and n — 2. 
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3.1.3 Hubble-Normalized Scalars 

In this section we introduce so-called Hubble-normalized scalars, which form a set of dimensionless 
variables that, together with the Hubble scalar H, describe the state of n-fluid FL cosmologies. These 
dimensionless scalars also play an important role in relating the models to observations (see section |3T2.2| ). 

Firstly, the density parameters Qi, which describe the dynamical effect of the matter densities, are 
defined by 

n< = 5§a- O-s i) 

The total density parameter Q, which describes the dynamical effect of the total matter content, is defined 
by 

n = 3^. (3.32) 

It follows from (l3~l3T ) and (I33TT) that 

n 

n = J2^i- (3.33) 

i=l 

Lastly, the dynamical effect of the spatial curvature is described by Q, k , defined by 
It follows from dBTTTb that 

n k = -j^. (3.35) 
In terms of the parameters f2 and O^, the Friedmann equation (13.211) assumes the simple form 

n + n k = i, (3.36) 

which may be further expanded in terms of fi, as 

n 

^^ + fi fe = l. (3.37) 
i=i 

We note that (13.361 ) implies that Q, determines the spatial curvature as follows: 

n < i open FL, (3.38) 

tt = 1 flat FL, (3.39) 

n > i dosed FL. (3.40) 
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There is an additional dimensionless variable that plays a fundamental role, namely the deceleration 
parameter q, defined by 

it 

q = . (3.41) 

(£) 2 

Clearly, if q < then the expansion is accelerating (£ > 0) and if q > then the expansion is decelerating 
(l < 0). In an n-fluid FL model, q can be expressed in terms of the density parameters by using the 
Raychaudhuri equation (13.221 ) with (13.91) and (13.31I ). We write 

n 
i=l 

The quantities fij, O, and q are dimensionless, whereas H has dimensions of (length)^ 1 . Together, 
these variables are fundamental quantities in the FL models and more general models (see, for example, 
Wainwright and Ellis (1997), section 5.2 and chapters 12 and 13). 

The evolution equations for Oj and can be obtained by differentiating (13.311) and (13.351 ) and apply- 
ing (13- 14b . (13.231) . (13.241 ) and (13.421) . On performing this procedure, we obtain 

^ = (2q - (37< - 2))SkH, and ^ = 2q£l k H. (3.43) 
3.1.4 A Compact State Space for n-fluid FL Cosmologies 

The one drawback of using H-normalized variables is that they cannot give a unified description of models 
that have both an expanding (H > 0) and a contracting (H < 0) epoch. In particular, when H = the 
density parameters (13.311) and (13.341 ) are undefined {i.e. Qi — ► +oo as H — > 0). In order to circumvent 
this difficulty, we introduce a set of variables which are bounded throughout the evolution of the model. 
For an n-fluid FL cosmology, with total energy density 

n 

li = J2m>0, (3.44) 

1=1 

we define the dimensionless matter variables, which describe the relative significance of the different 
matter components: 

Xi = ~, IM > 0. (3.45) 
A* 

Then by (13.441) we have 

n 

£>i = l, 0<x*<l, (3.46) 
i=l 
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i.e. the Xi are bounded. 

Following Wainwright (1996, pl24-5), we now introduce a curvature parameter (l so as to compactify 
the state space: 

V3H S 



It follows that SI satisfies 



n = arctan . (3.47) 



\ < & < \ , (3-48) 



and that 

sign(ft) = sign(F). (3.49) 

This variable thus describes the state of expansion of a model. In addition, it describes the spatial curva- 
ture, i.e. the deviation from flatness, as a consequence of (13.321 ) and (l3-38b - (13-40b - In particular, it follows 
from (13391 and (13321) that U = \ gives the flat FL models. 

We now derive evolution equations for Cl and the Xi- It i s convenient to introduce a quantity q that 
plays the role of the deceleration parameter q, but unlike q, is bounded when H = 0. Motivated by ( 13.421 ), 
we define 

n 
i=l 

Together, ( 13.191 ) and (13.501 ) imply that q is restricted by 

^(3 7 n-2) <g< i(3 7 i-2). (3.51) 
It follows from (13.151 ). (13.231) . (13.451 ). (13.461 ) and (13.501 ) that the total energy density satisfies 

^- = -2(q + l)H». (3.52) 

at 

The evolution equation for Xi then follows from (13.231 ). (13.451) . (13.501 ) and (13.521 ) and is given by 

^ = [2q-(3rfi-2)]H X i. (3.53) 
at 



Upon noting that ( 13.471 ) implies 
we can rewrite (13.531) as 



// J^tanO, (3.54) 



</V ' ' -[2g-(3 7i -2)](sinJ%;. (3.55) 



dt V 3 cos 
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The form of ( 13.55D suggests defining a new time variable f by 

= /-cosO. (3.56) 

CIT V 

It is essential to note that this change of variable is well-defined, since cos 0. > on account of (I3.48I ). In 
terms of f , (13.551 ) reads 

^ = [2q - (3 7i - 2)](smn)x«. (3.57) 

In order to obtain the evolution equation for Q we need the evolution equations for fi and H. First, it 
follows from d3~32l and (13361 ) that 

— = -2J3flH(q + 1) cos ft. (3.58) 

(XT 



Secondly, equations (l3T24l >. d3~45l >. (l3~30l) . (13341 and (13361) lead to 

(tan 2 n + q) cos f2. (3.59) 



dH l~Vu„. 2 



V 3 

We differentiate (13.541 ) with respect to f and use (13.581) . After simplifying using (13.54I ). we obtain 

— — = — q cos 2il cos 17. (3.60) 

Equations ( 13.571 ) and (13.601 ) are the desired evolution equations. The state space is the bounded n- 
dimensional subset of M n+1 described by (O, xi, . . . , Xn) an d subject to the restrictions (13.461 ) and (13.481 ). 
The state space is not compact due to the strict inequality in (13.481 ) and so we seek to compactify this 
space, since then we can apply theorems from the theory of dynamical systems. 

In order to compactify the space, we simply observe that the right hand sides of the ODEs (13.571 ) and 
( 13.601 ) are well-defined, in fact analytic, at the points = ±|, which form part of the boundary of the 
state space. We can thus extend the state space to include the points with £l = ±|, thereby compactifying 
it. Since 

lim n = 0, (3.61) 

|n|-f- 

as follows from (13.541) . we interpret the sets Cl = ±2 as representing FL cosmologies with \i = 0. It 
follows from ( 13.141 ) that the set Q = ±| represents the Milne universe, given by (13.28b - 
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State Space Representation of n-Fluid FL Cosmologies 




State Space: 




(n,xi,...,xn)em. n+1 






n 

/ j A.J 
8=1 


1 < Yi < 1 -S < f3 < S 


(3.62) 


Evolution Equations: 










dn 
~dJ 

dXi 
df 


= — Q COS iCl COS Cl, 

= [2g-(3 7i -2)](smO)xi, 


(3.63) 
(3.64) 


with 




i=l 


(3.65) 



Comments: 

i) The ODEs (13.631 ) and (13.641 ) describe the essential dynamics of n-fluid FL cosmologies in terms of 
the dimensionless variables (f2, xi, • • • , Xn) and the dimensionless time f . Since the state space, as 
defined by (13.621) is compact, the solutions of the ODEs are defined for all f £ R, and as a result 
the ODEs define a dynamical system on the state space (see, for example, Wainwright and Ellis 
(1997), p87, Corollary 4.1). The evolution of an n-fluid FL cosmology is thus described by an orbit 
of the dynamical system (a solution curve of the ODEs). As a consequence, we are guaranteed the 
existence of a past attractor A~ of the ODEs {i.e. an invariant set to which all orbits, except possibly 
a set of measure zero, are past asymptotic to as f — ► — oo) and a future attractor A + {i.e. as a past 
attractor, except as f — * +oo). The past attractor describes the asymptotic regime at early times and 
the future attractor describes the asymptotic regime at late times. 

ii) This description of the evolution of n-fluid FL cosmologies as a dynamical system on a compact 
state space is new. It generalizes the formulation given by Wainwright (1996) for the case n = 2. 

hi) The physical state of an n-fluid FL cosmology is determined by the variables (H,fj,i, . . . , /x n ), 
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which are related to the dimensionless variables via the total density fx, according to 




(3.66) 



(see (13.541 ) and (13.45I )). On using (13.661 ) to eliminate H, (13.581 ) assumes the form 



djj, 



2(q + l)(sinO)|U. 



(3.67) 



df 



Then (13.671 ) determines n as an integral in terms of Cl{f) and Xi(^)> U P to an arbitrary constant. 
This freedom implies that each orbit determines a 1 -parameter family of cosmological models, 
whose physical variables are related by a multiplicative constant. Since an n-dimensional dynamical 
system has an (n — l)-parameter family of orbits, it follows that the class ofn-fluid FL cosmologies 
is labelled by n parameters. We will specify these parameters explicitly in section 13.21 

iv) Clock time t is related to f via the ODE (13.561 ) and can thus be expressed as an integral involving 
0(f) and \jl{t ), the latter determined from ( I3.67I ). Clock time will take on values in an interval 



where ij is finite or — oo and tf is finite or +oo, depending on the model. 

v) The dimensionless variables Cl, x% an d <j are closely related to the Hubble-normalized scalars Qi 
and q when H ^ 0. From (13.321 ) and (13.471 ) we have 



(3.68) 



tan 2 il 



1 

tt 1 



| < O < f , O t^O. 



(3.69) 



Similarly, from (133TT ), d3"32l and d3l3T ) we have 




(3.70) 
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3.1.5 Invariant Sets and Equilibrium Points in the Compact State Space 

Our goal in the rest of this section is to describe the qualitative behaviour of the orbits in the compact state 
space, thereby understanding the dynamics of n-fluid FL cosmologies. Since the state space is compact, 
it follows that every orbit has a past and future attractor (see, for example, Wainwright and Ellis (1997), 
p99-100). 

In general the attractors of a dynamical system are complicated invariant sets. But in the present 
situation we will show that they are determined by the equilibrium points of the dynamical system (section 
13- 1 -7b . which we now present. 

The form of the ODEs ( 13.631 ) and ( 13.641 ) determines three physically important invariant sets, given as 
follows: 

1) Flat FL set. This set is the n — 1 dimensional subset given by Cl = ±f (which corresponds to 
O = 1, fifc = on account of ( 13.691 ) and (13.36b ). The orbits in this set describe flat n-fluid FL 
cosmologies, expanding if Cl = | and contracting if Cl = — j. 

2) Milne set. This set is the n — 1 dimensional subset given by Q = ±| (which corresponds to f2 = 
and Qk = 1 on account of ( 13.69b and (13.36b ). The orbits in this set describe the Milne vacuum 
solution. 

3) m-fluid sets (1 < m < n). These sets are determined by setting n-mof the \i variables to be 
zero, with the m remaining %i variables being non-zero. The orbits in these sets describe m-fluid 
FL cosmologies. 

The equilibrium points are determined by setting the right hand side of equations (13.631) and (13.641 ) to 
zero. If O 7^ 0, then ( 13.641 ) and ( 13.65b implies \i = 1» f° r one value of i, \j = 0> f° r 3 7^ an d 
q = I(3 7i - 2) / 0. Then (131)31 implies cos & = {i.e. Q = ±§ ) or cos 2^1 = (i.e. Q = ±f ). If 
Q, = then (13.63b implies q = 0, with q given by (13.65b . In this way, we obtain the following sets of 
equilibrium points: 

1) Flat FL points ( F^) are given by 

f} = ±vr/4, Xi = l, Xj = 0iij^i. (3.72) 

These equilibrium points lie in the intersection of the flat FL set and the various 1-fluid sets. 
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2) 



Milne points ( Mf ) are given by 




Xi = l,Xj = if j / i. 



(3.73) 



These equilibrium points lie in the intersection of the Milne set and the various 1 -fluid sets. 



3) 



Einstein static points (E) are given by 



n 



n = o 



^(3 7i - 2) Xl = 



(3.74) 



i=l 



and so constitute an (n — 2)-dimensional set of equilibrium points. These equilibrium points corre- 
spond to the n-fluid interpretation of the Einstein static model (13.291 ) and only exist when 7„ < 2/3. 

3.1.6 Classification of n-Fluid FL Cosmologies 

In this section, we will prove that the n-fluid FL cosmologies that have an epoch of expansion (i.e. H > 
in some time interval) form three qualitatively different generic^ subclasses, defined as follows: 

1) Ever-expanding models satisfy H > for all t. 

2) Recollapsing models satisfy H(t*) = 0, H > for t < i* and H < for t > t* for some t*. 

3) Bouncing models satisfy H(t*) = 0, H < for t < i* and H > for t > t* for some i*. 

The variables Cl and q are particularly important in this analysis. We begin with the following propo- 
sition on the behaviour of q. 

Proposition 3.1.2 In any n-fluid FL cosmology with n > 1, q is strictly decreasing (increasing) in an 
epoch of expansion (contraction). Further, ifn = l then q is constant. 



invariant sets of dimension n. 
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Proof: By differentiating (13-501 ) with respect to f and using ( 13.571 ), we obtain^] 



df 



4<f - ^(3 7i - 2) 2 Xl 



i=i 



sinO. 



By the Cauchy-Schwarz inequality and (13.461 ) we have 



(3.75) 



4<f 



Xi 



i=l 



< 



E( 3 ^ - 2 )' 



i=l 



with equality holding if and only if n = 1. @ 
On recalling (I3.49I ). it follows that 

sign 



dq 



Xi 



i=i 



Xi 



^(3 7l - 2) 2 Xi , (3.76) 



i=l 



-sign(F), 



(3.77) 



which establishes the result. □ 



This proposition leads to a second result, presented here: 

Proposition 3.1.3 For any non-static n-fluid FL cosmology, fi(f#) = is satisfied for at most one f*, and 
in this situation ^p(f*) 7^ 0. 

Proof: We first show, by contradiction, that f2(f*) = is satisfied for at most one f*. Suppose 
f2(fi) = = f2(f2) and O 7^ for f G (fi,T2). Without loss of generality, we assume that Cl > 
for f G (fi,?2). Then by proposition 13. 1.21 it follows that ^(fa) < g(fx) since g is strictly decreasing 
whenever $7 > 0. Further, by (13.601 ). we have that £1 satisfies 

n'(ri) = -?(n), and jV(r 2 ) = -g(f 2 ), (3.78) 



5 This equation is equivalent to (14.15) in Wainwright and Ellis (1997), since q = |(3u> + 1). The proposition is equivalent 
to <? s — w > 0, if n > 1. 

6 Equality holds in d3.76b if and only if there exists c such that 

(37; - 2)^/^7 = c^/xl, i=l,2,...,n, 



which is possible if and only if n = 1. 
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where the prime (/) denotes differentiation with respect to f. But 0(f) > for f 6 (ti,^) implies 
that fi'(fi) > (i.e. g(fi) < 0) and fV(f 2 ) < (i.e. g(f 2 ) > 0) and so q(f 2 ) > q{f\), which is a 
contradiction. 

Finally, it follows from (13.631 ) that if VL = 0, *p = is satisfied only when g = 0, which is disallowed 
for non-static models. □ 

The main result of this section now follows immediately. 

Theorem 3.1.1 There are exactly three generic classes of n-fluid FL cosmologies (n > 1) that have an 
expanding epoch, namely ever- expanding models, recollapsing models and bouncing models. Bouncing 
models occur if and only if ' 7n < |. 

Proof: Since H = if and only if (l = on account of (I3.47I ). proposition 13 . 1 . 3 1 states that the Hubble 
scalar can be zero at most once, and that H changes sign if it becomes zero. Thus models with more 
than one expanding epoch are excluded, which implies that the only possible classes are those given. If 
7„ > |, it follows from ( 13.191 ) and (13.501 ) that q > |(37„ - 2) > 0, and hence from (13.631) that bouncing 
models are disallowed. □ 

To the best of the author's knowledge, the classification theorem is new for n > 2, although it is a well 
known result for n = 2, when the matter content is a perfect fluid (71 > |) and a cosmological constant 
( 72 = 0) (see Rindler (1977), p234-8). 

3.1.7 Asymptotic Behaviour of Solutions 

In this section we determine the asymptotic behaviour of the three generic classes of n-fluid FL cosmolo- 
gies identified in the previous section. We will first describe the behaviour of the length scale £(t) as a 
function of clock time t in each of the three generic classes of FL cosmologies, and then use these results 
to determine the past attractor and future attractor. 

Important in this analysis is the set of orbits asymptotic to an Einstein static solution. In the following 
proposition, we demonstrate that orbits of this type are non-generic. 

Proposition 3.1.4 The set of orbits in the compact state space that are asymptotic to an Einstein static 
equilibrium point has dimension n — 1. 
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Proof: We can write the evolution equations ( 13.571) and ( 13.601 ) as a dynamical system of the form 

dXi 



df 



F i (X ,X 1 ,...,X n ), i = 0,...,n, (3.79) 



on R n+1 , where Xq = ti and Xj = Xj (j = 1) • • • ; n )- The physical state space is the hyperplane given 
by (I3.62I ). and hence is n dimensional. Consider the Jacobian matrix Jij = Upon evaluating Jij on 
the Einstein static set, given by ( 13.741 ), we obtain 

Joo = 0, J 0i = -|(3 7i - 2), J i0 = -(37i - 2) X i,E, Jij = 0, (3.80) 

where f, j = 1, . . . ,n and denotes the value of x% at the Einstein equilibrium point. The Jacobian 
matrix, when restricted to the physical state space, has two non-zero eigenvalues 



A - -± \ YP>li ~ 2) 2 Xi > 0, (3.81) 

and n — 2 zero eigenvalues]^] It follow^ that each Einstein static equilibrium point has a 1 -dimensional 
stable manifold and a 1 -dimensional unstable manifold. Thus the set of orbits future asymptotic to the 
Einstein static set has dimension exactly one larger than the dimension of the Einstein static manifold, i.e. 
dimension n — 1. □ 

In the next proposition, we establish the behaviour of the length scale fas a function of time f , for 
ever-expanding models. 

Proposition 3.1.5 Iffl> Ofor all f G R and & -/-> as f —> ±oo, then 

lim 1(f) = — oo, and lim 1(f) = +oo. (3.82) 

T — ► — OO T — > + 00 

Proof: The length scale I defines a time variable r according to 

I = £ e T . (3.83) 



7 The zero eigenvalues reflect the fact that the set of Einstein static equilibrium points has dimension n — 2. 
8 Using the Stable/Unstable Manifold Theorem. See, for example, Perko (1996, pl07). 
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It follows from COS) and (I3T831) that 

dt _ 1 
dr ~ H' 



(3.84) 



Then, using (13.54I ). (13.561 ) and (13.841) . and the chain rule, we obtain 



dr 
df 



(3.85) 



Since £1 > for all f G R and Cl -f* as f — > ±oo, (13.851 ) then implies 



lim r(f) 



do, and lim r(f) = +oo. 

f ^+00 



(3.86) 



This completes the proof, on applying (13- 83 b - □ 



Proposition 3.1.6 For a recollapsing model, q(f) > for all f, and for a bouncing model, q(f) < Ofor 



Proof: For a recollapsing model, it follows from the definition that there exists f* so that = 

and that &(f) changes sign from positive to negative at f*, which implies ^P(t*) < (strict inequality 
follows from Proposition 13. 1.31) . The evolution equation (13.631) now implies g(f*) > 0. By Proposition 
13.1.21 q has a global minimum at f*, i.e. q(f) > q(f*) for all f, which gives the desired result. The proof 
for bouncing models is similar. □ 

We now determine the asymptotic behaviour of clock time t (i.e. as f —> ±oo). Let 



where the dependence of t on f is determined via equation (I3.56I ). Then U and tf are finite or infinite 
depending on the type of universe and during the evolution ij < t < tf. Clock time is related to the length 
scale via the Friedmann equation. In particular, for any n-fluid FL cosmology, we can integrate (13.201 ) to 



all f. 



t{ = lim i(f ) 



tf = lim t{f) 



(3.87) 



T — > — OO 




iiit 2 = 3 ( V) 



-3 7i +2 



(3.88) 



9 See p46. 
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where Aj is an integration constant whose inverse can be interpreted as the characteristic length scale of 
fluid i. On substituting (13.881 ) into the Friedmann equation (13.211) . we obtain 

|) 2 =E(vr 37 * +2 -if. ( 3 - 89 ) 

Theorem 3.1.2 The following statements are true: 

i) For an ever-expanding n-fluid FL model that is not asymptotic to the Einstein static universe, U is 
finite, tf = oo and the length scale satisfies £{U) = Ofor some t{ and £ — ► +oo as t — > +oo. 

ii) For a recollapsing n-fluid FL model, t{ andtj are finite, l{ti) = and £(tf) = 0. 

Hi) For a bouncing n-fluid FL model, U and t f are infinite and £(t) satisfies £(t) — > +oo as t — > ±oo. 

Proof: Case i) In this case, we will make use of the following elementary results from calculus: 
(I) If f(a) = 0(a r ) as a — > 0, and r > — 1, then lim / f(a)da < oo, i.e. the integral converges. 

ra 

(II) If f(a) E 0(a r ), as a — > oo, and r > — 1 then lim / f(a)da = oo, i.e. the integral diverges. 

a-*ooj 1 

Recall that ever-expanding models satisfy H{f) > V f. Thus by Proposition 13.1.51 solutions which 
are not asymptotic to the Einstein static solution satisfy £ — >0asf— >— oo and £ — > oo as f — > oo. Let 

f(t) = 1 (3.90) 

T?=i(vr 37i+2 -^ 



For an ever-expanding model, this function is defined and continuous (for < a < oo) since H > 
implies that the denominator is strictly positive (see ( I3.89I )). Hence, we can integrate ( 13.891 ) to obtain 

t{£) = t(l) - / f(t)d£, for < £ < 1. (3.91) 



Further, the integrand satisfies 

/ 371 -2 \ 

f{l) = O [£—^- J , as £ -» 0. (3.92) 

Since 71 > | it follows from result (I) that there exists some finite time ij, defined by t{ = lim^^o t(£), 
such that £{ti) = 0. 
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Similarly, we can integrate ( 13.891 ) to obtain 



t{£) = t(l) + / f(£)d£, for £ 



> 1. 



(3.93) 



J i 



Further, the integrand satisfies 



f(a) = O 



as£ 



oo. 



(3.94) 



Since 7 n > it follows from result (II) t — * oo as £ — > oo. 

For the following two cases we make use of the fact that q and I are related by 



£ 
1 



(3.95) 



as follows from d522l . (I3T45T) and (f3750T> . 

Case ii) For a recollapsing model, Proposition 13.1.61 implies q(f) > for all f G R. Then ( 13.951 ) implies 
that ^(t) is concave down on its domain and so we can conclude there exists ti and tf such that £(ti) = 
and £(t f ) = 0. 

Case iii) For a bouncing model, Proposition 13.1.61 implies q(f) < for all f G R. Then (13.951 ) implies 
that is concave up on its domain and hence £ — > oo as f — > ±oo. Since sign(£) = sign(i^), it follows 
that £ is minimal at f*. If we choose £ m in = £(t*), we have that £(f) > £ m in > for all f G R. 

We must now show that ^ does not go to infinity in finite time. Choose £ so that H > for £ > £ . 
Then f(£) in (13.901 ) is defined and continuous for £ > 1. We can integrate (13.891 ) from £ = 1 to obtain 
( 13.931) . Since the integrand satisfies 



and 7 n > 0, it follows from result (II) that (13.931 ) diverges, and so these models satisfy £ — > oo as t — > oo. 
By symmetry, we have that £ — ► oo as i — > — oo. □ 

The past and future attractors of the dynamical system on the compact state space can now be deter- 
mined using theorem f3. 1.21 




as a 



oo 



(3.96) 
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Theorem 3.1.3 The past and future attractors in the state space of n-fluid FL cosmologies with an ex- 
panding epoch are 

A' = F+ U and A + = Ff U F+, for 7ri < §, (3.97) 

and 

A- = F+, and A + = F{ U M+ for 7 „ > §,. (3.98) 

Proof: We begin by obtaining algebraic expressions for \i an d ^- From (13.881 ) we have 

ml 2 = 3(V)~ 37l+2 . (3.99) 

Then d3~45l >. d3T88T> and (ET991 lead to 

(V)-^ +2 mnm 
" £™ = i(V)~ 37j+2 

Using ( I3.47I) , we write the Friedmann equation ( I3.25b in the form 



l-tan 2 ft = — T . (3.101) 



Then using (13.451 ) and ( 13.991 ) we have 



l-tan 2 f)= ^ n JL-—^ . (3.102) 



To determine the limits of ( 13.1001 ) and ( 13.1021 ) as I — > 0, we write 

Xi = — _L. +2 (3-103) 

and, using g33) , (13^991 and (l3~T00l) . 

1 - tan 2 = ^xi(Ai^) 371 " 2 . (3.104) 

It follows that 



since 71 — 7^ > and 71 > |. 
Similarly, by writing 



lim xi = 1, lim Q = ±f , (3.105) 

£-►0+ ^->0+ 



Xn = — -3-V-+2 7' (3.106) 
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and 

1 - tan 2 Cl = Kxn{\n£f ln ~ 2 , (3.107) 

we obtain 

f 0, if 7„ < | 

lim X n = l, lim fl={ I (3.108) 

t-++°o [if, if 7n > | and K = -1. 

It now follows from Theorem 13.1.21 and equations (13.721 ) and (13.731 ) that orbits describing models in 
the three generic classes approach the equilibrium points given in table l3.ll as f — ► — oo and as f — > +cxd. 
□ 



Comment: Since the orbits describing models in the three generic classes approach the equilibrium 
points in table 13. II as f — ► ±oo, it follows that the density parameters Qi and approach their values at 
the equilibrium points, which we give in table 13.21 
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Model Type f — » — oo f — > +00 
Ever-expanding models 

In < § F+ F+ 

ln >l,K = -l F+ M+ 

Recollapsing models F±~ F-f 

Bouncing models (-y n < | only) F~ F+ 

Table 3.1: The past and future attractors of generic n-fluid FL cosmologies. 



Equilibrium Point Behaviour 



Ff -> for i / 1, fix -► 1, O fc -► 0, 

F± fii -»• for i ^ Ti, ^ n 1, ft fc -> 0, 

ikf ± ft* -> Ofori = 1, ... ,n, -> 1. 



Table 3.2: Asymptotic behaviour of Hubble-normalized quantities in generic n-fluid FL cosmologies. 
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3.2 Parameterization of FL Models 

In this section we introduce the n essential parameters that label the family of n-fluid FL cosmologies, 
namely a conformal parameter A and n — 1 dimensionless mass parameters mi, . . . , m n _i, that we shall 
refer to as intrinsic parameters. In section |3~.2.2I we introduce an equivalent set of parameters, namely the 
observational parameters O^o and Hq and show how they determine the intrinsic parameters. Finally, we 
show that the evolution equations (13.571 ) and (13.601 ) admit n — 1 independent conserved quantities, which 
characterize the orbits in the n-dimensional state space. In particular, we will show that the values of the 
conserved quantities are determined in terms of the intrinsic parameters. 

3.2.1 The Intrinsic Parameters 

We now integrate the conservation equation and demonstrate that the constants of integration of this equa- 
tion can be used to parameterize n-fluid FL cosmologies. As we have seen previously in (I3.88I ). on 
integrating (I3.20I ). one obtains 



The constants of integration Aj have dimension (length) 1 and so each can be thought of as defining a 
characteristic length scale for each fluid via A^ 1 . When ( 13.1091 ) is substituted into (13.211 ). there results a 
first order ODE for £(t): 



In order to define a dimensionless length scale and time variable for n-fluid FL cosmologies, one needs 
to choose a preferred Aj that will play the role of an overall conformal factor. Typically we will choose A 
to be A n , which is a natural choice if 7„ = 0. In this case, it follows from ( 13.1091 ) that 



fiifi = 3 (Xd) 



-3 7i +2 



(3.109) 




(3.110) 




(3.111) 



The dimensionless length L and time T are then defined as 



L = XI 



and T = Xt. 



(3.112) 



Equation ( 13.1091 ) suggests defining a set of dimensionless parameters mi, . . . , m n _i vi; 




(3.113) 



10 



'For i = n, this equation states that m 
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so that 

fi i = 3X 2 m i L- 3 ' y \ (3.114) 
In dimensionless form, the Friedmann equation ( 13. 1101 ) then reads 

' i=l 

For any choice of parameters rrii > 0, i = 1, . . . , n — 1, along with the conditions 

L(0) = 0, and ^ > 0, (3.116) 
the DE ( 13.1151 ) has unique solution (as shown in section |C2l ), given by 

L = L(T, mi, ... , m n _i). (3.117) 
It then follows from ( 13.1121 ) that £(t) can be written as 

l(t) = X' l L{Xt, mi, ... , m„_i). (3.118) 
The line element then takes the form 

ds 2 = \~ 2 (- dT 2 + L 2 dY?), (3.119) 

where L is given by (13.1171 ) and dY? is a metric of constant curvature. From (13.91 ), (13.1121 ) and (13.1151) 
we obtain 

H 2 = X 2 mL~ 3 ^ - KL- 2 ^j . (3.120) 

Further, using (13.311 ), ( 13.1141 ) and ( 13.1201 ) the density parameters Qi are written as 

mi L-^ +2 

= r _3 7 . + 2 T? ' (3 " 121) 

2^j=i m j L 7j - K 

In summary, equations (13.1141 ) and (13.1 181) show that the ra-fluid FL cosmologies form an ra-parameter 
family, labelled by a conformal factor A and n — 1 mass parameters mi, . . . , m n ^\. We will refer to these 
parameters as the intrinsic parameters. 



"Equations ( 13 . 1 1 4-t and < |3 . 1 19t illustrate the well-known scale-invariance of Einstein's field equations: if gtj is a solution, 
then so is §ij — \~ 2 gij for any constant A, with the matter terms scaling appropriately, e.g. p, = A 2 /i. 
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3.2.2 The Observational Parameters 

In this section we introduce the observational parameters associated with ra-fluid FL cosmologies, based 
on the current epoch to, which represents the age of the universe. We denote the values of H, q, f2j and 
Q k at t by 

Ho = H(to), qo = q(to), fi^Q = Oj(to), ^fc,o = ^fc(*o)> 
where Hq is the Hubble constant. Then, the set of constants 

{t ,H ,qo,0, Afi ,n k) o} , A = l,...,n, (3.122) 

are collectively referred to as the observational parameters (see, for example, Wainwright and Ellis 
(1997), p55-8). 

In addition, we require the parameter £q, which denotes the length scale at the present time, i.e £q = 
£(to). This quantity is not an observational parameter, but if K ^ 0, £$ can be determined in terms of 
observational parameters using ( 13.351) and ( 13.361) , i.e. via 

£2 = ~ K = ~ K n (3.123) 

Hfakfl H^il - Ya=i ^i,o) 

In the case of zero curvature (K = 0), £q can be chosen arbitrarily in order to fix the multiplicative factor 
that determines £ (also see section |3~.l.l| ). 

We now present the observational formulation of the n-fluid FL cosmologies. In this formulation, the 
dimensionless length scale and dimensionless time variable are given by 

a = —, and T = H t. (3.124) 

£q 

From (13.881) we have 

IH = 3Hfo ifi a- 3 ^. (3.125) 
Using ( 13.1241 ) and ( 13.1251 ), the Friedmann equation ( 13.141 ) assumes the form 

= 5> >0 a-^+ 2 + (3.126) 
' i=i 

We use uniqueness of the solution to ( 13. 126b under some appropriate initial conditions (see section IC.2I ) 
to conclude that in an expanding epoch the length scale is given by 

£ = e Q a(H t,n l70 ,...,n nfi ). (3.127) 



3.2. Parameterization ofFL Models 



46 



It follows from ( f379b . ( 13.1241) and (13.1261 ) that the Hubble scalar assumes the fornix 



H 2 = Hi 



k,0 a 



i=l 



(3.128) 



and further that the density parameters ( 13.311 ) and ( 13.451 ) assume the form 

a,oa- 37 ' +2 



E-=i%,oa- 3 ^ +2 + ^, ' 



and 



Xi 



1^ g- 3 ^+ 2 



(3.129) 



(3.130) 



The parameters Qk,o, Qo, to in (13.1221 ) can be related to the parameters {fi^o, Ho}, as follows: First, 
Qk,o is determined in terms of Q^q via the Friedmann equation (13.371) : 



fc,0 



Second, go is determined in terms of f^o via (13.421 ): 

n 



(3.131) 



(3.132) 



8=1 



Third, knowing i?o, we can obtain to by integrating (13.1261) as 



1 



o a/FK^o) 



da, 



where 



F(y, Q it0 ) = 1 + ^,o(a 2 ~ 371 - 1) 



(3.133) 



(3.134) 



and Jlfc o nas t> een eliminated using (13- 131b - As a consequence, the only independent observational pa- 
rameters are the n + 1 quantities defined by 



(3.135) 



12 An expression of this form is given by Peacock (1999), eq. (3.8). 
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It is useful to note that the normalized scale factor can be expressed in terms of the redshift z, defined 

by 

l + z = ^ = -, (3.136) 

(see Wainwright and Ellis (1997), eq. (2.41)). One can thus use (13.1281) and ( 13.1291) to express H and fi 4 
in terms of z. 

We now present the link between the observational parameters and the intrinsic parameters. On sub- 
stituting m = 3H 2 Qi (see (13.311 )) into ( 13.1091 ) at t = t , we obtain 

Ar 1 = t^^Hlll)^. (3.137) 

3.2.3 Conserved Quantities 

We now show that the evolution equations ( 13.571 ) and ( 13.60b admit n — 1 conserved quantities. 
From (13.451 ) and (13.1141 ) we have that 

Xl = rriiL-^. (3.138) 

We now eliminate A 2 //i and L by taking products of powers of three distinct matter variables (we require 
n > 3). Indeed, it follows from (13.1381 ) that 

xr ik xf~ ll xT l3 = mf-^mf-^ml^. (3.139) 

In the case of K ^ 0, we also obtain a conserved quantity describing the curvature of a model. We use 
(13.91 ) and (13.471 ) to rewrite the Friedmann equation (13.141 ) as 

II - tan 2 n| = (3.140) 

Then from (13.1381) and (13. 140b . we have 

3%' -2 _3 7i 4-2 M 2013(7,-7,) 3 7j-2 -37^+2 cita-w 

Xi Xj |I — tan ll\ Us hJ =m i m- . (3.141) 

Since there can be only n — 1 independent conserved quantities, we must fix some of the values of i, 
j and k in ( 13.1391 ) and ( 13.1411 ) to obtain a minimal independent set. We choose i = 1, k = n in ( 13. 139b 
(with 2 < j < n — 1) and i = 1 and j = n in (13.141b and so defined 

Mjixu ■ ■ ■ ,Xn) = xl^X?- 7 ^^, (3-142) 



3 Observe that if we choose j = 1 or j — n in < I3. 142b we simply obtain Mi = 1 and 7Vf n = 1. 
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and 

3 7 n-2 (371-2) 

/C( X i, Xn, n) = xF^Xn (1 - tan 2 O). (3.143) 

Upon recalling that m n = 1, the conserved quantities are then determined by the intrinsic parameters 
mi, . . . , m n -i according to 

^( Xl ,..., X n) = m?~ 7 X n ~ 7l > j = 2,...,n-l, (3.144) 

and 

37n-2 

/C(xi,Xn,fi) = ^r^ 5 • (3-145) 

Comment: We have derived these conserved quantities algebraically, and it follows immediately from 
(131441 and (l3~T45l) that 

-7^ = 0, -jr = 0. (3.146) 

It involves a lengthy calculation to verify (13.1461) using the ODEs (13.631 ) and (13.64 1 ). 

If K = (i.e. £1 = j) then K, = 0. As described in section |3~.2.2I we can use the freedom to rescale 
£ to set 

Ai = X n , (3.147) 

which implies m\ = 1. Then (13.1441 ) gives 

Mj = mj" -71 , j = 2, . . . , n - 1. (3.148) 

The conserved quantities (13.1421 ) and (13.1431 ) can also be expressed in terms of the density parameters, 
as follows: 

M 3 (ni,...,n n ) = nf^nf-^nl 1 ' 7 ' , j = 2,...,n-i, (3.149) 

3 7 n-2 37i-2 

/ccni.fin.fi*) = -fi 1 3(7i - 7 - ) n n 3(7i - 7n) n fc , (3.150) 

The conserved quantities are thus related to the observational parameters by 

M J (n 1 ,...,n n ) = ^f ln) ^%"- ll) ^ J \ (3-151) 

and 

3 7 n-2 (371-2) 

£(ni,n n ,n fc ) = -n^-^'nn.^^nfc.o, (3.152) 

upon evaluating ( 13.1491 ) and (13.1501 ) at t = to. We note that conserved quantities of this type have been 
previously discovered, for example, by Ehlers and Rindler (1989) for RDCA FL universes, and more 
generally by Lake (2006) for n-fluid FL universes, using a different method. 



3.3. 2-Fluid FL Cosmologies - Qualitative Analysis 



49 



3.3 2-Fluid FL Cosmologies - Qualitative Analysis 

We now use the results of sections |3~T1 and [3721 to give a qualitative analysis of the dynamics of 2-fluid FL 
cosmologies. The results of this section, in particular figures |3~TT|3.3I are contained in Wainwright (1996), 
although the use of the conserved quantity /C is new. 

3.3.1 General Features 

The state space is two-dimensional, being described by the variables (O, xi, X2) subject to the restrictions 
(13.461) . We use (13.461) to eliminate xi vra 

Xi = l-X2- (3-153) 

The state space is then the rectangle in the (Cl, X2)-plane given by 

- vr/2 < < vr/2, and < %2 < 1- (3.154) 

In accordance with theorem 13.1.11 there are two qualitatively different classes of 2-fluid FL cos- 
mologies, depending on whether 72 > | or 72 < | (recall that 71 and 72 are assumed to satisfy 
< 72 < 71 < 2). We will discuss these two classes in sections 13.3.21 and 13.3.31 using the two most 
important examples from a physical point of view0 

i) Universes with radiation and dust (i?DC-universes, 71 = |, 72 = 1), 

ii) Universes with dust and cosmological constant (DCA-universes, 71 = 1, 72 = 0). 

An important difference between these two cases is that if 72 < | the Einstein static universe, repre- 
sented by an equilibrium point E, is included in the family of models. The point E is given bj^l 

X2,E= J 7l ~ 2 v « = 0, K = +l, (3.155) 
3(71 - 72) 

(see equation ( I3.74I )). A second related difference is the behaviour of the compact deceleration parameter q 
as given by (13.65 b - If 72 > I then q is positive on the state space (13.621 ) and all models are decelerating. If 



14 Other choices are qualitatively similar 

15 The restriction 72 < I ensures that X2,e < 1, as required. 
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72 < | then q is negative on a subset of the state space. Figure l3T2l shows that the orbits of ever-expanding 
models eventually enter this subset, and as a result the models are accelerating at late times. 

The portraits of the orbits are drawn by first sketching the various invariant sets and equilibrium points, 
as given in section |3~. 1.5 1 in particular the boundary of the state space (13.1541 ) and the flat FL invariant set 
Q = ±|. In the case 72 < |, the orbits that are past and future asymptotic to the Einstein static point 
E should be drawn. At this stage, the orbits that have been drawn represent exceptional models and form 
a "skeleton" for the state space. The orbits of typical models can now be drawn by noting the past and 
future attractors as given in theorem 13.1.31 In drawing the portraits it is helpful to note that the evolution 
equations (13.631) and (13.641) are invariant under the interchange 

(n,f)-(-n,-f), (3.156) 

i.e. changing the sign of Cl reverses the direction of time. In particular, the invariant set 
— § < < — f represents cosmologies that evolve in the same way as those with j < Cl < | , with the 
difference that the direction of time is reversed. These models are contracting throughout their evolution 
(Q < implies H < 0) and hence are not potential models of the real universe. 

In practice the state space can be sketched numericallj^l using the conserved quantity /C that is defined 
by (13.1431 ) for the family of 2-fluid cosmologies. On substituting (13.1531 ) into (13.1551 ) and choosing n = 2, 
we obtain 

3 72 -2 ( 3 -n - 2 > 

X2) = (1 - X2) X 2 3(71 " 72) (1 - tan 2 fi). (3.157) 

The orbits are given by 

JC(n, X2) = constant. (3.158) 
The value of /C at E is denoted /C e and is given by 

372-2 (371-2) 

JCe = {1-X2,e) 3 ^-^{X2,e) 3( ">. (3.159) 

It follows that the equation 

K(n, X 2)=JC E (3-160) 

describes the orbits that are past and future asymptotic to E, i.e. it describes the stable and unstable 
manifolds of E. 

The value of K, determines whether a model is ever-expanding or not, as shown in the following table. 



'For example, using the plots [ implicitplot ] command in Maple. 
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Case 72 > 2 



Restriction on /C Qualitative Behaviour 

2 

3 



— oo < K, < Ever-expanding 
< K, < +oo Recollapsing 

Case 72 < § 

—oo<K,<K,e Ever-expanding 

Ke < K. < +oo Recollapsing or bouncing 

Further, the sign of K, determines whether the model is open, flat or closed. We note that all generic classes 
of FL models are covered by the given values of K. 

In the case 72 > | the critical value of K, is fC C rit = 0, indicating that the flat FL models play an 
exceptional role as regards to future evolution. In the case 72 < |, JC cr it = K-E and the exceptional 
models are those that are past or future asymptotic to the Einstein static universe. This distinction is 
clearly indicated by comparing figure |3"TTl (-R-PC-universes) with figure |3T2l (£>C A- universes). 

We note that the conserved quantity K, can also be expressed in terms of the density parameters. It 
follows from ( 13.1501 ) with n = 2 that 

372-2 _ (371-2) 

/c(o l5 o 2 , n fc ) = -ni^-* 2) n 2 ^^n k , (3.i6i) 

with 

Q k = 1 - Six - Q 2 - (3.162) 
3.3.2 FL Cosmologies with Radiation and Dust (i?L>C-universes) 

We now consider 2-fluid FL cosmologies containing radiation (71 = |) and dust (72 = 1). The portrait of 
the orbits, shown in figure l3TT1 is drawn as described in the previous subsection. 

For the i?DC-universes, the expression ( 13.1571 ) for the conserved quantity /C becomes 

K = ^(1 -tan 2 ^). (3.163) 

Equivalently, the expression (13.1611 ) gives /C in terms of the density parameters as 

K = -^. (3.164) 
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FIGURE 3.1: The (T2, xi) state space for 2-fluid FL cosmologies containing radiation (71 = |j and dust 
(72 = 1), showing regions of expansion (H > 0, i.e. £1 > 0) and regions of contraction (H < 0, i.e. 
n < 0). 

The value of JC is infinite or zero on the skeleton of the state space: 

K, = +00 if xd = Oand < \, 

/C = if Xd = 1 or = f , 

/C = -oo if ft = ±§,orxd = Oand > f. 

Note that /C is indeterminate at the past attractor A~ = (xd = 0, fl = \ ) and at the future attractor 
A + = F r ~ U M^" (xd = 0, ft = j ; = 1, ft = | ). This indeterminacy arises since infinitely many 
orbits meet at and A + . 

Discussion: Typical models differ qualitatively depending on the value of /C, which determines where 
the orbit lies in the state space. Firstly, for ever-expanding models if the orbit passes close to M+, the 
dust component will never be significant dynamically (i.e. ft^ <C 1) and this will correspond to K, being 
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large and negative. On the other hand if the orbit passes close to Ft then will attain a maximum value 
close to 1, corresponding to /C being negative and close to zero. Secondly, recollapsing models differ 
quantitatively by how close the orbit comes to Ft, i.e. by the maximum value of Xd> which in turn will 
depend on K,. One can determine the dependence on JC explicitly, as follows: 

i) For ever-expanding models, we obtain 

n ^ = 1 + 2 1 > /zjc ' /C<0 - (3 " 165) 

ii) For recollapsing models an expression for Xr\min and Xd\max can be obtained by setting = in 
(13.1631 ) and solving the resulting second degree polynomial. We obtain 



2 

i + vT+ 4/c ' 



Xd\ m ax ~ -, . n — ; — n^> Xr\ m i n — 1 Xd\ m ax ' /C > 0. (3.166) 



3.3.3 FL Cosmologies with Dust and Cosmological Constant 

We now consider 2-fluid FL cosmologies containing radiation (71 = |) and dust (72 = 1). The portrait of 
the orbits, shown in figure [3^21 is drawn as described in section IT. 3. II 

For the DCA-universes, the expression ( 13.1571 ) for the conserved quantity K becomes 

K(Cl, X a) = X d 2 '\l l, \l ~ tan 2 f2). (3.167) 
Equivalently, the expression (13.1611 ) gives K, in terms of the density parameters as 

k = -n^ 2/3 n A 1/3 n k . (3.168) 

The critical value of xa, the value of xa at the Einstein static equilibrium point, is given by (13.1551 ) as 

Xa,e = i (3.169) 
It then follows from ( 13.1591 ) that the critical value of K, is 

JC E ={±)- l/ \ (3.170) 
The value of K, is infinite or zero on the skeleton of the state space: 
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K 
K 
K 



+00 if \Q\ < j and xa 



■ or 1, 



: if Q 
■ — 00 if (l 



±- 



:| or |Q| > f and xa = or 1. 



Note that K. is indeterminate at the past attractor A 
attractor ^4+ = deS+ U F~ (xa = 1, & = f ; xa = 0, 



F+ (xa = 0, fi = §) and at the future 

4^ 



Contracting 
H<0 



Xa 



Expanding 
H>0 




Flat FL 
K = 



Accelerating 
q<0 



Decelerating 
q>0 

£2 



Figure 3.2: The (£1, X2) state space for 2-fluid FL cosmologies containing dust (71 = 1) and cosmolog- 
ical constant (72 = 0), showing regions of expansion (H > 0, i.e. Cl > 0), regions of contraction (H < 0, 
i.e. £1 < 0), regions of acceleration (q < 0) and regions of deceleration (q > 0) in relation to the Einstein 
static solution E, which satisfies H = and q = 0. 
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FIGURE 3.3: The (Cl, x) representation of the state space for the FL models with matter, a cosmological 
constant and spatial curvature. The thumbnail plots depict the length scale £(t) of FL models correspond- 
ing to orbits in each region of the state space and to each equilibrium point. Cases (1 ), (2) and (3) represent 
the three generic families of FL models, i.e. their orbits fill an open subset of the state space (2-parameter 
families). They represent, respectively, ever-expanding models, recollapsing models and bouncing mod- 
els. On the other hand, cases (4), (5), (6a) and (6b) correspond to a single orbit in the state space (a 
1-parameter family), while the cases in the third row represent the unique solutions corresponding to the 
equilibrium points. 
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Discussion: Firstly, ever-expanding models differ qualitatively depending on the magnitude and the sign 
of the curvature scalar Qk- If the orbit stays close to flat FL (O = the curvature will never be significant 
dynamically (i.e. <C 1) and this will correspond to JC being close to zero. Secondly, recollapsing 
models (and bouncing models) differ qualitatively by how close the orbit comes to the Einstein static 
model, i.e. by the maximal (or minimal) value of xa- Since xa < XA,E in recollapsing models (see 
(13.1551) ). the cosmological constant will never be the dominant fluid. Similarly, since xa > XA,E hi 
bouncing models, there will be no dust dominant epoch, i.e. Xd ~ 1- One can determine the dependence 
on K, explicitly, as follows: 

i) For ever-expanding models, we obtain 

I^La* = F^nr' £<Ke. (3.171) 

ii) For recollapsing and bouncing models, an expression for XA\max an d XA\min can be obtained by 
setting Cl = in (13.1671 ) and solving the resulting third degree polynomial. 



3.4 3-Fluid FL Cosmologies - Qualitative Analysis 

We now use the results of sections I3.1H3.3I to give a qualitative analysis of the dynamics of 3-fluid FL 
cosmologies in which the matter content is radiation, dust and a cosmological constant (A > 0) (labelling 
the x-i as Xr, Xd and xa)- The analysis given in this section is new. 

3.4.1 General Features 

The state space is three-dimensional, being described by the variables (Q, x-r, Xd, Xa) subject to the 
restriction (13.461) . We use (13.461) to eliminate Xr via 

Xr = l-Xd-XA- (3.172) 

The physical state space is then a solid triangular prism in M 3 given by 

-|<^<|, 0<Xd,XA<l, Xd + XA<1- (3.173) 

We can construct a "skeleton" for the state space using various invariant subsets, as shown in figure 
13.41 The choice of Xi = for any one fluid gives a 2-fluid FL cosmology, identified as follows: 
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Xa = i?DC-universes, 
Xd = i?CA-uni verses, 
Xr = L>CA-uni verses. 

Further, the orbits in the planes £1 = ±j represent the flat 3-fluid universes. 

The Einstein static universe is represented by a line segment of equilibrium points denoted E r E<i, 
which is the intersection of the planes £1 = and q = 0. The stable manifold S + of the set of Einstein 
static points is a surface through the line segment E r E^ terminating at deS~ when xa = 1 and at the line 
segment F^F^ when xa = 0. This surface can be described by an equation of the form 

& = f(xd,XA), (3-174) 

with Xd and xa restricted as in (13.173I ). Because of the symmetry (O, f ) — > (— f2, — f ) (see (13.156I )). the 
unstable manifold S~ is the surface 

fi = -/(XAXA). (3.175) 

We know from theorem D.l.ll that there are three generic classes of RDCK FL universes. Their orbits 
are contained in invariant subsets bounded by the surfaces S^, as follows: 

Ever-expanding universes Q, > \f(xd, Xa)\ (the subset to the right of S + and S~), 
Recollapsing universes (the subset between S + and S~ and below the plane q = 0), 
Bouncing universes (the subset between S + and S~ and above the plane q = 0). 

The planes Cl = and q = 0, although not invariant sets, play an important role in describing the 
dynamics. The plane Cl = divides the state space into a region of expansion (Cl > 0) and a region of 
contraction (Cl < 0). The variable q is given by (13.651 ) which here specializes to 

Q = Xr + \xd - XA- (3.176) 

On using ( 13.1721 ) to eliminate Xn we find that q = on the plane 

Xd + 4XA = 2. (3.177) 

Equation ( 13.1771 ) then divides the state space into a region of acceleration (q < 0) and a region of deceler- 
ation (q > 0) as follows: 
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Region of acceleration xa > i(2 — Xd), 
Region of deceleration xa < z(2 — Xd)- 

It follows from the general discussion in section 13.2.31 that there are two conserved quantities for the 
RDCA FL universes, defined by 

KA ( XrXA\ 1/3 , y 1- tan 2 ft 

.M = — — , and /C = — j-. (3.178) 

V Xrf / {XrXA) l/2 

Using (13.1721) we can eliminate Xr from these equations and hence obtain J\A and /C in terms of ft, Xd and 
Xa- The values of M. and /C at the Einstein static equilibrium points E are denoted M. e and KZe and are 
obtained on substituting (13.1771) into (13.1781) . We have 



M 



(3 X A - 1) 3 XA 



1/3 



< M E < oo, (3.179) 



16(1 -2xa) 4 . 

/C B = [(3xa-1)xa]~ 1/2 , 2</C £ <oo. (3.180) 

The orbits of typical models are described by the equations 

K.(Cl, Xd, Xa) = constant, and M(xd, Xa) = constant, (3.181) 

i.e. the orbits are the intersections of these two families of surfaces in the state space. As a special case 
the flat universes satisfy K = {i.e. = |) and then the equation 

M(Q,Xd,XA) = constant (3.182) 
describes the orbits in the invariant set ft = f and can be used to numerically sketch them, as in figure 



Discussion: Knowing the past attractor A~ = F 7 f and the future attractor A + = F~ U deS + one can 
visualize the orbits of typical universes. For example, the orbits of ever-expanding models join F r + to 
deS + . The models will differ qualitatively depending on whether they come close to Ft (f^ ~ 1 in 
a neighbourhood of F^) and whether they come close to the Milne set (Cl = ^) or to the surfaces 5 ± , 
in which case there will be an epoch in which spatial curvature is dominant (i.e. |Ofc| « 1). One can 
determine the points in state space at which | | and Xd attain a maximum value, as follows. If K ^ 0, it 
follows from (13.361 ) and ( 13.691 ) that Cl and \Qk\ are simultaneously extremal. Since in an ever-expanding 
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model Q is extremal when q = 0, we conclude that is extremal on the plane given by ( I3.177I ). The 
matter variable Xd will attain a maximum when x'd = 0> which implies g = |, as follows from (13.641) and 
( 13.651 ). After substituting ( 13.1721 ) to eliminate the Xr term in ( 13.1761 ), we conclude that Xd will attain a 
maximum on the plane 

Xd + 4xA = l. (3-183) 

One can then determine the dependence of Xd\max on M. explicitly from (13.1721) . (13.1781 ) and (13- 1 83b . as 
follows: 

3 3/4 

Xdlmax ~ 4^3/4 + 3 3/4 ' (3.184) 

3.4.2 A Model of the Physical Universe 

In this section we discuss the RDCK FL universes as viable models for the large scale dynamical be- 
haviour of the physical universe. Current observations of distant galaxies, type I a supernovae and the 
cosmic microwave background radiation support the following conclusions: 

i) The universe is expanding in the present epoch (Hq > 0). 

ii) The universe is accelerating in the present epoch (qo < 0). 
hi) The early universe was radiation-dominated (hot big-bang). 

iv) Pressure-free matter (baryonic matter and dark matter) is dynamically significant in the present 
epoch. 

v) The spatial geometry is close to flatness in the present epoch. 

These conclusions suggest modelling the physical universe as an RDCA FL cosmology, with obser- 
vational parameters 

&r,o, ^d,o, ^fc,o> ^A,o> Qo, Hq, and to- (3.185) 

Equations (13- I31b - (I3- 133b limit the number of independent quantities in this set to four, as expected for 
an RDCK model. Values for these parameters have been determined using observations, such as high 
redshift galaxy surveys (for example, by 2dF and SDSS; see Spergel et al. (2006) for detailed references) 
and analyses of the power spectrum of the cosmic microwave background (for example, by WMAP). We 
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FIGURE 3.4: The skeleton of the state space for RDCA cosmologies, using Xd ar >d Xa as variables, 
with Xr = 1 - Xd ~ XA- 



Xd+4X A =1 




Xd+4X A =2 



FIGURE 3.5: The invariant setfl = \ in the state space for the RDCA universes, describing the expand- 
ing flat FL models. 
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refer the reader to Spergel et al. (2006) and Tegmark et al. (2004) for an analysis of recent observational 
data, leading to the following values: 



H = 71.0 ± 4Mm s' 1 Mpc' 1 , (3.186) 

n rfi = (4.9 ±0.5) x 10" 5 , (3.187) 

n dfl = 0.27 ±0.03, (3.188) 

O Ai0 = 0.73 ±0.04, (3.189) 

n kfi = -0.010 ±0.014. (3.190) 

Using the observational values (13. 1 86b -( f37T90l . one can apply d3.63N3.65l l and d3.69N3.71l > to nu- 
merically obtain the past and future behaviour of all matter quantities \i ar *d the curvature indicator O. We 
present the results of one such numerical simulation in figure [3761 along with computed epochs of note. 

The conserved quantities in the RDCA universe, as given by (13- 178b . can be expressed in terms of the 
density parameters as 

and 

K = s£*n^n k , (3.192) 

(see (13. 149b and d3.150b ). Then using d3.187b - d3.190D . we obtain 

M = (2.63 ± 0.07) x 10~ 3 , (3.193) 

and 

K = -1.72 ±1.91. (3.194) 

The values of these conserved quantities, which were fixed at the end of the inflationary epoch, determine 
a small subset of RDCA universes that can potentially describe the real universe. The value of M. 
ensures the occurrence of a matter-dominated epoch that is sufficiently long, but not too long (see equation 
d3.184b ). The value of K, ensures that is close to in the past and also into the future. 
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Figure 3.6: Past and future evolution of the RDCA model of the physical universe, depicting the matter 
variables (top), Hubble-normalized energy density (bottom) and epochs of primary interest. 
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3.5 Explicit FL Cosmologies 

In this section we give a unified description of the known solutions of the Friedmann DE that can be 
expressed in terms of elementary functions, for 1 -fluid and 2-fluid FL cosmologies. In these cases, the 
Friedmann DE can be transformed to the form 



-^j =Q(C,m,K), (3.195) 

where Q is a quadratic function in C, m is a dimensionless mass parameter and K is the curvature 
indicator. Here £ is a function of the length scale factor i and the variable T is either proportional to 
conformal time r\ (in the case A = 0) or proportional to clock time t (in the case A > 0). 

Recall that the n-fluid Friedmann equation can be written as a DE of the form (13.1 101 ), reproduced 
here for convenience: 

(^ 2 = E(^r 37l+2 -^ (3-196) 

^ ' 1=1 

We are primarily interested in expanding models, which begin at a big-bang singularity, i.e. the length 
scale satisfies £(t s ) = at some time t s , and > in some interval. Since the DE ( 13.1961) is autonomous, 
it is invariant under the change t — > t + C. We will use this freedom to set t s = 0, thereby fixing the 
constant of integration that arises in solving ( 13.1961) . We will thus solve ( 13.1961) subject to the requirement 
% > and the initial condition 

£(0) = 0. (3.197) 
Using conformal time (ID. lb . equation ( 13.1961 ) becomes a first order ODE for £{rj): 

^(^) 2 = twr 3 ^ 2 -K. (3.198) 

The solution £(tj) can in principle be expressed in terms of clock time t using the equation (ID. lb - As with 
( 13.1961 ), when solving ( 13.1981 ) for £(r/), we will impose the requirement ^ > and the initial condition 

£(0) = 0. (3.199) 

In the following sections we consider five classes of models in which either ( 13.1961 ) or ( 13.1981 ) can 
be transformed into ( 13.1951) . The solution of ( 13. 195b subject to the requirement ^ > and the initial 
condition £(0) = is given in Appendix In each case, we give the solution for £{t) or £(tj) and note 
that the corresponding expressions for the energy densities and pressure are given by equation ( 13.151 ), 
(I3T71) and (l3~TT4l) . 
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3.5.1 Single Fluid FL Cosmologies 



In this case we consider solutions to the Friedmann equation with n = 1, § < 7 < 2 and arbitrary curva- 
ture. We choose the conformal factor to be A = Ai. The Friedmann DE (13.1981 ) assumes the form 

2 



1 



We perform the change of variable 



with 



I 2 \dr) 
C = (A^) 1//3 




/3 ! 



3 7 



so as to obtain 



fdC 



- 



1 - KC\ 



Then from (IC.3b and (IC.4b . the solution is 



C = S K (T), 



where 



Sk(T) 



sinh(T), 
T, 

k sin(T), 



openFL (K = — 1), 
flatFL (if = 0), 
closed FL (K = +1). 



Thus, 



l{rj) = A 

Using ( 13.1 14b , the energy density is given by 

H = 3A 2 L- 37 



13 



L = XL 



(3.200) 



(3.201) 



(3.202) 



(3.203) 



(3.204) 



(3.205) 



(3.206) 



(3.207) 



This solution, for arbitrary 7, was first given by Harrison (1967) (see equations (152), (18) and (32) 
in Harrison (1967), for open, flat and closed models, respectively). Specific cases were discovered earlier, 
however. Friedmann (1922) gave this solution in the case of dust and positive curvature {i.e. 7 = 1 and 
K = +1) and Tolman (1931) gave this solution in the case of radiation and positive curvature {i.e. 7 = | 
and K = +1). 
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These solutions describe two qualitatively different classes of models: ever-expanding models if K = 
— 1 or 0, and recollapsing models if K = +1. In this case, it follows from (13.2051) and (13.2061) that 
recollapse (i.e. £ — ► as rf — ► 77/) will occur at 

77/ = pn. (3.208) 

For flat FL (K = 0) we can obtain £ = £(t) explicitly by integrating ( 13.1961 ) with n = 1 and K = 0. 
The result is 



2 

(3.209) 



The parameter A is not essential, since by rescaling I it can be assigned any value. Equation ( 13.209 1 ) thus 
represents the unique flat FL universe with a single fluidf^l 



In the case of radiation (7 = |, (3 = 1) one can integrate (ID. 1 b to obtain, 

A -1 (cosh 77 — 1), openFL (K = — 1), 

(3.210) 

A _1 (l - cost?), closed FL (K = +1). 

Substituting back into (13.2041 ) and rearranging yields 

£(t) = X- 1 [2Xt- K(Xt) 2 ] 1/2 . (3.211) 

We note that the second term here represents the contribution to the solution from spatial curvature, i.e. 
the drift from flat FL. If K = — 1, the limit A — > 00, gives the Milne solution (given in section [3. 1.21 
equation ( 13.281 )). 

In the case of dust (7 = 1, f3 = 2) one can integrate ( ID. II ) in closed form, but cannot invert the 
expression in terms of elementary functions. In this case t and 77 are related by 

|A _1 (sinh7? — 77), openFL (K = — 1), 

. (3.212) 
ttA (77 — sin 77), closed FL {K = +1). 



"That (13.209b does define a unique solution is confirmed by the fact that the physical quantities H and /1 do not depend on A, 
being given by 



ff=-t , and M =— i 

on using d3.9t , {332} and the fact that Q = l. 
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3.5.2 FL Vacuum Cosmologies 

In this case we solve the Friedmann equation where the only fluid component is a cosmological constant, 
i.e. n = 1, 7i = 0. We choose the conformal factor to be A = Ai, which leads to 



A = J±A, (3.213) 

(see equation ( 13. 11 11 1). The Friedmann equation (13. 1961 ) in this case is 

d£\ 2 

— J =(\lf-K. (3.214) 



We perform the change of variable 



so as to obtain 



C = \£, T = Xt, (3.215) 



2 



C? - K. (3.216) 



dTj 

This DE can only be solved subject to £(0) = and j& > in the case of open FL, with the solution 
given by (1C.3I ). For the case of non-negative curvature, we instead require that £(0) = 1 and that C is 
increasing for all T > 0. It follows that (13.2161 ) has a unique solution in each case, leading to 

A" 1 sinh(At), open FL (K = -1), 
£(t) = < A- 1 exp(At), flatFL (if = 0), (3.217) 
A" 1 cosh (At) , closed FL (K = +1), 
where A is given by (I3.213I ). The energy density is constant and given by 

fi = 3A 2 = A. (3.218) 

These solutions were first given by Robertson (1933), de Sitter (1917) and Lanczos (1922) in the case 
of negative, zero and positive curvature, respectively. They correspond to the orbits M + — ► deS + for 
K = — 1, deS + for K = and deS~ — > deS + for K = +1, along the 1-fluid cosmological constant 
manifold in figure 13.31 In the case of flat FL, (13.2171 ) gives the de Sitter solution, described in section 
section [3~L2] equation (13.271 ). 
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Often, in the case of the de Sitter solution, it is preferrable to use conformal time instead of clock time. 
In particular, when K = this choice of time variable is defined on an interval (—00,77^) (see Appendix 
IdT >. Upon defining 

fj = Vf -r], (3.219) 
the length scale ( 13.2171 ) for a flat background is simply 

i(v) = (3-220) 
3.5.3 Flat FL with a Single Fluid and A 

In this case we solve the Friedmann equation with n = 2, 71 = 7 arbitrary, 72 = and zero curvature (i.e. 
K = 0). We choose the conformal factor to be A = A2, which leads to d3.213l >. As described in section 
13.2.21 we can use the freedom to scale £ to set Ai = A. The Friedmann equation (13.1961 ) then takes the 
form 

-) =(\£r^+ 2 + (X£) 2 . (3.221) 
at J 



We apply the change of variable 
so as to write ( 13.1961 ) as 

Then using dC3l ), we obtain 
and so 



C = (X£)^, T = | 7 Ai (3.222) 
^) 2 = l + £ 2 - (3-223) 



£ = sinhT, (3.224) 



sijili ( ~^-^t 



2 

■■'•-1 



(3.225) 



£{t) = \- 1 

Using (13.1141) . the energy density is given by 

H = H f + H A = 3A 2 L~ 37 + A, L = XI. (3.226) 

These solutions were first given by Harrison (1967) in equation (23). They represent ever-expanding 
models and correspond to the orbit F + — > deS + along the K = invariant set, as in figure l3.3l 
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3.5.4 FL Cosmologies with Radiation and Dust (i?L>C-Universes) 

In this case we solve the Friedmann equation with n = 2, 71 = | and 72 = 1. We choose A = A2 as the 
conformal parameter. Then, by (13 - 1 13b . the single mass parameter is given by 

m = m r = I — ■ (3.227) 



A 



The Friedmann equation (13.1981 ) assumes the form 



We apply the change of variable 

C =-=()£), T = V, (3-229) 

Jm 



(3.230) 



which gives 
Using (1C.3I) and (IC.4I) . the solution is 

C = S K (T) + ^-S K {\T) 2 , (3.231) 



m 1 



and so 

l(n) = A" 1 (y/mS K {ri) + S K {\r,) 2 ) , (3.232) 
where Sk is given by (13.2051) . Using (13- 1 14b . the energy density is given by 

H = n r + n d = 3\ 2 mL~ 4 + 3A 2 L~ 3 , L = XL (3.233) 

In this case, one can integrate (ID. II ) so as to obtain an expression for the clock time t = t(rj), but cannot 
invert this expression to obtain I = £(t): 

A -1 (\/m(cosri r\ — 1) + ^(sinhr/ — 77)) , open [K = —1), 

A" 1 + T2V 3 ) , Sat (K = 0), (3.234) 

A -1 (\/m(l — cosr]) + 5(77 — sinjy)) , closed (K = +1). 

The solution (13.2321) contains the one-fluid dust and radiation solutions given in (13.2061 ) as special cases. 
One can see this by writing (13.2321 ) in terms of parameters A r = Ai and A^ = A2, 

i(v) = K'Sxiv) + X^SK^rj) 2 , (3.235) 



t = < 
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and letting A r —* oo or — > oo. 

These solutions were first giver0by Chernin (1966). They describe all orbits in the state space for 
iiDC-universes with H > in some epoch (see figure IBTTT ) and hence describe two qualitatively different 
classes of models: ever-expanding models if K = —1 (the orbits F x — > M + ) or if = (the orbit 
F x — ► Fg") and recollapsing models if if = +1 (the orbits F x + — > i^ - ). The dimensionless parameter 
m determines key physical properties of the models. First, it determines the epoch of matter-radiation 
equality (fi r = ^ according to 

L eq = m, (3.236) 

as follows from (13.2331) . Second, for recollapsing models, m determines the time of recollapse, according 
to 

r]f = 2(tt — arctan(2y^)), (3.237) 

as follows from ( I3.232I ). 

3.5.5 FL Cosmologies with Radiation and A (i?CA-Universes) 

In this case we solve the Friedmann equation with n = 2, 71 = |, 72 = and arbitrary curvature. We 
choose the conformal factor to be A = A2 so that ( 13.2131 ) holds. Then, by ( 13.1131 ), the single mass 
parameter is given by 

V- 2 



m = mx = I _± 1 . (3.238) 
The Friedmann equation (13.1961) then takes the form 

= m{X£)~ 2 + {Xlf - K. (3.239) 



dt 

Upon making a change of variable according to 



C = -^=(\£) 2 , T = 2Xt, (3.240) 



we obtain 

' -.KC + C 2 . (3.241) 



icy 1 



dT J y^m 



Harrison reports that in the case K = +1 this solution was given by Lemaitre (1927), but we have not verified this reference. 
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Using dC.3b , the solution is 

C = sinhT - -Lifsmh 2 (±T), (3.242) 

and so 

£{t) = A" 1 [^msmh(2Xt) - Ksmh 2 (Xt)] 1/2 . (3.243) 

Using (13. 1 14b . the energy density is given by 

H = n r + n A = 3X 2 mL~ A + A, L = XL (3.244) 

The solution (13.2431 ) contains the CA vacuum solutions given in (13.2171 ) and the one-fluid radiation solu- 
tion given in ( 13.2111 ) as special cases. One can see this by writing ( 13.2431 ) in terms of parameters A r = Ai 
and Aa = A 2 , 

£( v ) = [A^ 1 A A 1 sinh(2A A t) - KX^ 2 smh 2 (A A i)] 1/2 , (3.245) 

and letting A r — > oo or Aa — > 0. Also, in the case of zero curvature (K = 0), ( 13.2431 ) reduces to ( 13.2251 ) 
with 7 = |. 

This solution was first given by Harrison (1967) (see equations (59) and (38), for the case of negative 
(K = — 1) and positive (K = +1) curvature, respectively). It describes all orbits in the state space 
for i?CA-universes with H > in some epoch (see figure 13.21 ) and an initial singularity. Hence, this 
solution describes both ever-expanding models and recollapsing models. The parameter m determines 
key physical properties of the models. First, it determines the epoch of radiation-cosmological constant 
equality (pi r = fj,&) according to 

L eq = mJ, (3.246) 

as follows from ( 13.2131 ) and ( I3.244I ). Second, for recollapsing models, m determines the time of recollapse 
according to 

t f = A~ 1 arctanh(2^). (3.247) 
This result follows upon rewriting (13.2431) using standard hyperbolic identities as 

£(t) = A" 1 [i sinh(2At) (2 - K tanh(At))] 1/2 . (3.248) 

We note that at the critical value of m, denoted m* and given by m* = |, the solution (13.2431) 
simplifies to 

£(t) = A" 1 - exp(-2At))] 3 . (3.249) 

This solution is future asymptotic to the Einstein static solution and corresponds to the orbit F + — ► E in 
figure [372] 
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3.5.6 Explicit Solutions in Terms of Observational Parameters 

In sections I3.5.1H3.531 we have presented the explicit solutions using intrinsic parameters, namely, a 
conformal factor A, the curvature indicator K and, for the 2-fluid solutions with non-zero curvature, a 
dimensionless mass parameter m, which describes the relative significance of the two fluids. 

In order to link with observations it is desirable to write the solutions in terms of the observational 
parameters. In order to do so we need to use the formula that relates the intrinsic parameters A, to the 
observational parameters, derived in section 13.2.21 In the case of K ^ 0, we can use equations ( 13. 1231 ) 
and (13.1371) to write 



-2 _ a I fi »,o N 371 



K °\MJ ' (3 ' 250) 
One can think of (13.2501 ) as relating A^ , the length scale determined by the i th fluid, to £q, the length scale 
determined by the spatial curvature. In particular, we will need the following special cases of (13.2501 ): 

-2 _ B 2 f ^r,0 \ D„Ai^: n „ /„, _ ^ 



\- z = f* , Radiation ( 7r = %), (3.251) 

X d 2 = l l{^fl) ' Dust( 7d = l), (3.252) 

Note that the Friedmann equation (13.36b . evaluated at to gives 

fi fc>0 = 1 - fi , (3-253) 

where Q,q is the total density parameter evaluated at to. In the case of a cosmological constant (i.e. 7 = 0), 
we can obtain an expression for Aa directly from ( 13.1371 ), as follows: 

A " = dw = ^' <7A = 0) ' <3 ' 254) 

In the case of K = and n > 1, we use the freedom to scale I in order to set Ai equal to A n . Then, as 
a result of (13.1371) . the length scale l§ is determined explicitly in terms of O^o and Ho according to 

(4#o) 6(7l-7n) = ^TX ( o 37l ~ 2) • (3.255) 

The intrinsic parameters Aj are then obtained via (13. 137b - 

There is no work to be done in writing the densities in observational form, since they are given in 
general by equation (13. 1 25b . which we repeat for convenience: 

Hi = ml^a^\ (3.256) 
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where a = £/£q. 

In order to illustrate this process we give the length scale function for the i?DC-universes in terms of 
the observational parameters. We first rewrite i[rj) in ( 13.2321 ) in terms of A r and A^: 



i{r ] ) = \- l S K (r ] ) + \- 1 S K {\rif. 
Then (13.2511 ) and (13.2521 ) immediately give 



m = to 



r.O 



II -fin 



1/2 



Sk(v) + 



d,0 



li - n 



(3.257) 



(3.258) 



gives the i?C-uni verse. 



where Q,q = Q r fi + fi^o- Then O rj o = gives the DC-universe and Q^q 
Peacock (1999) refers to the terms 

firf.O , fir,0 

2|1 — O | ' ^ 2|l-O | 
as the "dimensionless masses" (p79), when discussing RDC-universes. Note that these quantities are 
closely related to those given in (13.2511 ) and (13.2521) . upon applying (13.253I ). These quantities also appear 
in the DC- and RC-universes as given by Coles and Lucchin (1995) (see p39, eq (2.4.2) for dust in terms 
of 7] and p41, eq (2.5.2) for radiation in terms of t). 



3.6 Discussion 

We summarize the principal features of the dynamics of n-fluid FL universes. 

i) n-fluid FL universes have at most one epoch of expansion (H > 0) and at most one epoch of 
acceleration (q < 0, or equivalently, w < — |). These results depend crucially on the assumption 
that the fluids are non-interacting. We refer to Tolman (1934), p402 and p429-431, and to Clifton 
and Barrow (2007) for a discussion of models in which this assumption is not made. 

ii) On approach to a singularity (either initial or final, I — ► 0) the 1-fluid is dominant (i.e. xi —>■ 1) and 
the matter content dominates the spatial curvature (fi — ► 1, O.^ — > 0). In a typical ever-expanding 
model (£ — > oo) the n-fluid is dominant (\n — > 1), and the matter content dominates the spatial 
curvature if and only if j n < |. 

iii) The n-fluid FL cosmologies admit n — 1 conserved quantities formed from the density parameters 
0,1, ... , $7 n . In particular, they admit one curvature quantity and n — 2 matter quantities. 
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We now comment on an important difference between n-fiuid FL cosmologies with n > 2 and 2- 
fiuid cosmologies. For a 2-fluid FL model, there are two essential parameters, a conformal parameter 
A and a mass parameter mi, or equivalently, two observational parameters fii ; o and f^o- There is also 
one conserved quantity K. that keeps track of the spatial curvature. For a 3-fluid FL model there is an 
additional mass parameter 777-2, and hence three observational parameters fi^o, ^2,0 an d ^3,0- There is 
also a second conserved quantity M. The new dynamical feature that arises in generalizing from a 2-fluid 
model to a 3-fluid model relates to the intermediate dynamics. In an ever-expanding 3-fluid model the 
additional degree of freedom is the maximum value attained by \2, which indicates the extent to which 
fluid-2 becomes dominant during the evolution (recall that fluid-1 is dominant at early times (xi — > 1) and 
fluid-3 is dominant at late times {\z ~^ !))• I n an ^^CA-model of the physical universe, the maximum 
value of Xdust is close to unity, which ensures that there is an extended matter-dominated epoch. 

The dynamical systems analysis of the n-fluid FL universes leads to an interesting conclusion regard- 
ing the behaviour of the spatial curvature, as follows. For any e > 0, there is an open set of n-fluid FL 
universes with 7« < | whose spatial curvature scalar satisfies < e throughout the evolution. In 
other words, if j n < |, in particular, if fluid n is a cosmological constant, the total matter-energy content 
"controls" the spatial curvature throughout the evolution. This result is a consequence of the fact that 
generic ever-expanding models satisfy fij. — > as £ — > and since j n < |, fi^ — > as £ — > +00. It 
follows that I I attains a maximum value e which depends on the orbit in question. By restricting to 
orbits that lie within a sufficiently small neighbourhood of the flat FL invariant set Qk = 0, the value e 
can be made arbitrarily small. The above result has been found using different methods by Lake (2005), 
for DC A FL universes. He indicates that the result will also hold for RDCK universes (see also Lake 
(2006)). Lake comments on the significance of this result in connection with the so-called flatness problem 
(Lake (2005)). 

One of the main results of this chapter is the analysis of the RDCK FL universes and the comparisons 
with observations. Discussion of these models dates back to the 1930's. Our analysis is novel in that it 
gives new insight into the dynamics of these models, by representing the evolution as orbits in a compact 
state space. We now give a brief survey of previous work on these models. 

Tolman (1934) gives the Friedmann equation for closed RDCK universes (equation (106.3), p408) 
and states that de Sitter has studied the integration of this equation (de Sitter (1930), (1931), not read- 
ily available). He gives the special solutions of this DE that represent universes that are past or future 
asymptotic to the Einstein static universe, in the form t = f(£) [equations (161.8) and (161.4)]. 
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Coquereaux and Grossman (1982) give a detailed qualitative discussion of the Friedmann equation for 
closed RDC A-uni verses using elliptic and Weierstrass functions. Debrowski and Stelmach (1986) extend 
the analysis to open i?DCA-universes. They give an explicit solution in terms of elementary functions 
when the parameters are suitably restricted. 

The work that is closest to ours is that of Ehlers and Rinder (1989), who give a dynamical systems 
analysis of i?DCA-universes, using Hubble-normalized variables. They give the explicit expressions for 
the density parameters as functions of the scale factor [equation (19)] and also show the existence of 
two conserved quantities [equation (29)]. Their state space representation is inevitably incomplete since 
the Hubble-normalized variables are undefined at the instant of maximum expansion in recollapsing and 
bouncing models. 



Chapter 4 



Perturbations of FL Models 



The high degree of isotropy of the cosmic microwave background (CMB), when combined with the Coper- 
nican Principle, provides strong support for the belief that the large scale structure of the observable 
universe is very well described by the Friedmann-Lemaitre universes, which are exactly isotropic and 
spatially homogeneous. The real universe is, of course, not exactly isotropic and spatially homogeneous, 
since there is complex structure associated with the observed distribution of galaxies, the overall expan- 
sion may not be exactly isotropic, and there may be primordial gravitational waves. But the current belief 
of most cosmologists is that the deviations from an exact FL cosmological model are sufficiently small 
that they can be described by considering linear perturbations of the FL models. 

4.1 Historical Development 

In this section we briefly discuss the history of the theory of linear perturbations of FL cosmologies. There 
are two main approaches, which we shall refer to as the metric approach and the geometrical approach. 

4.1.1 Metric Approach 

The metric approach to perturbations of the FL models was introduced by Lifshitz (1946) and subsequently 
discussed in greater detail by Lifshitz and Khalatnikov (1963). In the metric approach, one distinguishes 
a background spacetime, which is a FL model, and a perturbed spacetime, which represents the physical 
universe. The metric for the perturbed spacetime is written in the form 

9ab = 9(0)ab + $9ab, (4-1) 
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where <?(o) a 6 is the metric of the FL model and 5g a b is a small perturbation, i.e. \6g a f,\ <C \9(o)ab\- The 
fluid 4-velocity and energy density in the perturbed spacetime are likewise written in the form 

u l = u\ 0) +6u\ (4.2) 
M = M(o) + Sfi. (4.3) 

The density contrast is defined by 




(4.4) 



The problem of gauge 

In defining a perturbation, one is effectively establishing a one-to-one correspondence between points 
in the physical spacetime and points in the background spacetime. However, this correspondence is not 
unique: one may make an infinitesimal change in the correspondence keeping the background coordinates 
fixed. Such a transformation is called a gauge transformation. A quantity that is unchanged by an in- 
finitesimal gauge transformation is said to be gauge-invariant. We recommend Ellis and Bruni (1989) for 
a detailed explanation of gauge-invariance. 

The problem of gauge invariance has plagued the study of linear perturbations since the pioneering 
work of Lifshitz (1946), resulting in authors publishing contradictory predictions on the behaviour of 
perturbations of FL cosmologies. The reason is that the metric perturbation Sg a b, the velocity perturbation 
5u l and the density contrast 5 are not gauge-invariant quantities, and hence their time evolution can depend 
on the choice of gauge. 

A result due to Stewart and Walker (1974) gives a useful criteria for identifying gauge-invariant quan- 
tities: 

Theorem 4.1.1 (Stewart- Walker Lemma) Let T/ ) be a tensor field on a background spacetime and let 
T = T( ) + AT be the corresponding tensor on a perturbed spacetime. If Tt \ = then T is gauge- 
invariant. 

In 1980 Bardeen reformulated the metric approach using gauge-invariant variables. He did not make 
use of the Stewart- Walker Lemma in choosing the gauge-invariant variables. Instead, he wrote out the 
transformation laws for the perturbations 5g a b, <5u a and S under an infinitesimal gauge transformation 
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and then, by inspection, formed linear combinations of the perturbations and their derivatives that were 
invariant. 

This approach placed the analysis of perturbations of Friedmann-Lemaitre models on a sounder math- 
ematical foundation and has been used in much subsequent research^ Nevertheless, the metric approach 
still possesses several shortcomings: 

i) First, it is not clear what physical quantities Bardeen's gauge-invariant quantities correspond to. 
Indeed, a gauge-invariant quantity can have different physical interpretations in different gauges. 
The most important ambiguity concerns the density contrast 5, as defined by (14.4I ). Bardeen defines 
two gauge-invariant quantities in terms of 5, denoted by e m and e g . The first equals the density 
contrast 5 when one uses the so-called co-moving gauge, while the second equals S in the so-called 
Newtonian gauge (see Bardeen 1980, p22). 

ii) Second, as emphasized by Hawking (1966), the metric tensor is not a physically significant quantity 
since one cannot measure it directly. 

4.1.2 Geometrical Approach 

The geometrical approach has its origins in a paper by Hawking (1966). Instead of using the components 
of a perturbed metric as basic variables, he proposed using the evolution equations for the kinematic 
quantities, the Weyl curvature and the matter density to study how perturbations evolve. Although this 
paper initiated a significant new approach to studying linear perturbations of FL, the analysis of density 
perturbations was flawed and the paper had little impact. Subsequently, motivated by Hawking's paper, 
Lyth and Mukherjee (1988) used the evolution equations for \x and H to study density perturbations in an 
FL cosmology, although they did not do so in a fully gauge-invariant way. The decisive step was taken 
by Ellis & Bruni (1989), when they proposed using the spatial gradient of the matter density (V /z) as 
the basic variable to describe density perturbations. Unlike the density contrast 5, the spatial gradient is a 
gauge-invariant quantity by the Stewart-Walker lemma, since it is zero in any FL model (see proposition 
I2.1.4I ). The first comprehensive discussion of the geometrical approach was given by Bruni, Dunsby and 



'We refer to Mukhanov (2005) as a recent text that uses a simplified version of the Bardeen approach, restricting considera- 
tions to perturbations of flat FL (see chapter 7). 
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Ellis (1992B 

They also related the Bardeen approach to the geometrical approach, and gave clear physical 
interpretations of his gauge-invariant variables. 

In this chapter we give an introduction to the geometrical approach to linear perturbations of FL 
models. 

4.2 The Linearized Einstein Field Equations 

The geometrical approach is based on the conviction that one should use variables that have a direct 
physical or geometric meaning, and are both coordinate-independent and gauge-invariant. One thus uses 
tensorial quantities that are identically zero in an FL universe, and hence gauge-invariant, on account of 
the Stewart- Walker lemma. 

We shall see that a complete set of tensorial quantities is the following: 

i) The density and Hubble gradients: V a /i, V a H 

ii) The fluid kinematic quantities: a a b, w a 

iii) The Weyl curvature: E^, H a b 

These variables have two additional advantages: 

a) they directly describe the deviation of the physical universe from an idealized FL universe, and 

b) they are, in principle, observable, when normalized to be dimensionless using a power of H. We 
refer in particular to the fundamental paper of Sachs and Wolfe (1967), who analyzed the observa- 
tions of distant galaxies from this point of view, and Maartens et al. (1995b), who used the cosmic 
microwave background observations to bound these quantities. 

In the geometrical approach one thus replaces the EFEs for the metric tensor components by the 
equivalent system of evolution equations and constraints, satisfied by the tensorial quantities i)-iii), as 
given in Appendix lAl 

For simplicity we will assume the matter content consists of a perfect fluid with barotropic equation 
of state p = p(/d) and a cosmological constant A. We will use the standard notation, 

P 2 dp 

w = -, c s = — , (4.5) 
jj, djj, 



The name "geometrical approach" was introduced in this reference. This approach is also referred to as "the gauge-invariant 
and covariant (GIC) approach". 
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where c s is the speed of sound. 
4.2.1 The Linearization Process 

In order to derive the governing equations for linear perturbations of a FL cosmology, we must distinguish 
between quantities which are non-zero in the background (referred to as zero-order quantities) and those 
which are zero in the background (referred to as first-order quantities). For a FL background, the zero- 
order quantities are 

/i, w, (? s and H. (4.6) 
It follows from propositions 12. 1 . 1 H2. 1 .41 that the quantities 

V a /x, V a H, a ab , uj a , E ab , H ab , (4.7) 
are first order. For sake of brevity, we define 

X a = V„/x, (4.8) 

and 

Z a = 3X7 a H. (4.9) 

We note that the acceleration ii a is not an independent quantity, since it can be expressed in terms of 
the density gradient X a via dA.171 ), on using ( 14.51 ) and ( 14.81 ): 

u a = - n j s X a . (4.10) 
(1 + w)n 

The evolution equations for the zero-order quantities are simply the evolution equations for a FL 
cosmology with a single fluid and cosmological constant, given by (13.51 ) and (I3.6I ). with \i replaced by 
fx + A and p by p — A to make the cosmological constant explicit. We also include the evolution equation 
for w, which follows from (14.121) and (14.5b - 



Zero-Order Evolution Equations 






H = 


-H 2 - |(1 + 3w)/j, + |A, 


(4.11) 


A = 




(4.12) 


w = 


-3H(l + w)(c 2 s -w). 


(4.13) 
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Note: The Friedmann equation ( 13.81 ) plays an auxiliary role in that if K ^ 0, it serves to express £(t) in 
terms of the zero-order quantities. We use (13.111 ) to write (13.81 ) as follows: 



The evolution equations for the first order quantities are obtained from the general system of evolution 
equations and constraints in Appendix [A] by a process of linearization. There are three aspects to the 
linearization procedure: 

L\) Drop all products of the first order quantities (14.71 ) and their spatial derivatives in the general non- 
linear evolution equations. 

Li) The resulting equations form a coupled system of linear PDEs for the first order quantities, whose 
coefficients depend on the zero order quantities. For example, there are terms such as HX a , HE a b 
and fia a b- The linearization process involves replacing the zeroth order quantity by its value in the 
background FL model in each of these products. 

L3) The resulting equations also contain the projected spatial covariant derivative operator V a acting 
on first order quantities. The linearization process involves replacing V a , which is defined on the 
physical spacetime, by the corresponding operator on the background FL spacetime. 

4.2.2 The Linearized Evolution and Constraint Equations 

We now derive the governing equations for linearized perturbations by applying the linearization proce- 
dure L\ to the general system of evolution equations and constraint equations in Appendix lAl 

One aspect of the linearization process, namely dealing with the terms V/ um, curl(u a ) and V a (V b Ub) 
involving the spatial derivative of ii a , requires particular attention. First, it follows from ( 14.101 ) and L\ 



H 2 = W + V- W 



(4.14) 



thai 



,1 




(4.15) 



Since p = p(fx), we can use J4.5t to write 



(1 + w)fi 



F(n)V a » = F(n)X, 



where F(fi) is a function of /i whose specific form is unimportant. Then the term F(jj,)X a Xb is dropped. 
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Second, it follows in a similar manner from ( 14.101 ), ( 14.8D and L\ that 

-c 2 

curl(<i a ) = s curl(VaM)- ( 4 -!6) 
(1 + 

Applying (1E.25I ) with / = \i and using (IA.16I) gives 

cur\(u a ) = 6c 2 Hu a . (4.17) 
Finally, it again follows from (14.101 ) and repeated application of L\ that 



V a (V%>) = C ° V a (V». (4.18) 
(1 + w)n 

Then the commutation property (IE. 181 ) with f = (i and the evolution equation (IA. 16b for \i leads to 

— c 2 / 2K\ 
V (V 6 « 6 ) = V 2 - X a - 6c 2 FcurlK). (4.19) 



2 



(l + to)A* V J 

The evolution equations for the first order quantities are now obtained from (1A.9I) (IA. 10b . (IA. 18b . 
dA. 19b . (lA.25b and (IA.26b on applying the linearization procedure L\ and using (14.15b . (14.17b and (14. 19b 
to eliminate terms involving spatial derivatives of u a . 



First-Order Evolution Equations 




X(a) 
Z(a) 


= -4HX a -{l + w)fiZ a , 

-■ -3HZ a - \X a - f' fv 2 + £)x a - 6c 2 F curIK), 
(1 + w)n \ i A ) 


(4.20) 
(4.21) 


<?(ab) 


= 2Ha a b M s V/ a X b \ E ab , 
= (3c 2 -2)^ a , 


(4.22) 
(4.23) 


E {ab) 


= -3HE ab + cur\(H ab ) - ±(1 + w)/j,a ab , 


(4.24) 


H(ab) 


= -3HH ab -cur\(E ab ). 


(4.25) 



Similarly, constraints on the first-order quantities are obtained on applying the linearization procedure 
Li to (CATTi) . (IAT21 . (TAT31 . (COOl and (CA2TT) . 
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First-Order Constraint Equations 




\Z a + curl(a; a ), 



(4.26) 



a 







(4.27) 



V b E ab 



H ab - cur\(a ab ) 



(4.28) 




(4.29) 



(1 + w)(JLU, 



(4.30) 



There are two additional constraints involving the curls of the spatial gradients X a and Z a . Indeed, 
using (14.101) . (14.161 ) and (14.171) leads to 



This equation is in fact a consequence of the main system of evolution equations and constraints and hence 
we list it as an auxiliary equation. There is an analogous equation for curl(Z a ) which we shall not need. 



i) The evolution equations d4.20| )-( |4~2"5T ) and constraints d4.26l )-( |4~3"0l ) form a closed system of linear 
first order PDEs for the first order quantities (14.7I ). The coefficients (jl, w, c 2 s and H are zero-order 
quantities and are evaluated in the background FL model. Note that if K ^ the term K/i 2 in 
( 14.211 ) is given in terms of zero order quantities by the constraint ( 14.141 ). 

ii) It can be shown that the system of evolution equations and constraints is consistent: that is, the time 
derivatives of each constraint equation (14.26l) - (14.30l ) are identically satisfied as a consequence of the 
other equations. As a result, once a set of initial data is chosen to satisfy the constraint equations, 
these equations will hold at all later times (see Bruni et al. (1992), p42). 

iii) Equations (14.20l) - (14.25l) and (I4.26b - (l4.30b form the basis for the geometrical approach to the analysis 
of linear perturbations of FL models with a barotropic perfect fluid and cosmological constant as 
source. They can be extracted from Bruni et al. (1992) by specializing their more general equations 



curlpQ = 6(1 + w)fiHu, 



(4.31) 



Comments: 
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to the case of a barotropic perfect fluid (set q a , -k^ and $ to zero in equations (50)-(54), (57)-(62)). 
Most authors, however, do not make use of the full system of evolution and constraint equations. 

4.3 Harmonic Decomposition 

Before proceeding, we must discuss the fundamental role of spatial harmonics in the geometrical theory. 
The harmonic decomposition is analogous to a Fourier decomposition of a quantity that is defined along 
the background 3-spaces. Hence, harmonics are defined so as to be constant along flow lines (i.e., inde- 
pendent of proper time) and orthogonal to the fluid 4- velocity u. We refer the reader to Appendix |F| for a 
complete discussion of spatial harmonics in the geometrical approach. 

The key idea behind the harmonic decomposition is that a scalar, vector or tracefree symmetric rank 2 
tensor orthogonal to u a can be written as a linear combination of scalar, vector and tensor harmonics with 
purely time-dependent coefficients. In the literature it is tacitly assumed that the harmonics are complete 
in some appropriate space of tensor fields. 

The expansions are as follows: 

T = J2T m Q(°\ (4.32) 

k 

T a = ^T (0)A) Q(°)+^T (1)fe) Q«, (4.33) 

k k 

T ab = E^o^+E^oS+E^oS- (4 - 34) 

k k k 

This decomposition has two consequences: Firstly, it decouples the temporal and spatial dependence of 
each field, since all coefficients are purely functions of time and harmonics are constant along flow lines, 
i.e. Tu t k) = T(d,k)(t) an d Q = for all harmonics Q (see Appendix |F]). Second, in the framework of 
linear perturbation theory, it removes any co-dependence between scalar, vector and tensor harmonics. As 
in Bruni et al. (1992, p51), the summation in d4.32| )-( |4~3"4l ) can be over a discrete set or an integral over a 
continuously varying index. 

We note that the scalar decomposition ( 14.321 ) is analogous to the 3-dimensional Fourier transform 
for scalar fields. The vector decomposition (14.331) is similarly analogous to the Fourier transform, ex- 
cept applied to the curl-free (scalar) component and the divergence-free (vector) component of the vector 
field. The tensor decomposition (14.341 ) separates the curl-free divergent (scalar) component, the remaining 
divergent (vector) component and the divergence-free (tensor) component of the vector field. 
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Using this decomposition, it is natural to consider the kinematic quantities and Weyl tensor compo- 
nents of the real universe to be the fundamental description of the linear perturbations (since they are each 
gauge-invariant, each of these components can be expanded using the harmonic basis functions). 

The decoupling of the perturbation types is extremely important in the linearized theory since it allows 
us to consider scalar, vector and tensor perturbations independently. For each perturbation type, we will 
proceed as follows: First, we first define a basic variable. Second, we use the linearized evolution equa- 
tions (14. 1 1 b - d4~30T > to derive a second order differential equation which describes the evolution of the basic 
variable, which we shall call the governing DE. Finally, we use the linearized evolution and constraint 
equations to express all kinematic quantities and Weyl tensor components in terms of the basic variable, 
hence providing a complete solution to the linearized equations for each type of perturbation. 



4.4 Scalar Perturbations 

We first focus on scalar perturbations, which are defined as solutions of the linearized equations (14. 11b - 
( 14. 301 ) that can be expanded in terms of the scalar harmonics Q(°\ and Q^, defined by equations 
(1F3I ). (IF 51 ) and (1F6I ). Following Ellis et al. (1989), we choose the basic quantity for scalar perturbations 
to be the fractional density gradient V a , defined as 

V a = (4.35) 
A* 

or equivalently, using (14.81) . 

V a = ~X a . (4.36) 
We also rescale Z a according to Z a = £Z a , i.e. 

Z a = 3£V a H. (4.37) 



4.4.1 The Governing DE 

In this section we derive the governing DE for scalar perturbations, following Ellis et al. (1990) (also see 
Wainwright and Ellis (1997, p290-3)). Since we are considering scalar harmonics, we can write 

T> a = J2' D (k)Qa' k) - ( 4 - 38 ) 
k 
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Since (l4~38l) and dEE) imply curl(£> a ) = 0, and hence curl(X a ) = 0, (|43TT> gives 

= 0. (4.39) 

It follows that equations (14.201 ) and (14.211 ) decouple from the remaining linearized evolution equations. 
Expressing these equations in terms of the rescaled variables V a and Z a leads to 

V {a) = 3wHV a -(l + w)nZ a , (4.40) 
Z {a) = -2HZ a -l,V a - J ^(v* + §)v a . (4.41) 

These equations can be combined to give a second order evolution equation for V a . On differentiating 
(14.401 ) and using (14.11I )- (I4.13I ). (14.141 ) and (14.411) . we obtain 

V {a) + AHV {a) + BH 2 V a - c\ (V - ??pj V a = 0, (4.42) 

where 

A = 2-3w-3(w-c 2 s ), (4.43) 

and 

B = — 1(1 — w)(l + 3u»)0 - 6wCl A + 3(w - eg) (-^ - . (4.44) 

We now expand each term in terms of scalar harmonics. The derivatives of V a along the fundamental 
congruence are obtained by differentiating (14.381 ) and noting Q^a' k ^ = (see (IF. 101 )). as follows: 

A«> = EA*M 0,fc) > %) = EA*M ' fc) - ( 4 - 45 ) 

k k 

In order to remove the Laplacian term from the evolution equation (14.42b . we use the identity (IF. 131 ). which 
can be written as 

V 2 -^)Qi 0) = -^i 0) - (4-46) 
On applying (14.451 ) and (14.461 ) to (14.42b . we obtain 

V {k) + AHV {k) + (b + ^) H 2 V (k) = 0, (4.47) 
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where A and B are given by ( 14.431 ) and (14.441 ). This is the governing equation for scalar perturbations, in 
terms of clock time t. This equation is equivalent to (73) in Bruni et al. (1992), upon specializing to the 
case of a perfect fluid, and (214) in Ellis and van Elst (1998), in the case of A = 0. 

4.4.2 Kinematic Quantities and Weyl Curvature 

We now write the kinematic quantities and Weyl curvature in terms of the density gradient coefficients 
and the scalar harmonics Q(°\ Since we are considering scalar harmonics, Z a and a ab can be 
expanded as 

Z a = ]T Z {k) Q^ k \ a ab = £ a {k) Q^\ (4.48) 

k k 

We can use ( 14.401) to express Zi k \ in terms of T>i k \ and and since u> a = it follows from ( 14.261 ) and 
(IRT6l) that 

Z {k) = (4-49) 
Next, it follows from (14.481 ) and (lE8l) that curl(cr afe ) = 0. Hence, using (14.281 ) we have that 

H ab = 0. (4.50) 

Expanding E ab in terms of scalar harmonics via 

E ab = Y,E {k) Q ( ^\ (4.51) 

k 

and using ( IF. 161 ) leads to 

£V b E ab = Y^^(k 2 - 3K)E (k) Q^ k \ (4.52) 

k 

The coefficients Ei k \ can then be expressed in terms of X?^) via the constraint (14.29b - We thus obtain 
the following expressions for the harmonic coefficients of the non-zero Hubble-normalized first-order 
quantities: 



4 Bruni et al. (1992) make use of the scalar quantity A= \7 a T> a instead of T>^, denned in this manner so as to restrict 
considerations to scalar perturbations without it being necessary to resort to a harmonic expansion. This choice of scalar quantity 
automatically removes all vector harmonic components from D a since vector harmonics are divergence-free by definition dF.22b . 
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(1 + w) 
k 



k 2 



k 



2(k 2 - 3K) 



(4.53) 
(4.54) 
(4.55) 



At this stage, we have satisfied all of the linearized evolution equations and constraint equations except 
for (14.221 ) and (14.241) . We can now verify that D a , a ab and E ab , as given by (14.38b . (14.48b . (14.511 ) and 
d4.53ft - d4.55l) satisfy ( |4.22| > and ( |4.24| > (via Maple). We have thus established the following result: 



Proposition 4.4.1 For any solution XV^n of the governing DE A4.47\) , the first order quantities given by 
H4.38\) , ( 14.481 ), ( |4.5iD f |4.5J| )-( |4.55D satisfy the complete set of linearized evolution and constraint 
equations i4.20h$4JPb . 



4.5 Vector Perturbations 

In this section we focus on vector perturbations, defined as solutions of the linearized equations ( |4.201 >- 
(14.30b that can be expanded in terms of the vector harmonics Qa and Q Jl. Vector perturbations are 
closely associated with vorticity, and so we choose the basic quantity for vector perturbations to be the 
co-moving vorticity vector, defined by 

W a = lu a . (4.56) 

4.5.1 The Governing DE 

In this section, we derive the governing DE for vector perturbations. On differentiating (14.56b along the 
fundamental congruence and using d4.23b , we obtain 

W {a) - (3 C 2 - l)HW a = 0. (4.57) 

We expand (14.57b in terms of the associated vector harmonics P^' k ^ as 

Wa = Y, W (k) P a 1,k) - ( 4 - 58 ) 

k 
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On substituting (14.581 ) into ( 14.571) , we obtain 

W (k) - (3c 2 s - l)HW [k) = 0. (4.59) 
This equation is the governing DE for vector perturbations, in terms of clock time t. 

4.5.2 Kinematic Quantities and Weyl Curvature 

We now write the kinematic quantities and Weyl curvature in terms of the co-moving vorticity coefficients 
W(fc) and vector harmonics By the definition of vector perturbations, all quantities can be expanded 
in terms of vector harmonics and associated vector harmonics, as follows: 

k k k 

E ab = £ E (k) Q% k \ H ab = H W P£*\ (4.60) 

k k 

From (14.3 IK (14.361 ) and (IF. 23 1) , we can express T>(u\ in terms of Wm, and hence obtain Zi k \ in terms 
of W(fc) from (14.20b . (14.361 ) and (14.371) . Further, using (IF34I) . the divergences of a ab , E ab and H ab can be 
written as follows: 



W b *ab = £^(^-2if)a (fc) Q , (4.61) 

k 

lV b E ab = ^^(k 2 -2K)E (k) Q a , (4.62) 

k 

£V b H ab = J2jr(k 2 -2K)H {k) Q a . (4.63) 

k 

We can now use ( 14.261 ), ( 14.291 ) and ( 14.301 ) to determine the coefficients era.), Eqa and Hn^ i n terms of 
W(fc). We thus obtain the following expressions for the harmonic coefficients of the first-order quantities: 

V (k) = 6(1 + w)H£W (k) , (4.64) 

Z {k) = -U- 1 [2K-{l + w)^ 2 ] W (k) , (4.65) 
2k 

a (k) = h 2 o rr r 1 [( fc2 - 2if) + 2(1 + w)fx l 2 ] W (k) , (4.66) 



2k 


k 2 


- 2K 




4k 


k 2 


- 2K 




2k 


k 2 


- 2K 



E {k ) = k2 _ 2K (l + w)»H£W {k) , (4.67) 
2k 

H (k) = + (4 - 68) 
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At this stage, we have satisfied all of the linearized evolution equations and constraint equations except 
for (14.221 ). (14.24b . (14.25b and (14.28b . We can now verify that V a , Z a , a ab , u a , E ab and H ab , as given by 
d4~60l ) and d434l - d438l ) satisfy d432l . d4~24l ), ( 14251 ) and d428l (via Maple). We have thus established the 
following result: 

Proposition 4.5.1 For any solution W^) of the governing DE A4.59\) . the first order quantities given by 
H4. 601 ) and H4.64\) - 0.68\) satisfy the complete set of linearized evolution and constraint equations H4. 20\ - 

EM- 

4.6 Tensor Perturbations 

In this section we focus on tensor perturbations, defined as solutions of the linearized equations (14. lib - 
(14.30b that can be expanded in terms of the tensor harmonics Q ab . Tensor perturbations are interpreted as 
gravitational waves, described by a coupling of the electric and magnetic Weyl curvature. 

As regards to the choice of basic variable, we shall see that for tensor perturbations the only non-zero 
first order quantities are a ab , E ab and H ab . On the grounds of mathematical simplicity we choose the 
comoving shear tensor as the basic variable, denoted X ab and defined by 

Xab = iOab- (4-69) 

4.6.1 The Governing DE 

In this section, we derive the governing DE for tensor perturbations]^ It is a consequence of applying 
tensor harmonics that rank two tensors have zero spatial divergence, i.e. 

V h a ab = ^ V b E ab = 0, V b H ab = 0. (4.70) 

It then follows from the constraints ( 14.26b , ( 14.29b and ( 14.30b that all vector quantities are zero, i.e. 

V a = 0, Z a = 0, and co a = 0. (4.71) 

Using (14.71b , the shear evolution equation ( 14.22b can be rewritten in terms of X ab as 

x (ab) = ~HX ab - £E ab . (4.72) 



5 A similar derivation was given by Dunsby et al. (1997). 
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We now use ( 14.241 ) to calculate {£E ab )\, noting that £ = H£. We use the constraint ( 14.281 ) to eliminate 
H ab and then apply (1E.31I ) to the "curl curl X ab " term. The result is 

(iE ab Y ± = -2H(£E ab ) + + w)fi + jpj X ab - V 2 X ab . (4.73) 

We can now combine these equations to give a second order DE for X ab (differentiate (14.721) using (14.111) 
and (14.73b . and then use (14.721 ) to eliminate E ab from the resulting equation). After expressing the coeffi- 
cients of X ab in terms of the density parameters Q and Q\ (see !3.3 lk we obtair 



{ab) 



-(3u> + l)ft + 2ft A + 



2K 
HM 



H X ab — V X ab = 0. 



We now expand X ab in terms of tensor harmonics according to 

Xab = S ^ J X (k)Qg' h 



' ab 



(4.74) 



(4.75) 



Substituting ( 14.751 ) into (14.741 ) and using dF.381 ) yields 



(*) 



k 2 + 2K 
H 2 £ 2 



H 2 X^ = 0. 



(4.76) 



This is the governing equation for tensor perturbations, in terms of clock time t. 



4.6.2 Kinematic Quantities and Weyl Curvature 

We now write the kinematic quantities and Weyl curvature in terms of the dimensionless shear coefficients 
X(u\, tensor harmonics Q^jj and associated tensor harmonics pj£ '. Since ( 14.701 ) confirms that a ab , E ab 
and H ab can be expanded in terms of tensor harmonics, we can write 

O-ab = ^2 a (k)Qab k) > E ah = ^ E {k)Qa b ' k \ H ab = ^ H (k) P ab' k ^ ( 4 " 77 ) 
k k k 



6 It should be noted that the evolution equation (22) given in Dunsby et al. (1997, pl219) contains two errors. The original 
equation should instead be written as 

Aa ab + fe<T a{ , + [|e 2 - ± m (9t - io)] (j ab = o, 

where Aa a t, = a± a i, — V 2 u a f,, so as to be consistent with the remainder of that paper. This equation is then equivalent to l |4.74t 
for A = on recalling that X a b = £<Jab and Q = 2>H. 
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where |3| 

P ab = £cur\(Q ab ). (4.78) 

Then using (14.691 ), (14.221 ) and (14.281 ) we can write a ab , E ab and H ab in terms of Xha. We thus obtain the 
following expressions for the harmonic coefficients of the non-zero first-order quantities: 



V(k) = r x X {h) , (4.79) 
E {k) = -r l (x {k) + HX {k) ), (4.80) 
H [k) = r 2 X {k) . (4.81) 

At this stage, we have satisfied all of the linearized evolution equations and constraint equations except 
for (14.241 ) and (I4.25I ). We can now verify that a ab , E ab and H ab , as given by ( 14.771 ) satisfy ( 14.241 ) and (14.251 ) 
(via Maple). We have thus established the following result: 

Proposition 4.6.1 For any solution Xha of the governing DE H4. 76\ , the first order quantities given by 
A4. 77\ and A4. 79D - rf?T571) satisfy the complete set of linearized evolution and constraint equations \4.20\ - 

a. 



4.7 Dynamics of the Linear Perturbations: General Features 

In this section we discuss the dynamics of the linear perturbations. We first discuss the importance of 
conformal time in analyzing perturbations of FL and rewrite the governing DEs for each perturbation type 
in terms of this time variable. Finally, we examine the dependence of the first-order quantities on the 
solutions to the governing DEs. 

4.7.1 Conformal Time and the Particle Horizon 

It is a consequence of propositions 14.4.11 14.5.11 and 14.6.11 that the spatial gradients, kinematic quantities 
and Weyl tensor components can be determined directly from the set of solutions of the governing DEs. 
As a consequence, we turn our attention to the behaviour of solutions to these equations. 



7 See definition in section |R3l 
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We make the choice of conformal time 77 as the independent quantity, defined in terms of the length 
scale i and clock time t by 

§4 



Under this choice of time variable, the Roberson-Walker line element d 1 - 31 ) can be written in the form 

ds 2 = £ 2 (7])(-dri 2 + g a/3 dx a dxP), (4.83) 

where g a p is a 3-metric of constant curvature. 

Conformal time is directly linked to the notion of the particle horizon, defined as the boundary of the 
part of the universe that is visible to us. In terms of conformal time, the distance to the particle horizon 
(see figure |4~T| ) is given by 

d H = ^(77)77. (4.84) 



A(r?) = ( ^ ) £(v), (4-85) 



(4.86) 



Since the wavelength A for a perturbation of mode k at some fixed time 77 is given by 

/27i 

then from (14.841 ) we obtain the key relation 

A _ 2vr 
dn krj 

Hence, if krj <C 1 the wavelength of the perturbation mode is large compared to the distance to the particle 
horizon (since A 3> du )■ In this situation it is customary to say that the perturbation is outside the horizon. 
Conversely, if krj 3> 1 the wavelength of the perturbation mode is small compared to the distance to the 
horizon (since A <C du)- Similarly, in this situation it is customary to say that the perturbation is inside 
the horizon. 

The range of the conformal time variable varies depending on the characteristics of the FL model, as 
shown in Appendix [D] A summary of the behaviour of conformal time and clock time for each of the 
generic classes of FL model with an initial singularity is given below: 



FL Model Range of r\ Range of t 

Ever-expanding models, A = < 77 < +00 < t < +00 

Ever-expanding models, A^O 0<r/<r/j 0<t<+oo 

Recollapsing models, A arbitrary < 77 < i]f <t < tt 



4.7. Dynamics of the Linear Perturbations: General Features 



93 



i 


I 

dn(r[2) 






dn(r[ i) 

















11 = 111 



FIGURE 4.1: The relation between the distance to the particle horizon, denoted dnirj) an d the conformal 
time parameter r\. 

4.7.2 Governing DEs Using Conformal Time 

We now present the governing DEs for scalar, vector and tensor perturbations in terms of conformal time 
and specialized to a 7-law equation of state, i.e. p = (7 — l)p,. It follows from (14.51) that 

w = c? s = 7 - 1. (4.87) 

The three governing DEs, given by (14.47b . ( 14.59b and (14.761) in terms of clock time t, assume the following 
form in terms of conformal time. The change of variable can be performed simply by using these identities: 

t± = ±, e—^ — -Hl- i (4.88) 
dt dry' dt 2 dr] 2 dr]" 1 

which follow from (14.821 ). 
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Governing DEs for Perturbations using Conformal Time 








Scalar 








V( k) + (4 - 3 7 )H£V{ k) + (-|(2 - 7 )(3 7 - 2)Q - 6(7 - 


1)«A + 


r 2 k 2 \ 

j%) *&> w = 0. 


(4.89) 


Vector 








W (k) -(&y-4)H£W {k) 


= 0. 




(4.90) 


Tensor 

X[ k) + 2H£X' {k) + (-(3 7 - 2)fl + 2n A + - 


2 + 2K N 
H 2 £ 2 , 


) H 2 £ 2 X {k) = 0. 


(4.91) 



In the case of vector perturbations, the governing DE ( 14.901 ) can be integrated directly, on applying^] 



H = j r (4.92) 



W (fc ) = C(\t)^~\ (4.93) 



We then obtain 

' A/ (fc) 

In the case of scalar and tensor perturbations, the second-order governing DEs can be transformed into 
normal form by making a change of dependent variable. In both cases the resulting DE§| has the following 
form, in terms of a constant frequency u and a potential U (rj): 

Yfc + [co 2 - U( V )] Y {k) = 0. (4.94) 

The wave number k enters only into the frequency, which also depends on the curvature index K, and, in 
the case of scalar perturbations, on the equation of state. The potential U(rj) depends on the background 
FL model (H, £, O and Oa). The specific DEs, given in the table below, are obtained by making the 
indicated change of variable in the governing DEs ( 14. 891 ) and ( I4.91I ). In order to simplify the subsequent 
analysis we use the notation 

2 2 
[5=- -, or, equivalent^, (3 = — , (4.95) 

37 — 2 3w + 1 



8 This expression is obtained from l |3.9t on making the change of variable l !4.88t . 

9 To the best of the author's knowledge, this similarity between tensor and scalar perturbations has not been noted before. 
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as inBardeen (1980). 



Normal Form of Governing DEs for Scalar and Tensor Perturbations 




Scalar 




1 l " 

Change of variable: VnA = (\l)P £Vfc), 

Governing DE: X>'( k) + [uj 2 - U(r])} V {k) = 0, 


(4.96) 
(4.97) 


with 




uj 2 = c 2 k 2 + ^((3-l) 2 K, 

u(v) = ^(i + f3)(pn + n A )H 2 £ 2 . 


(4.98) 
(4.99) 


Tensor 




Change of variable: X/q = (A^) -1 Xn^, 

Governing DE: X'fa + [uj 2 - U(rjj\ X {k) = 0. 


(4.100) 
(4.101) 


with 




uj 2 = k 2 + 3K, 
U(rj) = ±(l + /?)Otf 2 £ 2 . 


(4.102) 
(4.103) 



Note: Here A is the constant conformal parameter introduced in section l3T2~Tl having dimension (length) 

Remark: Writing the governing DEs in normal form has several advantages. Firstly, it gives insight into 
the qualitative behaviour of the solutions. In particular, if lo 2 3> U(rj) the solutions will be approximated 
by oscillatory solutions of frequency u, while if uj 2 <C U (??), the solutions will be approximated by power 
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law solutions. Secondly, if the background FL model has A = 0, the general solution can be written in 
terms of Bessel or Legendre functions, as we now show. 

4.7.3 General Solutions for Scalar and Tensor Perturbations with A = 

In the case A = 0, the potential U (ry) in the governing DEs is given by 

U(rj) = ^-{l + p)nH 2 e 2 (4.104) 

P 

for both scalar and tensor perturbations (see equations ( 14.991 ) and (14.1031) ). The background solution is 
(PX2061 ): 

£(rj) = \- 1 [S K (T)f , T = ^, (4.105) 

where 

S K (T) = (T, smh T, sin T) , (4. 106) 
for K = (0, -1, +1). On noting that Q,H 2 l 2 = \ixi 2 , it follows from (ET207T) that 

QH 2 t 2 = (\t)~k (4.107) 

Thus, using (14.1051 ) we obtain 

flH 2 £ 2 = [S K {T)]~ 2 . (4.108) 

Observe now that the governing DEs in normal form (equations ( 14.971 ) and (14.1011 )) coincide with the 
Bessel DE or the Legendre DE, as given in Appendix [Hi ( see equations (1H.1I) - (1H.6I) ). with the parameters 
determined as follows: 



Curvature DE Type Parameters 

K = Bessel DE fj, = /3 + \ a = uj 

K = — 1 Legendre (Toroidal) DE /x = /3 + i v = — i + ifito 

K = +1 Legendre (Conical) DE /j, = j3 + \ v = -\ + (3uj 



The value of to, as given by ( 14.981 ) and (14.1021 ) distinguishes between scalar and tensor perturbations. The 
solutions of the three DEs are given by (1H.2I ). (IH.4I ) and (1H.6I ). in terms of Bessel and Legendre functions. 
The resulting expressions for T>^ and obtained from ( 14.961 ), (14.1001 ) and ( 14.1051 ) are listed in the 
following table 



4.7. Dynamics of the Linear Perturbations: General Features 



97 



General Solution of the Governing DEs for Scalar and Tensor Perturbations When A = 



Flat FL (K = 0) 



Scalar 
Tensor 



C + J p+ i(ari) + C-Y g+ i(ari) 



where 



a = (c s k, k) 
for scalar and tensor perturbations, respectively. 



Open FL (K = -1) 

Scalar V^(ji) = (sinh3 



Tensor (rj) = ( sinh 4 



C+Q$ ( cosh 3 ) + C-Pif fcosh 3 



where 



(4.109) 



(4.110) 



(4.111) 



(4.112) 



and the frequency u is given by ( 14.981 ) and (14.1021 ) respectively. 
Closed FL (K = +1) 

Replace hyperbolic functions by trigonometric functions in the solutions (14.1111 ) for open FL, and 

use 

v + \ = (5uj. (4.113) 



Remark: The physically important cases, namely when the matter content is dust (i.e. 7 = 1, /3 = 2) 
or radiation (i.e. 7 = |, f3 = 1) require comment. For dust, the scalar solution is not valid since 
= 7 — 1 = 0. We will derive the solution in this case in section 14.81 The tensor solution is valid, 
however, and the relevant Bessel and Legendre functions are in fact elementary and are given in section 
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14.8.41 For radiation, both the scalar and tensor solutions are elementary. We will give them in section [ 

One can use asymptotic approximations of the Bessel and Legendre functions (Appendix|H]) to obtain 
approximate solutions in restricted regimes. 

4.7.4 Functional Dependence of First-Order Quantities 

We are interested in describing to what extent a perturbed FL cosmology deviates from an exact FL 
cosmology. 

It is known that cosmological observations place bounds on the first order quantities, specifically 
on dimensionless quantities formed by dividing the first order quantities by an appropriate power of the 
Hubble scalar. The specific quantities are 

Kinematic quantities: -£jr,-£, (4.114) 

Spatial gradients: * , (4.115) 

Weyl tensor: §^§£- (4-116) 

The first analysis of this type was due to Sachs and Wolfe (1967) who showed that observations 
of distant galaxies could, in principle, place bounds on these quantities. More recently, Maartens et al 
(1995a,b) have shown that the high isotropy of the CMBR (^Z ~ 1CT 5 ) places strong bounds on these 
quantities!'"! 

It is natural to use these Hubble-normalized quantities as a measure of how close the physical model 
is to an idealized FL model, and we will thus calculate their time dependence for the three perturbation 
types. A growing mode is an indication of an instability in the background FL model. 

We now briefly mention the connection between the geometrical approach and Bardeen's gauge- 
invariant metric approach to cosmological perturbations (Bardeen (1980)). In Bardeen's approach, the 
basic variables are gauge-invariant linear combinations of the metric perturbations. However, it turns out 

(2) 

that Bardeen's metric potentials <&h, * and Hf have the same time dependence as certain expressions 



10 The analysis is complicated and subject to assumptions being imposed on the time derivatives of the multipoles of the 
CMBR. 
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formed from the basic variables in the geometrical approach. It can be shown thaiU 

$ H - nH 2 £ 2 V {k) (scalar), (4.117) 

* - nH 2 £ 2 W {k) (vector), (4.118) 

H^ 2) - 2H£X {k) + X' {k) (tensor). (4.119) 

We refer the reader to Bruni et al. (1992, p45-9) for an overview of the relation between the geometrical 
framework and Bardeen's metric approach. 

Using the results in sections 14.4.21 14.5.21 and 14.6.21 we can obtain the dependence of the Hubble- 
normalized quantities (14.1 14| )- (|4.1 161 ) on the basic variables ~^(k) and Xr k \ in terms of conformal 
time. In order to enable comparisons to be made, we also include Bardeen's metric potentials (14.1171 )- 
( 14.1191 ). These results are provided in table 14.11 



"This result relies on comparing the governing DEs for V^, W(k) and X^ with the governing DEs defined in Bardeen's 
approach (Bardeen equation (4.9), (4.13) and (4.14)) and noting that the paired quantities (X>(fc), e m ), (W(fc), V c ) and (X^, 
satisfy the same evolutions equation. Equations J4.117t -( f4.118b then follow from Bardeen (4.3) and (4.12). 
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Functional dependence of first-order quantities. 



Quantity 

&ab 

H 
H 

H 2 

Egb 

m 

Hab 
H 2 



Scalar Perturbations 



-3( 7 - l)P( fe ) + 




(k) 



Hi 



V 



(fc) 



Hi 



1 

H£ 



3( 7 - l)V {k) - 

nv (k) 





V, 



(fc) 



Hi 



Vector Perturbations 



k 2 -2K 
H 2 i 2 



+ 67^ 



H£W, 



(k) 



(k) 



Hi 



2K 



H 2 i 2 



(fc) 



-3 7 ft 



(fc) 



(fc) 



(fc) 



Tensor Perturbations 



(fc) 








(fc) + X [k) 



Hi H 2 i 2 



(fc) 



H 2 i 2 



if. 



(2) 



tlH 2 i 2 V 



(fc) 



nH 2 fw, 



(fc) 



2H£X^ + Af/j 



L (fc) 



ii 3 



[| 7n _ 3(7 _ i)] 



(fc) 



V' 



+ 



(fc) 



iW ii 2 ^ 2 



K 







H 2 i 2 H 3 i 3 



Table 4.1: Expressions for the first order quantities in terms of the basic variable for each perturbation 
type. 



4.8. Perturbations of FL Models with Pressure-Free Matter 



101 



4.8 Perturbations of FL Models with Pressure-Free Matter 

For the case where the background FL model contains dust and possibly a cosmological constant, the 
analysis of scalar perturbations simplifies significantly since one can avoid performing the harmonic de- 
composition. The evolution equation (14.421 ) for the fractional density gradient V a simplifies to 

V {a) + 2HV {a) - lQH 2 V a = 0. (4.120) 

Since (14.101 ) implies u a = for a dust background, a short calculation using (11.51 ) and (11.181 ) reveals 
h a b = 0, and hence that 

V (a) =V a , and V {a) =V a . (4.121) 

Namely, we use the fact that the fluid 4-velocity is a geodesic in the case of dust. Thus, we can drop the 
spatial projections from ( 14.1201 ) and so obtain 

V a + 2H V a - \SlH 2 V a = 0. (4. 122) 

Although it follows that V a does not depend on the perturbation mode k, a harmonic decomposition is 
required in order to obtain expressions for the shear and electric Weyl tensor (see, for example, (14.541 ) and 
(14351) ). 



4.8.1 The General Solution in Integral Form for Scalar Perturbations 

A second simplification is that the evolution equation (14.1221 ) has a first integral, which we derive as 
follows, using the spatial gradients V a , Z a and 1Z a . 

In the case of dust, the evolution equations, given in (14.40I) . (I4.41I ) and (IB .13b simplify significantly, 
and can be written as 

Va = -Z a , (4.123) 

(£ 2 Z a y = -^£ 2 V a , (4.124) 

(£ 2 TZ a y = —AKZ a . (4.125) 

Note that these DEs are not independent, since the 3 spatial gradients are linearly dependent. This result 
can be seen upon taking the spatial gradient of (11.431) and linearizing so as to obtain 



Tl a = —AHZ a + 2/jlD, 



(4.126) 
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Observe that the quantity C a , defined 



C a = fll a - 4KV a . (4.127) 

satisfies 

^-C a = 0. (4.128) 
at 

We can rewrite (14.1271) in terms of V a t> a using (14. 123b . (14.1261) . (13.311 ) and (13.35b . so as to obtain 

c a = u\Hv a + (§n + n k )H 2 v a ). (4.129) 

This equation is in fact a first integral of the evolution equation (14. 122 1 ), since (14.1221 ) is a consequence of 
(141231) and d4T24l) . 

On noting that, in the case of dust, (14.1 II) can be rewritten as 

h = -(^n + n k )H 2 , (4.130) 

we can further simplify ( 14.1291 ) to read 



d fV n \ C a 



dt\H) AH 2 £ 2 ' 

The solution of (14.1311) is thus 



(4.131) 



V a (t) = CW H J* JL + cH|, (4.132) 
for time-independent vectors d ' and Cd , with Ca = C a . For simplicity, we define 

and so write ( 14.1321 ) as 

V a (t) = C«X>«(t) + C^V^\t). (4.134) 
In terms of conformal time r/, the growing mode X)( + ) (ry), as given by (14.1331) . reads 

V^( V ) = H r^L. (4.135) 



12 Also see Bruni et al. (1992), eq (38). 
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Equation ( 14.1351) can often be used to gain insight into the asymptotic behaviour of dust perturbations. 

If H —> or H — > oo as rj — > rf (where rf = ±oo is allowed), we can rewrite X)( + ) as 

(ft) drj\ 

#> ~ (4.136) 
On applying L'Hopital's rule, in the limit as rj —>■ rj*, we obtain 

v{+) ~ ~Wi' (4 - 137) 

Then (14.1301) gives an expression for the asymptotic behaviour of the growing mode in terms of zero-order 
quantities: 

VM{r,) ~ l a Jp-K - aS "^" - - (4 " 138) 
Although the general solution (14.1331 ) cannot always be obtained in closed form, one can analyze the 
asymptotic behaviour of these solutions. We are primarily interested in three asymptotic epochs, namely 
flat FL (77 w 0), Milne (77 S> 1, w 1) and de Sitter (rj « rjf, S7a ~ 1)- The general behaviour of 
perturbations within these epochs shall be discussed in detail in section l4~T0l 

4.8.2 Scalar Perturbations with A = 

We now use the general solution (14.1351 ) to obtain the solution for a dust background with A = 0. The 
length scale (14. 105b . when specialized to the case of dust, is given by 

\£=(S K (v/Z)) 2 , (4-139) 
and the Hubble parameter by H = £'/£ 2 and (14.1391 ). 

where Ck{T) is defined analogous to Sk(T) to be 

C K {T) = (1, cosh T, cos T), (4.141) 

for K = (0,-1, +1). 

We now consider the case of K ^ separately. In this case, the growing mode can be determined 
from (14.1351 ). ( 14.1391 ) and (14.1401 ). on integrating 
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This integral is elementary and can be expressed in terms of H and £ after a short calculation: 

X>(+) = l-{2-lnH£) + H 2 e. 

At 



It follows from (14.1071) with (3 = 2 that 



£lH 2 i 2 



which allows us to rewrite the Friedmann equation (13.361) as 



(4.143) 



(4.144) 



H 



2„2 



-K + 



We can now eliminate the H 2 £ 2 term in ( 14.1431 ), obtaining 
The solution of (14.1221) for K / is then 



V =C (+) 



(1 - ±77^) - K 



+ 



where i7 and I are given by ( 14.1391 ) and ( 14. 140b . 

In the case of K = 0, the integral (14.1421) evaluates to 



£>(+) = 

^ 10'' ' 



(4.145) 



(4.146) 



(4.147) 



(4.148) 



where the coefficient is chosen so as to be asymptotic to the growing mode solution for K ^ when 
7] —>■ 0. The solution for a flat FL background is then 



V a = C^(^ V 2 )+Ci-\8r,- 3 ). 



(4.149) 



We depict the growing and decaying mode solutions in figure [ 

We note that the solution (14.1471 ) can also be derived from the governing DE (14. 122b . upon applying 
the general solutions dH.21 ), (IH.34I) and dH.351) . 
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FIGURE 4.2: A depiction of the qualitative behaviour of the growing and decaying mode solutions of 
&4.122\) when A = 0. From left to right, the plots depict negative (K = —1), Hat (K = 0) and positive 
(K = +1) curvature. 



4.8.3 Scalar Perturbations with A > 0, K = 



In the case A > with flat spatial geometry (K = 0), the general solution ( 14.1331 ) can be given in terms of 
toroidal Legendre functions with clock time as the time variable, as follows. We transform the governing 
DE (14.1221 ) to normal form by means of the change of variable 



(\l)~ l V a . 



The result is 



T>a ~ H 2 V a = 0, 



Vn 



(4.150) 



(4.151) 



The length scale of the background solution is given by (13.225I ). We calculate H using H = £/£. The 
results are 

t = A-^sinhT) 2 / 3 , # = AcothT, (4.152) 



where 

T = |At, and A = J^. (4.153) 

We substitute (14.1521) into (14.1511) and use T as the independent variable. After applying a standard 
identity, we obtain 

d 2 



dT 2 U " 9 



1 + 



1 

sinh 2 T 



V a = 0, 



sinh ( iV3At 



-2/3 



(4.154) 
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This DE is the toroidal Legendre DE dH.31 ) with \i = ~ and u = i whose solution is given by (IH.4b . Thus, 
we have 

2? a (t) = (sinhT)-i rc( + )Q^(coshT) + C^P^ (cosh T) . 



'1/6 V ; ^ - 1 1/6 

The second term of this solution is in fact an elementary function due to the identity 



(4.155) 







i 



(smhT)-tP i 5 / / 6 6 (coshr) = Acoth(T), (4.156) 
for A constant. It follows from (14.1521 ) that (14.1551 ) can be rewritten as 

V a (t) = CW(sinhT)-HQ^(coshT) + C^j, (4.157) 

where T is given by (14.1531) . Observe that the second term corresponds to the second term in the general 
solution ( I4.132I ). However, the growing mode in ( 14.1571 ) does not correspond to the growing mode given 
in the general solution (14.1321) . since it diverges as t — > + (see dH.391 )). Thus the growing mode £>(+)(» 
in (14.1321 ) must be a linear combination of the two solutions in (14.1571) : 

v(+) = H f Q Jf% = ai(sinhT)-sQ^(coshr)+a 2 y . (4.158) 
The growing mode defined in this way will satisfy 

lim£> (+) =0. (4.159) 
We depict the growing and decaying mode solutions in figure 03] 

4.8.4 Tensor Perturbations for A = 

The solution for (rj) is obtained in terms of Bessel and Legendre functions by choosing (3 = 2 in the 
general solution ( 14.1091) , (14.1111) and (14- 113b - Since the index has the value \i = |, these solutions are 
elementary functions, given by (IH.20b - (IH.25b . Upon comparing the functional dependence in each case, 
we obtain the following unified form of the solution: 

*(*)07) = y £ [c+* ( t)M + , (4-160) 



13 This result is obtained from Abramowitz and Stegun (1970), equations (8.5.1) and (8.6.16). 
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Figure 4.3: A depiction of the qualitative behaviour of the growing and decaying modes in an FL 
cosmology for density perturbations with dust and cosmological constant. 

where 

X+Jn) = {An 2 - K - 3H 2 £ 2 ) sm{nn) + 6nH £ cos(nr/) , (4.161) 
XZJrj) = {An 2 - K - 3H 2 £ 2 ) cos{nn) - 6nH£sm(nn), (4.162) 

and n 2 = k 2 + 3K. Here, the coefficients c + and c_ are constants that depend on k. 

Tensor perturbations describe the propagation of gravitational waves, which are represented by the 
electric and magnetic Weyl tensor, E ao and H^. We specialize the solution (14.1601 ) to K = —1 and 
use ( 13.206b and H = £'/£ 2 to determine £ and H. On applying ( 14.801 ) and ( 14.811 ), we can determine the 
asymptotic behaviour of and H a t, in three asymptotic regimes: The first regime (n <^ 1 and krj -C 1) 
describes long wavelength perturbations that occur outside the particle horizon close to the singularity; 
the second regime (77 <C 1 and kn S> 1) describes perturbations within the particle horizon close to the 
singularity; lastly, the long term behaviour of the perturbation, when spatial curvature is significant, is 
described in the third regime (n S> 1). 





n <C 1 and kn <^ 1 


7] <C 1 and kn ^> 1 


7? > 1 


E oh IFF 2 


c + n 2 + c_?7~ 3 


A cos {nr/ + </>) 


Ae 


~ v cos (nr/ + (f>) 


H a b/H 2 


c_|_7/ 3 + C-rj~ 2 


Asin{nn + 0) 


Ae 


~ v sin{nn + 0) 



Here c±, A and </> are arbitrary constants that depend on k. In the first regime we can see that the 
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tensor perturbations do not represent gravitational waves since they have a power law dependence on r\. 
We depict the behaviour of E^/H 2 in figure l4~4l for three given choices of n. 

4.9 Perturbations of Radiation-Filled FL Models 

We now derive the general solution of the governing DEs for a radiation fluid with zero cosmological 
constant (A = 0) in terms of elementary functions. 

The solution for (ry) is obtained by choosing f3 = 1 in the general solution ( 14.1091) , (14.1111) and 
( 14.1 13b . Since the index has the value /j = |, these solutions are elementary functions, given by (IH. 14b - 
dH. 19b . Upon comparing the functional dependence in each case, we obtain the following unified form of 
the solution: 

V {k) (ri) = c + pW(ry) + cJD^\n), (4.163) 



where 



V {+ \ri) = H£sm(c s krj)-c s kcos(c s krj), (4.164) 
2?(-)(r ? ) = H£cos(c s ki]) + c s ksm(c s kr]), (4.165) 



and 



coth 7] Open FL (if = -1), 
IK - i 7]- 1 Rat FL (if = 0), (4.166) 
cot 77 Closed FL (if = +1). 

The coefficients c + and c_ are constants that depend on k. 

The solution for Xqa (77) is obtained by choosing j3 = 1 in the general solution (14.1091) , (14.1111) and 
(14.1 131) . Since the index has the value /i = |, these solutions are elementary functions, given by (IH. 14b - 
dH. 19b . Upon comparing the functional dependence in each case, we obtain the following unified form of 
the solution: 

*(k)(v) = Yli c + X U^ + c - X (k)(v)], (4-167) 

Xqs{ti) = Hism{ni]) — ncos(nr/), (4.168) 
X^{r]) = H£cos(nr]) + rasm(rar/), (4.169) 

(4.170) 



where 
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Figure 4.4: Tensor modes of the Hubble-normalized electric Weyl tensor E^/H 2 for an open FL 
background (K = — 1) with dust. From top to bottom, the eigenmodes in each plot are n = 10, n = 100 
and n = 1000. 
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and n 2 = k 2 + 3K. The coefficients c + and c_ are constants that depend on k. 

As with perturbations of a dust background, we specialize the solution (14.1631 ) to K = — 1 and use 
(13.2061) and iT = I' /I 2 to determine ^ and H. On applying (l4~80l > and P~80l> . we can determine the 
asymptotic behaviour of E a b and H a b in three asymptotic regimes. 





r\ -C 1 and fcr/ <C 1 


77 <C 1 and ki] I 


r] » 1 


Eab/H 2 


c+r/ 2 + c_?7 _1 


Arj cos(nrj + 0) 




_?? cos(n?y + (/>) 


H a b/H 2 


c + r/ 3 + c_ 


j4?7 sin(nr7 + (p) 




_,? sin(nr/ + (/>) 



Here c±, A and are arbitrary constants that depend on k. We depict the behaviour of E^/H 2 in 
figure R31 for three given choices of n. 
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Figure 4.5: Tensor modes of the Hubble-normalized electric Weyl tensor E^/H 2 for an open FL 
background (K = —1) with radiation. From top to bottom, the eigenmodes in each plot are n = 10, 
n = 100 and n = 1000. 
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4.10 Perturbations in Asymptotic Epochs of FL Models 

In this section we give approximate solutions of the governing ODEs for scalar and tensor perturbations 
in three asymptotic epochs for FL cosmologies. We also give the behaviour of the spatial gradients, 
kinematic quantities and Weyl tensor components for these approximate solutions. Further, to facilitate 
comparisons, we give the Bardeen (1980) gauge-invariant metric potentials for the approximate solutions. 
The three asymptotic epochs are as follows: 

Flat FL. The flat FL background is described by 

n = i, n A = o, n k = o, \£=(^L] , m = -. (4.m) 

\PJ V 

This solution approximates the dynamics of any FL cosmology with an initial singularity with cosmolog- 
ical constant and arbitrary curvature at early times, since these models satisfy 



Iim(n,nA,fifc) = (1,0,0). (4.172) 

Tj— >0 



(see table 



Milne. The Milne background is described by 

n = o, n A = o, n k = i, \i = e r >, m=i. (4.173) 

This solution approximates the dynamics of any open FL cosmology with zero cosmological constant at 
late times, since these models satisfy 

lim (fi, J2 Aj flfc) = (0, 0, 1). (4.174) 

Tj— >00 

(see table E2>. 

de Sitter. The de Sitter background is described by 

n = o, n A = i, n k = o, xt = -, m=~, (4.175) 

T) fj 

where fj is defined by ( 13.2191 ): 

fj = rj — r]f. (4.176) 
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This solution approximates the dynamics of any ever-expanding FL cosmology with cosmological con- 
stant and arbitrary curvature at late times, since these models satisfy 

lim (n, n A , n k ) = (o, 1, o). (4.177) 

V^Vf 

(see table O. 

Remark: The approximate solution for the case of flat FL with krj <C 1 (Table 14.21 ) have been given by 
Goode (1989). To the best of the author's knowledge, the other results are new. 
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4.10.1 The Flat FL Asymptotic Epoch 



The solutions of the governing equations in the fiat FL background are obtained by substituting ( 14.1711 ) 
into ( 14. 109b and ( 14.93b and are given by 



>V (fe )(r/) 



c + Jp + i(c s krj) + c_F^ + |(c s fery) 
c + J p+ i(kri) + C-YpMkr]) 



--0 
772 ^ 



(4.178) 
(4.179) 
(4.180) 



In the expressions for T>n : \(rj) and X^(tj) we can identify two asymptotic regimes, namely kr] <C 1 
and fery S> 1. The restriction krj <tl 1 means that the perturbation mode corresponding to k is outside 
the particle horizon, while kr] S> 1 means the perturbation mode is inside the particle horizon (see sec- 
tion 14-7- 1 b - The approximate solutions subject to these restrictions can be obtained from the asymptotic 
expressions for the Bessel functions in Appendix IH.4I 

These approximate solutions will also be valid for an open or closed FL background, for the following 
reason. If K = ±1 and kr] <C 1, it follows that r] must satisfy r] <^ 1 since k > 1 - see (IF. lib . (IF 32b and 
dF43l . It thus follows from ( 14. 172b that 

(4.181) 



— «1 



i.e. the spatial curvature is dynamically negligible in the regime krj -C 1. 

The approximate solutions and the resulting time dependence of the first order quantities are given in 
table \42\(kr] < 1) and table |431(fcr/ > 1). 
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Perturbations of FL models subject to kr] <C 1 and Qk/^m *C 1. 

Quantity Scalar Perturbations Vector Perturbations Tensor Perturbations 



W (fc) - - {r/ 2 - 2 ^} 



0~ab 

on 

H 




-2/3j 




•2/3| 




y- 


-2/3 1 


on 









-2/3| 









R]i 


{77 3 ,r7 2 - 


-2/3 1 




•2/3| 









V a H 
H 2 


{r/ 3 ,r/ 2 - 


-2/3j 


{r/ 2 ~ 


•2/3| 









on 

H 3 


{r/ 3 ,ry 4 - 


-2/3}f 












H 2 


{r/ 2 ,^ 1 " 


-2/3| 




•2/3j 


OA 


Y" 


-2/9 1 


H a b 

On 

H 2 







W' 


-2/3| 


{^ 3 . 


Y" 


-2/9 1 



* - - {r?- 2/3 } 

4 2) - - - {iy- 2/3 }* 



Table 4.2: This table describes perturbations of any FL model subject to kr] <C 1 and f2£./r2 m <C 1. For 
open and closed FL models the second restriction requires that which is necessarily satisfied on 

account of the first restriction, since k > 1. Recall that (3 = 2/ (37 — 2), where 7 is the equation of state 
parameter (see j4.95\» . 

' The second mode in this case is obtained directly from the exact solution and is zero in the case of dust. 
* The second mode in this case cannot be obtained from the asymptotic behaviour of . 
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Perturbations of FL models subject to krj » 1 and f^/^m <C 1. 



Quantity 




Scalar Perturbations (7 7^ l)t 


Vector Perturbations 


Tensor Perturbations 






r/ 1 ^ 13 sm(c s ki] + </>) 


- 


- 








^2-2/3 




W 








r]~P sm(kr/ + 0) 


Cafe 




?7 2 ~^ cos(c s kri + 0) 




sin(fc?7 + 0) 


~H 












ff/i 




ry 2- ^ sm(c s kr] + <^>) 







V a F 
H 2 




r/ 3- ^ cos(c s krj + 0) 


^2(3 





V a ^R 

H 3 




t? 3 "^ sin(c s /cry + (ft) 








E a b 

H 2 




r/ 1 ^ 13 sm(c s krj + 




rf~^ cos(kr) + (ft) 


H a b 
H 2 







^2-2/3 


r] 2 ~@ sin(fc?7 + (ft) 



rj P cos(kr] + (ft) 

Table 4.3: This table describes perturbations of flat FL with kr\ 3> 1, and open and closed FL with 
kr] ^> 1 and r] <C 1. For open and closed FL models the second restriction requires that ^ <C 1. The first 
restriction then requires that k is sufficiently large. Recall that (3 = 2/(37 ~~ 2), where 7 is the equation 
of state parameter (see H4.95\> ). 

' For the case of 7 = 1 see table \4.2\ 



^ Tj sm(c s ki] + (f>) 

* - - T]~ 213 
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4.10.2 The Milne Asymptotic Epoch 

The governing DEs in the Milne background are obtained by substituting ( 14.1731 ) into ( 14.971 ), ( 14.901 ) and 
(14.1011) and are given by 

V'[ k) +uj 2 V {k) = 0, V [k) = e- (1 -^73 (fc) , (4.182) 

W ( , fc) + |(/3-l)W (fc) = 0, (4.183) 
X'( k) + n 2 X {k) = 0, X {k) = e-«X {k) . (4.184) 



where 



holds for scalar perturbations and 



to 2 = c 2 k 2 - (1 - i) 2 (4.185) 
n 2 = k 2 -3 (4.186) 



holds for tensor perturbations. For scalar perturbations there are two cases, namely lo 2 > and lo 2 < 0. 
The second case includes dust (i.e. c 2 = 7 — 1 = 0, (3 = 2), while the first case includes radiation (i.e. 
7 = |, (3 = 1). In addition, the condition lj 2 > is satisfied for all perfect fluids other than dust (i.e. 
1 < 7 < 2) provided that the wave number k is sufficiently large, i.e. k > k cr i t . We will restrict our 
considerations to the cases lo 2 > and dust. 
The solutions are 

_ , , f Ae~^ 1 ~^ 7? sin(a;77 + if 7 > 1, k > k crit , 

T> {k) ( v ) = { K 1 Vh 1 (4.187) 

[ c + + c_e if 7 = 1, 

W (fc) fa) = Ce- 2{1 -3)", (4.188) 

Af (fe )(T7) = Ae -7? sin(nr/ + ^), (4.189) 

where A is an arbitrary constant amplitude. 

In order to calculate the approximate expression for E a b we need an approximate expression for 
subject to T] S> 1. A simple calculation using (13.2061) and (13.2071) yields the approximate asymptotic 
behaviour 

= exp('-|7 ? Wl. (4.190) 

The approximate solutions and the resulting time dependence of the first order quantities are given in 
table f 
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Perturbations of Milne. 



Quantity 



V 



(*) 



(k) 



(k) 



Scalar Perturbations 



Vector Perturbations Tensor Perturbations 



7=1 



7 > 1, k 2 > k 2 crit 



{l,e- r ?} e" (1 "^ ),? sin(a;r ? + 



-2(l-i)r? 



e v sin(nr/ + 0) 



H 
H 

Vg/t 

Hfi 

# 3 

Egb 

H 2 

Hnh 



H 2 



- 

{l.e - "} e 



^ sin(a;r/ + 




{r/e v ,e V Y e (1 sin(o;r/ + 



-»7 „-2»n 



e ^ 0^ sin(ujr] + 







e -2(l-> 

e - 2 a-^ 



-2,, 



-2,, 



e v sin(nry + < 








e _,? cos(nr] + 



if, 



(2) 



{e-^,e- 2r? } e 



sin(u>7/ + 



-27, 



e ^ sin(nr/ + 



Table 4.4: TMs tabie describes perturbations of open FL models with rj S> 1 and r2fc/f2 m 3> 1 and 
hence krj 3> 1. Recall that (5 = 2/(37 ~~ 2), where 7 is the equation of state parameter (see H.95\l ). 

' In the case of 7 = 1 the growing mode is obtained from the growing mode for dust in an open FL background, given by 
f4lM . 
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4.10.3 The de Sitter Asymptotic Epoch 



The governing DEs in the flat de Sitter background are obtained by substituting ( 14.1751 ) into (14.971 ), ( 14.901 ) 
and (14.1011) and are given by 



dk 2 - \ (1+ i 



ft 



ft J fj 2 



%0 



Vtf M + 2[4-l 



'(*) 



ft 



(fc) 



1— - A 



(4.191) 

(4.192) 
(4.193) 



where the prime (') denotes differentiation with respect to fj and n 2 = k 2 + 3 - see (1F.43I ). The solutions 
to these equations are given bvF^l 



fj 2 f> 



C-Ji,i(c s kf]) + c + Yi,i(c s kfj) 

2^/3 2 



c^fj z + c+, 



if 7 = 1 



-2-4 



c ?7 , 
c_7)sin(nf/) + c+?? cos(nfj). 



(4.194) 

(4.195) 
(4.196) 



As in the case of the flat FL background, we can distinguish two regimes, depending on the wave 
number, i.e. kfj <C 1 and kfj 3> 1. These results are valid also for K = ±1, but require that fifc/Jly *C 1, 
which holds for 77 <C 1 on account of (14.1751 ). 

In order to calculate the approximate expression for E a b we need an approximate expression for U 
subject to fj <C 1. A simple calculation using ( 13.2251 ) and (13.2261 ) yields the approximate asymptotic 
behaviour 

ft~f/ 37 . (4.197) 

The approximate solutions and the resulting time dependence of the first order quantities are given in 
table @3] (kfj < 1) and table 1431 (kfj > 1). 



4 Observe that ( |4.191| > is Bessel's DE l lH.lt in normal form with a = c B k and /i 2 — ~ = 4(1 + 4), i.e. ^1 = | + 4. 
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Perturbations of de Sitter subject to kfj -C 1 and Qk/^A <C 1. 

Quantity Scalar Perturbations Vector Perturbations Tensor Perturbations 

^(fc) ~ - - {v,v 2 } 



°~ab r~3-| ~2t f ~3-| r ~2 ~3i 

s 



iJ 2 
H 2 



{fj ?,ff} 7? /3 



H 

Hp, 

- ^ 4 

H 3 

^ /*3 S 2(2+ih ~3 r~4 ~3 



{?7 ^jf} 



rf 



2 . 



Table 4.5: This tabie describes perturbations of any FL model that is future asymptotic to de Sitter, 
subject to kfj <C 1. Recall that (3 = 2/(37 ~~ 2), where 7 is the equation of state parameter (see if4.95D ). 

' Due to cancellation of higher order terms, the growing mode must be obtained from the exact solution (4.194\) . 
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Perturbations of de Sitter subject to kfj ~> 1 and Qk/^A *C 1. 

Quantity Scalar Perturbations (7 7^ 1) Tensor Perturbations 

£>(fc) ^ ?? 3 sin(c s /c?? + 4>) 

X {k) - - fjsm(nfj + 4>) 



fj 2 sin (nfy + > 








fj 3 cos(nfj + 
ff sin(nfj + < 



Cafe 




~2-i 
7/ 8 


H 










H 










2-i 

7) 13 






















~3-i 






E a b 




if** 


H 2 






H a b 






H 2 






<Z> H 






tt{2) 







7ycos(nr/ + 0) 

Table 4.6: This table describes perturbations of any FL model that is future asymptotic to de Sitter, with 
kfj ^> 1 and fj <C 1. In the case of scalar perturbations with 7 = 1 or vector perturbations, the behaviour is 
identical to that given in table \4.5\ Recall that f} = 2/(37 ~~ 2), where 7 is the equation of state parameter 
(see K95\> ). 
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4.11 Discussion 

We now summarize the principal features of linear perturbations of FL cosmologies. 

i) For each perturbation type the complete system of linearized evolution and constraint equations 
is reduced to a single ODE for the harmonic modes, together with expressions for the Hubble- 
normalized quantities ( 14. 1 14| >-( l4. 116b - Solutions of the governing DE for vector perturbations can 
always be determined explicitly in terms of £ according to ( 14.93b . 

ii) When written in normal form using conformal time, the governing DEs for scalar and tensor per- 
turbations are essentially the same if A = 0, i.e. they are governed by the same potential U with a 
different frequency to, but differ significantly if A > 0. 

iii) We have given a unified derivation of the known explicit solutions of the governing DEs for scalar 
and tensor perturbations in (14. 109b - d4. 1 131 >. Solutions are given in terms of Bessel functions if 
A = and K = with 7 arbitrary. Solutions are given in terms of Legendre functions if A = 
and K = ±1 with 7 arbitrary. In the case of dust (7 = 1) and radiation (7 = |), solutions for all 
perturbation modes can be written in terms of elementary functions. 

We now give a brief overview of the behaviour of scalar and tensor perturbations when the matter 
content is dust or radiation and possibly with a cosmological constant. 

Scalar perturbations, as described by the spatial density gradient V a , have a growing mode and a 
decaying mode in the regime krj <C 1. For pressure-free matter these two modes persist for all rj. The 
growing mode is bounded in the Milne regime, i.e. if A = and K = —1, and in the de Sitter regime, i.e. 
if A > 0. 

For radiation (or more generally if the barotropic index satisfies w > 1), scalar perturbations are 
wavelike in the regime c s kr] 3> 1 (interpreted as sound waves with speed of propagation c s ). The spatial 
density gradient is wavelike with constant amplitude in the Milne regime (unbounded if w > |) and in 
the de Sitter regime (with kfj <C 1) it has a growing and unbounded mode and a decaying mode (more 
generally this occurs if w > 0). 

Tensor perturbations, as described by the electric and magnetic Weyl curvature, have a growing mode 
and a decaying mode in the regime krj -C 1. In the regime krj S> 1 where £7&/O m -C 1 the behaviour is 
sinusoidal, with power-law amplitude. 



4.11. Discussion 



123 



In tables |4~2ll4.6| we summarize the time-dependence of the Hubble-normalized kinematic quantities, 
Hubble-normalized Weyl tensor components and metric potentials in three asymptotic epochs for FL 
cosmologies. Tables |4~2l and 14. 3 1 show clearly that the flat FL model is unstable, both with respect to long 
wavelength and short wavelength perturbations. This instability is not revealed by the gauge invariant 
metric potentials. Tables |4~4l - 14.61 show that spatial curvature has a stabilizing effect, as does the presence 
of a cosmological constant. 



Evolution and Constraint 
Equations 



As discussed in section 1 1 .2.61 when the Ricci identities (11.501) and Bianchi identities (11.511) are applied to 
the Einstein field equations dl.2| ), we obtain a set of evolution and constraint equations for the kinematic 
quantities and Weyl tensor components. For an outline of the derivation of these equations see Ellis (1973, 
p24-30) and Ellis and van Elst (1998, plO-12). 

Projected Field Equations. The Einstein field equations (11.21 ) can be rewritten as 

Rab = T ab ~ \Tg ah . (A.l) 

We contract dA.ll) with h ab and u a and use ( 11.26b in order to obtain the following relations: 

R = fx-3p, (A.2) 
RabU a u b = |(M + 3p), (A.3) 
R ab u a h b c = -q c , (A.4) 

R {ab) = ^ab, (A.5) 

(see Ellis (1973), p24). Using (11.301 ) we also obtain two important relations relating the Riemann curvature 
tensor and Weyl tensor components: 

Racbdu c u d = E ab - \R {ab) + \h ah R cd u c u d , (A.6) 
±rj ac ef R efbd u c u d = H ab + ± Vabcd R c e u d u e . (A.7) 
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Ricci Identities 




H = -H 2 + UV a u a + u a u a ~2a 2 + 2u 2 ) - i(u + 3p), 


(A. 8) 




(A.9) 


U/ a \ = —2Hu) a + a ab oj b - icurl(u a ), 


(A. 10) 


= VV afe -2V aJ ff-curl(cu a )-2e afec nV + 9a , 


(A. 11) 


= V a a; a -u a c; a , 


(A. 12) 


= Hob- 2u {a u; b) - V( a u b) - curl(cr afe ). 


(A. 13) 


The twice-contracted Bianchi identities (11.531) lead to the following evolution equations: 




fi = -3H( f i+p)-V a q a -2u a q a -a a b TT b a , 


(A. 14) 


q{a) = -^nq a - V a p - (fl+p)Ua - V -K ab 




-u b TT ab - a ab q b + e abc to b q c . 


(A. 15) 


We are primarily interested in the case of a perfect fluid, where these equations reduce to 




Twice-Contracted Bianchi Identities (Perfect Fluid Source) 




A = -3H(n + p), 


(A. 16) 


= (fi + p)u a + V a p. 


(A. 17) 



For simplicity, the equations obtained from the full Bianchi identities (11.511 ) are presented here in the 
case of a perfect fluid source. 
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Bianchi Identities (Perfect Fluid Source) 

E {ab) = -3HE ab + 3<7 (a c £ 6)c - |(M + PWb + curl(i7 a6 ) 

+V(a^ e <-^))' (A-18) 

tf (a6) = -3^^ + &r (a c # 6)c - curl^) 

"W^f + ^O' (A - 19) 

= V 6 ^, + 3^tf a6 + e abc a be H e c - \V a ^ (A.20) 

= V b H ab -3co b E ab -e abc a be E e c -(fi + p)u; a . (A.21) 



Spatial Gradients of Zero-Order Quantities. In order to derive the governing equations for linear 
perturbations of FL, it is necessary to derive evolution equations for the spatial gradients of the matter 
density // and the Hubble parameter H, defined by (14.81) and d4.9b . By propositions 12. 1 .41 and 12. 1 . W these 
quantities are zero in the background and hence gauge-invariant by the Stewart-Walker lemma. 
We make use of the fact that the Raychaudhuri equation dA.81 ) can be rewritten as 

H =\T -\(ix + p), (A.22) 

where 

T = -3H 2 + fi + A + V a u a + ii a u a - 2a 2 + 2uA (A.23) 

In order to proceed, we require an identity which permits the interchange of the ' derivative and the 
spatial derivative V a when applied to a scalar field /. A straightforward calculation using d 1 - 18b yields 

(V a /)i = -HV a f + (V a + u a )f - (a b + io a b )V b f, (A.24) 

for any scalar / (also see equation (92) in Ellis et al. (1990)). 

We now derive evolution equations for X a and Z a , following Ellis and van Elst (1998, p56). Substi- 
tuting / = \i into (I A. 241 ) and using (14.8b . d4~9l ), (1A.16I ) and (1A.17I) yields the evolution equation for X a . If 
we instead choose / = H in dA.241 ) and use d48l ), ( 14.91 ), dA.221 ), dA. 16b and dA.17b we obtain the evolution 
equation for Z a . The results of these derivations are given as follows: 
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Spatial Gradients of Zero-Order Quantities 






X {a) = -AHX a -(a b + oj b )X b - 


- (1 + w)/j,Z a , 


(A.25) 


Z {a) = -3HZ a -(a b + u; b )Z b - 


\x a 




+T u a + V a (V b ii b + u a u a 


- 2a 2 + 2u 2 ). 


(A.26) 



Spatial Curvature 



In this appendix we introduce the 3-curvature when the fundamental congruence has non-zero vorticity. 
This analysis leads to a generalization of the Gauss equation (I1.37I ). the 3-Ricci tensor dl.391 ) and 3-Ricci 
scalar. 

B.l The Generalized Gauss Equation 

We note that the projected covariant derivative V a , defined in (11.111) . is a non-commutative operator on 
the space-time manifold when vorticity is present, i.e. u> a / 0. After a short calculation using (11.51 ). (11.1 II ) 
and dl.181 ), we obtain 

V a V b f - V b V a f = -2u ab f. (B.l) 

We now derive a similar identity for tensors with one or more index. 
On projecting (11.181 ) orthogonal to u a and applying (11.111 ). we obtain 

V b u a = k ab , (B.2) 

where k ab is defined by 

Kb = u ab + <r ab + Hh ab , (B.3) 

(see Ellis et al. (1990), eq. (76)). 

Upon performing a similar computation as in (IB . 1 b and using (11.501 ). we obtain the expression 

V a V 6 X c - V b V a X c = -2co ab X {c) R abcd X d , (B.4) 
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B.2. The Spatial Gradient of the 3-Ricci Scalar 
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where ^'R a bcd is defined by 

^'Rabcd = {Rabcd)± ~ k ac k bd + k ad k bc , (B.5) 

(see Ellis et al. (1990), eq. (79)). We note that this definition reduces to the 3-Riemann tensor given in 
(11.371 ) when u a is irrotational. A similar identity to (IB .41) can be obtained for tensors with two or more 
indices, generalized according to 

V[ a V b ]5 c rf... e / = -UJ a bS±cd-ef + \{ ^R S cbaT S d ... e f H h ^R s fb a T cd ..J). (B.6) 

On contracting (1B.5I) . we can obtain an expression for the spatial Ricci scalar and trace-free spatial 
Ricci tensor in terms of the spacetime quantities from dB.51) , incorporating the Weyl tensor through d 1.30b 
and CE32]>: 

(3) R = R + 2R bd u b u d -m 2 + 2a 2 - 2cu 2 , (B.7) 
^S ab = E ab + \(h^h b — \h cd h ab )R cd — Ha ab — Huj ab 

- 1 (a 2 - u 2 )h ab + a a c a cb + io^iv cb + ^a a cb + o^ ch . (B.8) 

Using the field equations (11.21) in the form (1A.2I) - (1A.5I ) to eliminate R and R ab , we obtain 

(% = -6H 2 + 2fi + 2a 2 - 2u 2 , (B.9) 
( % a6 = E ab + lir ab - Ha ab - Hu ab -l(a 2 - u 2 )h ab 

+0-aCTcb + ^a^cb + ^a a cb + ^a U cb ■ (B • 1 0) 

B.2 The Spatial Gradient of the 3-Ricci Scalar 

The spatial gradient of the 3-Ricci scalar is defined by 

TZ a = £V a ^R. (B.ll) 

An evolution equation for 1Z a is obtained in the same manner as for V a and Z a , defined in ( I4.35I ) and 
(14.37I) . respectively. If we choose / = ^R in (1A.24I ) and use d4~9b . (IA.22I ). (IA.16I) and (IA.17I ) we obtain 

7t (a> = -2^ a -( C j a f ' + u; a f ')^-|TZ a 

-4£Tu a + £V a (V b u b - 2a 2 + 2uj 2 ). (B.12) 



B.2. The Spatial Gradient of the 3-Ricci Scalar 
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Further, we can obtain an auxiliary equation describing the evolution of 1Z a from dB. 12b and ( 14. 19b . 
presented here for completeness: 

4K c 2 ( -n K \ 

K {a) = -2HK a - -pZ a + 4^^^y [V 2 + -jv a + 2Ac 2 H 2 £ curl(cu a ). (B.13) 



Solutions of the Friedmann 

Equation 



In this appendix we list the solutions of the three first order quadratic DEs that arise by transforming the 
Friedmann DE. In the first two cases, the solution is uniquely determined by requiring ^ > and the 
initial condition 

L(0) = 0. (C.l) 

In the third case we require that 

L(0) = 1, and L'{T) > for T > 0. (C.2) 



Differential Equation 
Casel f^Y = l + 2aL + L 2 , 

Case 2 (j^Y = 1 + 2aL - L\ 
CaSe3 iff-*- 1 ' 



Solution 



L = sinh T + 2a sinh 2 ( \T) . (C.3) 



L = sinT + 2asin 2 (iT). 



L = coshT. 



(C4) 



(C5) 
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C.l. Solutions with a Stiff Fluid 
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There are four known cases where the Friedmann equation can be transformed into the quadratic form 
( 13.1951 ) via a change of variable: 

i) 1 -fluid plus spatial curvature. 

ii) 2-fluid flat FL. 

iii) 2-fluid FL plus arbitrary curvature. The equation of state parameters must satisfy 

7i = | + 2A 7 , 72 = | + A 7 , (C.6) 

where A7 is constant. 

iv) 3-fluid FL. The equation of state parameters must satisfy 

7i > 72 > 73, 7i - 72 = 72 - 73 = A7, (C.7) 

where A7 is constant. 

C.l Solutions with a Stiff Fluid 

There are three known cases where we can obtain a solution to the Friedmann equation where one fluid is 
a so-called stiff fluid, i.e. 7 S = 2. In this section we briefly present each of these solutions. 

SD -Universes. We choose the conformal parameter to be A = Xd, which, on using (13.1131) . leads to a 
single mass parameter defined by 

/ A., x ~ 4 
m s = I 

In this case the Friedmann equation reads 



m s =(y) (C.8) 



, dt 

We apply the change of variable 



^ =(\ s £)- 4 + (\ d £)- 1 . (C.9) 



C = ^£ 3 , T = 3^t, (CIO) 
A A 



obtaining 



d£_ x 2 
dT 



l + C. (C.ll) 



C.l. Solutions with a Stiff Fluid 
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Subject to the IC £(0) = 0, the solution is 

£(T) = ±T 2 + T, (C.12) 

or equivalently 

t(t) = \~ 1 [f(At) 2 + 3V^(At)] 1/3 . (C.13) 

S'-RC-Universes. We choose the conformal parameter to be A = X r , which, on using (13.1131 ). leads to a 
single mass parameter defined by 



In this case the Friedmann equation reads 

■di\ 2 

dri) 



m,= (y] • (C.14) 



(X s £y A + (X r £)- 2 - K. (C.15) 



We apply the change of variable 
obtaining 



dT ) \/rn 
Subject to the IC £(0) = 0, the solution is 



£ = A 2 £ 2 , T = 2r ) , (C.16) 



<ICV \ + -^C-KC 2 . (C.17) 



£(t) = A" 1 yW s S K {2r,) + S K ( V ) 2 ] 1/2 . (C.18) 

S'DA-Universes. We choose the conformal parameter to be A = Aa, which, on using (13.1 131) . leads to 
two mass parameters defined by 

™s=(jY\ md =(^y\ (C.19) 
In this case the Friedmann equation reads 

(ID 2 = (a ^ } ~ 4 + + {xi)2 - (c - 2o) 

We apply the change of variable 

C = X 2 S X£ 3 , T = 3Xt, (C.21) 



C.2. Existence-Uniqueness for the Friedmann Equation 



134 



obtaining 



=1 + 



dry 



Subject to the IC £(0) = 0, the solution is 

£{t) = X- 1 [^Isinh(3At) +m d sinh 2 (|At)] 1/3 



(C.22) 



(C.23) 



C.2 Existence-Uniqueness for the Friedmann Equation 



In this section we prove existence-uniqueness of solutions to the Friedmann DE (13.1 151) . given by 

2 n 



y i=l 



37i+2 _ K 



for any choice of parameters > 0, i = 1, . . . ,n — 1 and subject to the conditions 



L(0) = 0, and ^ > 0, 
dT 



We can rewrite (1C.24I ) as 



—\ = r-( 3 ^- 2 ) 
dT J 



i=i 



Since 71 > 72 > • • • > 7n and 71 > |, this DE is of the form 

2 



with 



F(0) = mi > 0. 



Since ^ > 0, we can rewrite (1C.27I ) as 



dL ~ F(L)V2 ' 
The general solution to dC.29l > is then, for all > 0, 



L > 0. 



T = G(L) + C, 



(C.24) 



(C.25) 



(C.26) 



(C.27) 



(C.28) 



(C.29) 



(C.30) 



where G(L) is the antiderivative of the RHS of dC.291 ). The requirement that T = when L = fixes C, 
and hence gives a unique solution. 



CONFORMAL TIME 



The conformal time variable rj is defined in terms of clock time t and the length scale £ via 

£4 

Unlike clock time, conformal time is not necessarily unbounded into the future for ever-expanding 
models. In fact, the presence or absence of certain types of fluids determines the behaviour of the confor- 
mal time variable. We present the following theorem: 

Theorem D.l Let (M, g, u) be an ever-expanding n-fluid FL cosmology whose length scale satisfies 
£(0) = and where each fluid satisfies a linear equation of state A3. 1 71 ) with equation of state parameters 
7i ordered according to 71 > 72 > • • • > J n . It follows that r] satisfies the following constraints: 

(a) rj is bounded below iff 71 > |. 

(b) rj is bounded above iff '72 < |. 

Proof. It follows from the definition of an ever-expanding model (see section [J. 1.6b that 




(D.2) 



Using the chain rule we expand ( ID. lb so as to obtain 



di 



-1 -1 



di 



£ 



(D.3) 



It then follows from (ID.2b and (ID.3I) that rj is a strictly increasing function of £ for t > 0. 
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Let T)o and 77* be the value of conformal time at £q and I*, respectively. We now integrate dD.3b over 
(t/o i V*) so as to obtain 

'dt 



r]*-r]o = F(£*) 



in 



dt 



(DA) 



Since 77 is strictly increasing, we note that 77* is bounded above if and only if F(£*) is bounded as — > 
+00. Similarly, we note that 77* is bounded below if and only if F(£*) is bounded as I* — > + . We now 
prove the two results separately: 

(a) Assume that condition (a) is not met. It then follows from ( 13.1961 ) that 3 a constant C > such that 

d£ 
dt 



< C V t < to. The integral dD.41 ) then takes the form 



V* = % 



d£\ 
dt) 



< Vo 



C 



(D.5) 



which is unbounded. 

If condition (a) is met it then follows from (13.1961 ) that 3 constants C > and e > such that 
§ > C£~ £ V t < t - The integral (E3> then takes the form 



V* = m 



dt 



- 1 ■ 1 eft 1 

de>m ~cJ 



(D.6) 



which is bounded. 



(b) Assume that condition (b) is not met. It then follows from (13.1961) that 3 a constant C > such that 



^ < C V f > t - The integral dD.41 ) then takes the form 



d£ > r] + lim — / — 



(D.7) 



which is unbounded. 

If condition (b) is met it then follows from (13.1961) that 3 constants C > and e > such that 
^ > C£ e y t < t . The integral (ID.4I ) then takes the form 



77* = 770 + lim 



V eft 



-1 



d£ > 770 + lim 



1 



(D.8) 



which is bounded. 

Thus, we conclude that the constraints hold, completing the proof. □ 
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We note that this theorem also holds in the case of closed FL as long as ^ > C, where C is a constant 
satisfying C > 0. This constraint prevents models which recollapse or whose length scale tend to a 
constant as t — > +00. 



Appendix E 



Covariant Differential Identities 



In this appendix, we present a list of useful differential identities involving the derivative operators used 
in this thesis. This list represents an extension of the list given in Ellis et al. (1990) Appendix A by 
incorporating identities involving the Laplacian and curl operators. 

Zero-order 3-curvature quantities for FL models: Before proceeding, we must first obtain zero order 
expressions for the background curvature. It quickly follows from (13.111 ). (11.401 ) and (11.471 ) that, to zero 
order, the background 3-curvature quantities satisfy the following relationships: 



i 2{3) Rabcd = K (h ac h bd -h ad h bc ), (E.l) 

£ 2{3) Rac = 2Kh ac , (E.2) 
f = 6K. (E.3) 

We now present a list of linearized identities, derived using (11.7I ). (11.111) . dl.16b - dl.17l ). (IB. II) and 
dB.6b . In all of the following relations, we shall assume that V a /, X a and T ab and all their derivatives are 
first order quantities. Further, we assume that T ab is symmetric and trace-free. 

First-order time derivatives: The following set of linearized identities are used when interchanging 
the derivative along the fundamental congruence (') with V, curl(-) or V 2 : 
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V a (/) - (V a /)± = HVaf-fUa, (E.4) 

V a (X 6 ) - (V a X 6 )-j_ = #V a X 6 , (E.5) 

V a (T 6c ) - (V a T 6c )l = tfV a T bc . (E.6) 

£V a (/) = (£V a fy ± -£fu a , (E.7) 

£V a (X 6 ) = (£V a X 6 )i, (E.8) 

^V a (T 6c ) = (^V a r fec )- ± . (E.9) 

curl(X a ) - (curlpf a ))- = #curl(X a ), (E.10) 

curl(T a6 ) - (curl(T a6 ))- = Hcur\(T ab ). (E.ll) 

^curl(X a ) = (£cur\(X a ))\, (E.12) 

£cur\(f ab ) = (tcur\(T ab ))\. (E.13) 

V 2 (X a ) = (V 2 X a )\ + 2HV 2 X a , (E.14) 

V 2 (T a6 ) = (V 2 T ah y ± + 2#V 2 T ab . (E.15) 

£ 2 V 2 (X a ) = (£ 2 V 2 X a )\, (E.16) 

e 2 v 2 (f ab ) = (e 2 v 2 T ab y ± . (E.n) 



Laplacian identities: The following linearized identities are used when interchanging the Laplacian 
operator V 2 with V or curl(-): 



£ 2 V a X7 2 f 
£ 2 V a V 2 X b 
£ 2 V a V 2 T bc 



= l 2 V 2 V a f-2KV a f-2fcur\(ta a ), 

= £ 2 V 2 V a X b - 2K [V a X b + V b X a - h a b V c X c ] , 

= £ 2 V 2 V a T bc - 2K \v a T bc + V b T a c + V c T b a - h a b V s T sc - h a c V s T bs 



(E.18) 
(E.19) 

(E.20) 
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£ 2 V a V 2 X a = e 2 V 2 V a X a + 2KV a X a , (E.21) 
e 2 V b V 2 T ab = £ 2 V 2 V b T ab + 4KV b T ab , (E.22) 

V 2 curlX a = curl(V 2 X a ), (E.23) 
V 2 curlT afe = curl(V 2 T afe ). (E.24) 

Curl identities: Finally, the following linearized identities are used when interchanging the curl operator 
curl(-) with V: 



curl(V /) = -W, (E-25) 

curl(V a X 6 ) = V (a curl(X 6) ), (E.26) 

cur\(V {a X b) ) = IV (a curl(X b) ), (E.27) 

V a curl(X a ) = 0, (E.28) 

curl(V b T ab ) = 2V 6 curl(T ab ), (E.29) 



Further, we obtain the following linearized identities for the curl of a curl, in the linearized context: 



^ 2 curlcurlX a = 
^ 2 curl curl T ab = 



-£ 2 V 2 X a + £ 2 V a V b X b + 2KX a , 
-£ 2 V 2 T ab + |£ 2 V (a V c T 6>c + 3KT ab , 



(E.30) 
(E.31) 



Appendix F 



Covariant Harmonics 



In this appendix, we present an overview of general identities and results related to the covariant harmon- 
ics. We extend the results of Bruni et al. (1992) Appendix B by introducing the notion of associated vector 
and tensor harmonics. 

(n) 

As in Hawking (1966), we define the spatial harmonics Q ab ... c as eigenfunctions of the spatial Laplace- 
Beltrami operator: 

V 2 Qi5.. c + ^Q<£.. c = 0, (F.l) 

where k is known as the wavenumber of the harmonic. Further, using (IF II ) we impose that the 
quantities are constant in time: 

Q ( aL = 0- (F-2) 

See Bruni et al. (1992, Appendix B) for the relation between the metric harmonics Y and the covariant 
harmonics Q. Hereafter, as we are working in the geometrical formalism, the use of the term harmonics 
will refer to the covariant harmonics defined by (IF II) and (IF2I) . For the problem at hand, we only require 
scalar harmonics, vector harmonics and two-index tensor harmonics (which, for purposes of brevity will 
be referred to as tensor harmonics). 

We briefly discuss notation for harmonics and quantities derived from harmonics. The superscript of 
the harmonic quantity Q(°> k ) denotes that it is a scalar harmonic. Similarly, we use 1 and 2 to denote vector 
and tensor harmonics, respectively. It is common in the literature to drop the wavenumber k when writing 
a harmonic quantity, so writing Q(°' k ^ as instead. The wavenumber of the harmonic being used is 
instead implicitly defined by the equation where the harmonic appears. Further, here and elsewhere, we 
shall use Q to denote the whole set of harmonics and the standard quantities constructed from Q: we have 
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F.l. Scalar Harmonics 
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scalar harmonics Q(°\ , vector harmonics Qa , Q ab and tensor harmonics Q ab - 

We now give the definitions and properties of each of these three classes of harmonic types. 

F.l Scalar Harmonics 

Scalar harmonics are defined as solutions of 

V 2 Q = ~Q, (F.3) 

that satisfy the time-independence constraint 

Q = 0. (FA) 
Each scalar harmonics Q defines a vector and a trace-free symmetric tensor via 

Q a = -^V Q, (F.5) 

and 

Qab = p V fe V a Q + \h ab Q. (F.6) 
Since it follows from (1F.4I ) and (IB. II ) that Q satisfies 

V [a V 6] Q = 0, (F.7) 

it can be quickly shown that Q a and Q a b satisfy the curl-free constraint, 

curl(Q a ) = 0, cur\(Q ab ) = 0, (F.8) 

congruence-orthogonality, 

u a Q a = 0, u a Q ab = 0, u b Q ab = 0, (F.9) 

and time invariance, 

Qa = 0, Q ab = 0. (F.10) 

We are restricted in our choice of k 2 to the eigenvalues determined by Lifshitz (1946), which in turn 
are determined by the spatial curvature according to 

9 9 f n > if if = -1,0, 

k = n - K < (Fll) 
n = 3, 4, . . . if A" = +1. 



F.2. Vector Harmonics 
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We note that it is common to include n = 2 as a possibly harmonic mode in the case of K = +1, but it 
can be shown via the constraint equation ( 14.291 ) and identity (IF. 161) that this mode is unphysical. In fact, 
this mode is gauge-dependent and associated with the freedom of choice of the velocity field u a . 

Some important identities involving the vector quantity Q a can be derived using the differential iden- 
tities in Appendix |E] and dF.3b - dF.5D : 

lV a Q a = kQ, (F.12) 
£ 2 V 2 Q a = -(k 2 -2K)Q a , (F.13) 
£V b Q a = -k{Q ab -\h ab Q), (F.14) 
^ 2 V [a V fe] Q c = \K{h ac Q b -h bc Q a ). (F.15) 

Similarly, some important identities involving the tensor quantity Q ab can be derived using the differ- 
ential identities in appendix |E] and dF.3b - dF76b : 

tV b Q ab = lk-\k 2 -3K)Q a , (F.16) 
£ 2 V a V b Q ab = Uk 2 -3K)Q, (F.17) 



t 2 V h V c Q ac = -|(fc 2 -ZK){Q ab -\h ab Q), (F.18) 



£ 2 V 2 Qab = -(k 2 -6K)Q ab , (F.19) 

^ 2 V[ a V b] Q cd = \K[(h ac Q bd - h bc Q ad ) + (h ad Q bc - h bd Q ac )] . (F.20) 

F.2 Vector Harmonics 

Vector harmonics are defined as solutions of 

k 2 

V 2 Qa = —pQa, (F.21) 

where Q a is divergence-free and time-independent, 

V a Q a = 0, Q a = 0. (F.22) 
Every vector harmonic Q a defines an associated vector harmonic P a via 



P a = £cur\(Q a ), 



(F.23) 



F.2. Vector Harmonics 
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which can be shown to satisfy dF.2 lb and ( IF. 221 ) using dE.231 ), dE.28b and ( IE. 12b . On taking the curl of 
(IF. 23b and applying (IF.21b and (IE.30b . we obtain Q a in terms of P a , 

Q a = (k 2 + 2K)- 1 £ curl(P a ). (F.24) 

The vector harmonics Q and associated vector harmonics P then define trace-free symmetric tensor quan- 
tities via 

Qab = -^(aQb), (F-25) 

and 

Pab = ~V {a P b) . (F.26) 
These tensor-valued vector harmonics are defined in order to satisfy zero second divergence, 

V a V b Q ab = 0, V a V b P ab = 0, (F.27) 

congruence orthogonality, 

u a Q ab = 0, u b Q ab = 0, u a P ab = 0, u b P ab = 0, (F.28) 

and time invariance, 

Qab = 0, P ab = 0. (F.29) 

The tensor-valued vector harmonics are related using (IF. 25b . (IF.26b and (IE.24b . as follows: 

^curl(P afe ) = l(k 2 + 2K)Q ab , (F.30) 
£cur\(Q ab ) = \P ab . (F.31) 

As with scalar harmonics, we are restricted in our choice of k 2 to the eigenvalues determined by 
Lifshitz (1946), according to 

9 n f n > if K = -1,0 

k =n — 2K < ' . (F.32) 

[ n = 3,4, ... ifiv" = +l 

Again, it is common to include n = 2 as a possibly harmonic mode in the case of K = +1, but it can be 
shown via the constraint equation (IA.2 lb that this mode is unphysical. 



F.3. Tensor Harmonics 
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Using the differential identities in Appendix |E] and (!F.21| )-( fF2"3l we can derive the following identity 
for the vector harmonics Q a and associated vector harmonics P a (in this case, we may replace Q a with 
P a to obtain a second identity): 

^ 2 V [a V b] Q c = \K{h ac Q b - h bc Q a ). (F.33) 

Some identities involving the tensor- valued vector harmonics Q ab and associated tensor- valued vector 
harmonics P ab can be derived using the differential identities in appendix |E] and (IF.2 lb - (IF.26b - In each of 
the following identities, the pair (Q ab , Q a ) may be replaced with (P ab , P a )'- 

iV b Q ab = \k-\k 2 -2K)Q a , (F.34) 

e 2 V {a V c Q a)c = -\{k 2 -2K)Q ab , (F.35) 

i 2 V 2 Q ab = -{k 2 -AK)Q ab , (F.36) 

Z 2 ^[aSb]Qcd = \K [(h ac Q bd - h bc Q ad ) + (h ad Q bc - h bd Q ac )] . (F.37) 

F.3 Tensor Harmonics 

Tensor harmonics are defined as solutions of 

V 2 Qa6 = ~ Qab, (F.38) 

where Q ab is divergence-free, trace-free and time independent: 

V b Q ab = 0, Q a a = 0, Q ab = 0. (F.39) 

As with vector harmonics, every tensor harmonic Q ab defines an associated tensor harmonic P ab via 

P ab = £cur\(Q ab ). (F.40) 

It can be shown that P ab satisfies the definition of a tensor harmonic (IF.38b - <IF.39b using (1E.24I ). (IE.29I) 
and (IE.13I) . It then follows immediately from (1F.38I) that the quantities Q ab and P ab satisfy congruence 
orthogonality: 

u a Q ab = 0, u b Q ab = 0, u a P ab = 0, u b P ab = 0. (F.41) 
Upon taking the curl of (IF.40b and applying (IF.38I) and (IE.31b . we obtain Q ab in terms of P ab : 

Qab = (k 2 + ?,K)-Hcur\{P ab ). (F.42) 



F.3. Tensor Harmonics 
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As with scalar and vector harmonics, we are restricted in our choice of k 2 to the eigenvalues deter- 
mined by Lifshitz (1946): 

, o f n > if K = -1,0 

k =n — 3K < ' . (F.43) 

[ n = 3,4,... ifiT = +l 

We can obtain an important identity for and P^ using ( IB. 61 ) and dE.ll) , as follows (the identity for 
P a b is identical, except with Q a b replaced with P a b)'- 

P V [a V 6] Q cd = \K [(h ac Q bd - h bc Q ad ) + (/i ad Q fec - /i& d Qac)] • (F.44) 



The 3-Cotton-York Tensor 



In this appendix we introduce the 3-Cotton-York tensor, which is useful in the analysis of tensor pertur- 
bations. Specifically, we are interested in using the 3-Cotton-York tensor in place of the co-moving shear 
X ab in developing a governing DE for tensor perturbations. It is of interest in this regard due to its relation 

(2) 

with Bardeen's metric potential Hj, . We will further show that the coefficients of the shear and Weyl 
tensor components can be rewritten in terms of Cn^. To the best of the author's knowledge, this analysis 
of tensor perturbations in terms of the 3 -Cotton- York tensor is new. 

In terms of the trace-free spatial Ricci tensor (I1.40I ). the 3-Cotton-York tensor is defined as (see, for 
example, van Elst and Uggla (1997, p2680-l)) 

( 3 )C ab = curl((% a6 ), (G.l) 

and hence is symmetric and tracefree. In the linearized context, we can use ( IB. 101 ) and ( 14.281 ) to obtain an 
alternative first-order expression for ^C ab in terms of the Weyl tensor components, namely 

^C ab = -HH ab + cur\(E ab ). (G.2) 

We note that an important property of this quantity is that it has zero spatial divergence in the linear theory. 
Indeed, upon taking the divergence of (1G.2I ) and applying (14-3 lb . (14.291 ), (14.301 ) and (IE.29I ), we obtain, to 
first order, 

V b ^C ab = 0. (G.3) 
The evolution equation for ^C ab then follows upon taking the time derivative of (IG.2I ) and applying 
63D , (1^241 . €29 and (IeTTTT) . giving 

2K 

®C ab = -3H ^C ab + ^Ha, - V 2 H ab + 2V (a V c ^ )c . (G.4) 
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G.l Evolution Equation for the 3-Cotton-York Tensor 

Instead of choosing the co-moving shear to be the basic quantity for tensor perturbations (as in section 
14.61 ), we can instead choose the basic quantity to be the dimensionless 3-Cotton-York tensor, defined by 

Cab = t 3 (3) C ab . (G.5) 

Upon recalling from dG.31 ) that the 3-Cotton-York tensor has zero divergence, it quickly follows that 
Cab also has zero divergence, i.e. 

V b C ab = 0. (G.6) 
Hence, C a b can be expanded purely in terms of associated tensor harmonics, according to 

k 

The derivatives of C a b along the fundamental congruence are obtained by differentiating (1G.7I ). giving 

C(ab) = E d W P i 2 ' fc) ' C W = J2^) P ab k) - (G-8) 

k k 

Upon differentiating (1G.5I ) and using (13.91) and (1G.4I ) we obtain 



' 2K 

Cab = £ 3 ( " V 2 ) H ab . (G.9) 



Differentiating again and applying (IG.5I ) and ( IE. 171) yields 

C ab + 3HC ab + - V 2 ^) C ab = 0. (G. 10) 

We now expand each term of this equation in terms of tensor harmonics using (1G.7I ) and (IG.8I ). applying 
(tF. 381) to remove the Laplacian term, obtaining 

C\ k) + 3HC\ k) + ^{k 2 + 2K) C {k) =0. (Gil) 

This is then the second-order DE governing the evolution of the 3-Cotton-York tensor. In terms of confor- 
mal time, this DE reads 

C'( k) + 2HlC{ k) + (k 2 + 2K)C (k) = 0, (G. 12) 

(2) 

which is identical to Bardeen's DE for the metric potential Hj, , given in Bardeen (1980) equation (4. 14). 
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G.2 Tensor Perturbations 

We now show the relation between (1G.11I ) and the governing DE for tensor perturbations (14.76b - We 
note that it follows by dG.31 ) that the 3-Cotton-York tensor is purely tensorial, i.e. the scalar and vector 
contributions to ^CVfc) are exactly zero. Expanding ^'Ct k \ in terms of associated tensor harmonics yields 

(3) C a , = E (3) ^)^ 2) - (G.13) 

k 

The coefficients ^Ci k \ can then be determined in terms of and the tensor harmonics via (IG.2I ) upon 
applying (147771 (l4~80l . (1478Tb and (1PT401 We obtain 

"7^2 (*(*) + 2HX (k)) ■ (G.14) 
If we expand dG.91 ) in terms of tensor harmonics using ( 14.771 ), dF.381 ) and dG.81 ), we obtain 

C (k) =t{k 2 + 2K)H {k) , (G.15) 

and hence using (14.8 lb that 



1 

/c 2 + 2K 

Now using (IG.16b . we can rewrite (I4.79I) - (I4. 8 1 1) as 



*W = 1,2 ■ oA ' (G " 16) 



= ~ fc 2 + 2K F (fc) + 2ff ^wJ ' (G - 18) 



H ^ = WTm 1 ' 1 ^- ,G - I9) 



Thus, we obtain the coefficients of the shear and Weyl tensor components in terms of the coefficients of 
the 3-Cotton-York tensor C( k y 



Appendix H 



Bessel and Associated Legendre 

Functions 



In this appendix we present solutions of the Bessel DE, toroidal Legendre DE and conical Legendre DE, 
which are used in solving the governing DEs for scalar and tensor perturbations. In section IhTTI we present 
the normal forms of the Bessel and associated Legendre DEs and their general solutions in terms of Bessel 
and associated Legendre functions. As discussed in section IH.21 these solutions can be written in terms 
of elementary functions in certain special cases. Further, resonant solutions of the Legendre DE appear 
in certain cases, as discussed in section HOI Finally, the asymptotic behaviour of the general solutions is 
given in section ITT741 



H.l Normal Forms of Bessel and Associated Legendre DEs 

In this section we present the normal forms of Bessel's DE, the toroidal Legendre DE and the conical 
Legendre DE along with their general solution in terms of Bessel functions and associated Legendre 
functions. Formally, the toroidal and conical Legendre functions are defined as the composition of the 
associated Legendre functions with cosh and cos, respectively. 



Bessel's DE Normal Form: 



dT 2 



+ 



(P 2 



rp2 



X = 0. 



(H.l) 
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Solution: 

X(T) = T5 [dJ^aT) + c 2 Y fl (aT)} . (H.2) 



Legendre DE (Toroidal) Normal Form: 



d 2 X 
dT 2 



+ 



> + *) 



2 ' sinh 2 T 



1\2 



X = 0. 



(H.3) 



Solution: 



X(T) = (sinhT)2 [ciP^(coshT) + c 2 Q^(coshr)] 



(H.4) 



Legendre DE (Conical) Normal Form: 



d 2 X 

~dr 2 



+ 



2 _ I> 



sin 2 T 



X = 0. 



(H.5) 



Solution: 



X(T) = (sinT)2 [ciP^(cosT) + c 2 Q£(cosT)] 



(H.6) 



H.2 Special Solutions of the Governing DEs 

For certain values of 7, ^ and u, the Bessel functions and associated Legendre P functions reduce to 
closed-form trigonmetric expressions. Further, in each of these cases we can derive a second linearly 
independent solution Q of the associated Legendre DEs dH.31 ) and dH.51 ) that is defined so as to reduce to 
(bz) in the limit z — > 0. This solution is a preferred choice, since it emphasizes the relation between 
solutions for K = ±1 and solutions with K = in the early time limit. In this section, we present 
the first few instances of this behaviour and its importance in modelling perturbations with a given fluid 
background. 

For brevity, we shall make use of the notation 

C(z) = cos(z), S(z) = sm(z), (H.7) 

and use the tilde ^ to denote dependence up to a multiplicative constant. 
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\x = \,v = — i + \/^& Describes evolution of (C, 7 = I). 



z^-Y^bz) - C(6z), (H.8) 

z^-J^bz) - S , (6z), (H.9) 

(sinhz) 1 / 2 -P^coshz) - C(6z), (H.10) 

(sinhz) 1 / 2 • Q^(coshz) - S(6z), (H.ll) 

(sinz) 1 / 2 -P^(cosz) - C(6z), (H.12) 

(sinz) 1 / 2 -Q^(cosz) - (H.13) 

// = !, 1/ = — 1 + v^fr Describes evolution of (C, 7 = 1), (A', 7 = |), (P, 7 = |). 

z l / 2 -Y^{bz) - z^Cibz) + 65(6z), (H.14) 

z 1/2 -J^(bz) - -z _1 5(6z) +6C(6z), (H.15) 

(sinhz) 1 / 2 -P^(coshz) - (cothz)C(bz) + bS(bz), (H.16) 

(sinhz) 1/2 • Q^(coshz) - -(cothz)5(6z) + 6C(6z), (H.17) 

(sin z) 1/2 (cos z) - (cotz)C(te) + (H.18) 

(sinz) 1 / 2 -Q^(cosz) - -(cotz)S(6z) + bC{bz). (H.19) 

//=§, 1/ = — i + Describes evolution of (C, 7 = A), (A\ 7 = 1). 



z 1/2 • y„(6z) 


- (6 2 


- 3z~ 2 )C(&z) - 36z _1 5(6z), 


(H.20) 


^ 1/2 • Mbz) 


- (6 2 


-3z" 2 ),S(6z) + 36z~ 1 C(6z), 


(H.21) 


(sinhz) 1 / 2 ■P^(coshz) 


- (6 2 


- 3 coth 2 2 + l)C{bz) - 36(coth z)S(6z), 


(H.22) 


(sinhz) 1 / 2 -Q^coshz) 


- (b 2 


- 3coth 2 z + l)S{bz) + 36(coth z)C(6z), 


(H.23) 


(sin z) 1/2 (cos z) 


- (6 2 


- 3 cot 2 z - l)C{bz) - 36(cot z)S(bz), 


(H.24) 


(sinz) 1 / 2 -Q£(cos z) 


- (6 2 


- 3 cot 2 z - l)S(bz) + 36(cot z)C{bz). 


(H.25) 
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H.3 Resonant Solutions of the Legendre DE 

For certain special cases, the associated Legendre functions P and Q are no longer linearly independent 
(see Abramowitz and Stegun p333, eq. (8.1.8)). One such example is the case of scalar perturbations with 
a dust background (7 = 1), where a second solution must be determined through alternative means. We 
now give the first five examples of resonance in the Legendre DE and give the second resonant solution 
of the Legendre DE for each case, denoted by Ru{z)- The toroidal solutions are given below. Conical 
solutions are obtained by replacing hyperbolic functions with trigonometric functions. 



(sinh z) 1/2 -P±0 2 (cosh z) - 1, (H.26) 
(sinh z) 1/2 -R]!u 2 (cosh z) - z. (H.27) 



(sinh z) 1/2 • p5i/ 2 ( cosh z ) ~ cothz, (H.28) 
(sinh z) 1/2 -R 3 Jy 2 (cosh z) - zcothz-1. (H.29) 



(sinhz) 1 / 2 -P^Jcoshz) - 2(cothz) 2 + (sinhz)- 2 , (H.30) 
(sinh z) 1/2 • R 5 Jy 2 (cosh z) - z (2(coth zf + (sinh z)- 2 ) - 3(coth z). (H.31) 



" 2' 2' 2 



(sinh z) 1 / 2 -P^ 2 (cosh z) - (sinhz)- 1 , (H.32) 
(sinh z) • R 3 Jy 2 ( cosn z ) ^ z(sinhz) -1 — coshz. (H.33) 



M ~~ 2' ^ ~~ 2' 2 



(sinh z) 1/2 • P V 3 2 /2 (cosh z) - (cothz)(sinhz)- 1 , (H.34) 
(sinh z) 1/2 • ^3/2 ( cosh z ) ~ 3z(cothz)(sinhz)" 1 - ((coshz) 2 + 2)(sinh z)" 1 . (H.35) 



H.4. Asymptotic Behaviour 
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H.4 Asymptotic Behaviour 

We now present the asymptotic behaviour of the Bessel functions and associated Legendre functions, both 
for early times and late times. 



Bessel Functions 



J u (t -► 0) 



u integer 
t u u noninteger 



Y u (t - 0) 



(H.36) 



J u (t -► +oo) ~ t~ 1/2 sin(/ + 0) 



Y u (t -► +oo) ~ t -1 / 2 COS (t + 0) 



(H.37) 



Associated Legendre Functions (i — > 1) 



1 // = o 

(t-l)H/2 ^integer 

(i — 1)~H/ 2 ^ noninteger 

ln(t - 1) m = 

(t _ i)-|H/ 2 h integer 

(i — 1)~M/ 2 i_l noninteger 



(H.38) 



(H.39) 



Toroidal Legendre Functions (t — > oo) 



For the case 



^ = — ^ + z6, & > 0, 
the asymptotic behaviour of the toroidal Legendre functions is given by 



i>£(cosh(t -► oo)) 
Q^(cosh(t -»• oo)) 



exp(-i/2) cos (6i + </>), 
exp(-t/2)sin(6t + 0), 



(H.40) 

(H.41) 
(H.42) 



as i — ► oo. 
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